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ABSTRACT 


( 

I  I 


In  the  past  the  lateral  wave  has  been  investigated  for  the  cai?  of  a 
sharply  bounded  transition  layer  and  a  layer  with  a  linear  velocity  variation. 

The  interpretation  of  the  lateral  wave  given  in  these  instances  cannot  be  extend¬ 
ed  to  an  arbitrary  smooth  layer,  and  in  many  cases,  the  exact  nature  of  the 
lateral  wave  contribution  becomes  unclear.  It  is  our  purpose  to  clarify  these 
matters  and  to  present  the  characteristics  of  lateral  wave  propagation  on  a 
layer  of  arbitrary  variation.  The  models  employed  can  simulate  an  inhomoge¬ 
neous  plasma  having  a  number  density  profile,  n(z),  which  varies  continuously 
between  two  homogeneous  half  spaces.  An  integral  representation  for  the 
scattered  field  in  the  optically  denser  half  space  is  found  and  evaluated 
asymptotically  in  the  high  frequency  limit.  This  asymptotic  evaluation  is  car¬ 
ried  out  in  two  parameter  ranges:  first,  when  the  layer  is  thick  compared  with 
wavelength;  and  second,  when  the  layer  thickness  is  arbitrary  but  the  observa¬ 
tion  point's  distance  along  the  interface  is  large  compared  with  layer  thickness. 

When  the  layer  thickness  is  large  compared  with  wavelength,  the  asymp¬ 
totic  analysis  of  the  scattered  field  shows  that  the  interpretation  of  the  lateral 
wave  depends  markedly  upon  the  gradient  of  n(z)  at  the  junction  with  the  optically 
rarer  homogeneous  half  space,  ft  is  found  that  when  a  finite  gradient  of  n(z) 
exists,  the  conventional  interpretation  of  the  lateral  wave  contribution  is  correct; 
however,  the  lateral  wave  mecftanism  is  different  in  the  case  of  a  zero  gradient. 
For  observation  points  situated  at  a  large  distance  along  a  layer  of  arbitrary 
thickness,  the  asymptotic  expression  for  the  lateral  wave  contribution  has  an 
amplitude  dependence  on  distance  identical  with  that  for  the  lateral  wave  on  an 
abrupt  interface.  In  addition,  the  lateral  wave  expression  reduces  to  the  thick 
layer  result  for  large  layer  thickness  and  it  reduces  to  the  abrupt  interface 
result  for  thin  layer  thickness. 
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INTRODUCTION 

This  report  is  a  study  of  the  fields  reflected  from  diffuse  plasma 
interfaces  with  particular  emphasis  on  lateral  wave  excitation  and  guid¬ 
ing.  It  was  motivated  to  a  large  extent  by  the  general  lack  of  knowledge 
concerning  lateral  wave  behavior  on  diffuse  transitions.  From  an  exam¬ 
ination  of  the  literature  it  became  apparent  that  Nakamura^  was  the  only 
investigator  to  treat  lateral  waves  excited  on  a  diffuse  transition  and,  in 
his  case,  only  a  linear  velocity  profile  was  considered.  The  diffraction 
effect  in  question  is  relevant  for  such  applications  as  the  scattering  of 
waves  by  inhomogeneous  dielectric  or  plasma  ducts  and  the  radiation 
from  antennae  in  the  presence  of  ionospheric  irregularities.  In  the  latter 
case,  the  presence  of  a  magnetic  field  may  introduce  additional  compli¬ 
cations  ;  however,  the  results  obtained  here  for  the  isotropic  problem 
should  provide  a  basis  for  future  <*tudy  of  lateral  wave  effects  when  a 
magnetic  field  is  present. 

The  transition  layers  to  be  considered  have  a  monotonically 
stratified  number  density  profile,  n(z),  which  varies  continuously 
between  two  homogeneous  half  spaces.  The  width  or  average  width  of 
these  layers  is  proportional  to  the  parameter  L;  as  L  becomes  small 
compared  with  wavelength,  the  transition  layers  tend  toward  an  abrupt 
transition.  The  equivalent  dielectric  constant  for  the  medium  is  given  by 

r.(z)  =  1  - (u)  /uu)2  ,  w2  =  n(z)  e2/m€  , 

P  P  ° 

and  is  representative  of  a  cold  electron  plasma  with  a  background  of 
positively  charged  immobile  particles.  In  the  above  formulae  uu  is  the 
applied  frequency,  oj  is  the  plasma  frequency,  e  and  m  are  the  charge 
and  mass  of  an  electron,  respectively,  and  is  the  free  space  permit¬ 
tivity.  The  field  incident  upon  the  layer  is  produced  by  an  electric  line 
current  source  placed  parallel  to  the  transition  in  the  optically  denser 
half  space. 
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Before  proceeding,  a  review  of  the  past  research  concerning 
reflected  fields  from  stratified  transition  layers  will  be  useful.  The 
ma^or  portion  of  this  research  can  be  divided  into  two  sections.  These 
are  :  first,  the  formulation  and  investigation  of  the  reflection  coefficients 
from  transition  layers  ;  second,  the  investigation  of  the  reflected  fields 
due  to  a  point  source  located  above  these  layers.  The  general  problem  of 
relating  the  dielectric  profile  of  the  layer  to  its  reflection  coefficient  has 
been  dealt  with  in  a  number  of  ways.  Brekhovskilh^  has  derived  two 
representations  for  the  reflection  coefficient.  One  of  these  converges 
rapidly  when  the  layer  is  thin  compared  with  wavelength,  while  the  other 
makes  use  of  the  geometrical-optics  approximation  and  converges  rapidly 
when  the  layer  is  thick  compared  with  wavelength.  Since  an  explicit  form 
of  the  reflection  coefficient  from  an  arbitrary  layer  is  difficult  to  obtain  , 
many  investigators  have  studied  particular  profiles.  Hartree^  has 
studied  the  linear  layer,  and  Epstein'  has  devised  a  layer  which  is 
completely  continuous.  He  then  studied  the  reflection  properties  from  it. 
Heading^  has  recently  generalized  Epstein's  results,  and  Wait^  has 
summarized  the  results  of  many  other  studies. 

Although  the  above  reflection  coefficient  formulae  are  useful, 
they  can  only  be  employed  to  calculate  the  reflected  field  when  a  plane 
wave  is  incident  upon  the  layer.  When  the  excitation  is  in  the  form  of 
a  line  source,  a  complete  spectrum  of  plane  waves  is  excited.  The 
_efle<'+ed  field  in  this  case  can  be  represented  as  a  continuous  sum  or 
integral  over  the  plane  wave  spectrum  weighted  by  the  appropriate  ampli¬ 
tude  coefficients  for  each  spectral  component.  These  amplitude  coeffi¬ 
cients  will  depend  on  the  reflection  coefficient  of  the  layer  and,  as  a 
result,  previous  investigations  of  reflection  coefficient  properties 
become  useful.  It  will  be  the  main  purpose  of  this  report  to  approximate 
asymptotically  the  integral  representation  for  the  reflected  field  in  the 
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bigh  reequency  limit.  The  asymptotic  approximations  will  then  be 
interpreted  geometrically  as  rays,  and  the  results  will  be  co  npared 
with  the  classical  theory  of  geometrical-optics.  This  comparison  will 
clearly  distinguish  the  diffracted  rays  from  those  which  are  predicted 
solely  by  calssical  geometrical  optics. 

As  has  been  mentioned  previously,  all  the  layers  considered 
here  depend  on  the  width  parameter  L  As  L  decreases,  the  layer 
approaches  an  abrupt  transition.  Since  this  is  the  case  we  should 
expect  all  of  our  results  to  approach  the  field  reflected  from  an  abrupt 
transition  as  k^L  becomes  small  (k^  is  the  free  space  wavenumber)  . 

The  requirement  that  k^L  he  small  implies  that  the  field  does  not  vary 
appreciably  over  the  layer's  width.  Lateral  waves  excited  on  abrupt 

(7) 

transitions  were  first  observed  and  studied  by  seismologists.  Jeffery 

(8) 

and  Muskart1  '  investigated  the  observed  waves  theoretically  by  evalu- 

(9) 

ating  the  transient  field  from  a  point  source.  Later  Ott  and  Brekhov- 

(10) 

skikh  investigated  the  time  harmonic  point  source  problems.  Other 

investigators  who  contributed  to  the  understanding  of  lateral  waves  on 

(11)  (12)  (13) 

abrupt  interfaces  were  Kruger  ,  Gerjouy  ,  and  Tamir  and  Felsen 

Tamir  and  Felsen  considered  the  lateral  wave  excited  by  a  line  source 

and  they  found  that  the  lateral  wave  had  essentially  the  same  behavior 

as  the  lateral  wave  excited  by  a  point  source. 

All  of  the  above  investigators  found  that  the  reflected  field  from 
an  abrupt  interface  consisted  of  two  contributions  :  a  reflected  and  a 
lateral  wave.  In  Fig.  I  the  ray  interpretation  of  both  of  these  contribu¬ 
tions  is  shown.  We  see  that  a  ray,  emitted  from  the  source  toward  the 
layer,  is  reflected  from  the  interface  at  the  angle  of  incidence  and  then 
proceeds  to  the  observation  point.  The  ray  contribution  for  a  line  source 
excitation  has  an  amplitude  dependence  of  0(k  ^  )  .  The  lateral  wave, 
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on  the  other  hand,  can  be  interpreted  in  terms  of  a  ray  that  is  emitted 
from  the  source  at  the  critical  angle  (8  =  sin  Jt~ )  .  This  ray  is 

refracted  along  the  interface,  and  it  then  sheds  energy  into  the  reflected 
field.  Its  contribution  to  the  reflected  field  for  a  line  source  excitation 
is  given  by 


WVb1 

a  e 
a 


(*  L) 
o  p 


3/2 


where  the  constant  a  is  the  excitation  coefficient,  and  L,  ,  L  ,  L  are 

a  1  p  2 

defined  in  Fig.  L  We  >  '  ould  note  that  the  wave  has  an  algebraic  decay 
with  distance  along  the  interface.  This  decay  results  from  the  continual 
shedding  of  energy  into  the  reflected  field  as  the  lateral  ray  progresses 
along  the  interface. 


Fig.  I 

Reflected  Field  From  an  Abrupt  Interface 
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An  examination  of  the  two  contributions  to  the  reflected  field  chows 

that  the  lateral  wave  contribution  is  smaller  by  a  factor  of  k~ *  than  the 

o 

reflected  wave  contribution  ;  nevertheless,  in  many  physically  meaningful 
situations,  the  lateral  wave  effect  is  of  importance.  For  instance,  when 
the  transient  field  is  observed  at  points  located  a  large  number  of  wave¬ 
lengths  along  the  interface,  the  lateral  wave  contribution  is  seen  as  a 
first  arrival.  For  this  reason  it  is  sometimes  known  as  a  head  wave.  The 
lateral  wave  also  becomes  important  in  the  time  harmonic  case  when  the 
medium  containing  the  source  has  a  slight  amount  of  loss.  The  wave  spends 
most  of  its  time  in  the  lower,  lossless  medium  and,  as  a  result,  becomes 
the  dominant  effect  for  observation  points  located  a  large  distance  along 
the  interface.  In  our  study  of  lateral  waves  we  shall  restrict  ourselves  to 
time  harmonic  problems  ;  however,  the  results  can  be  related  to  transient 
phenomena. 

In  an  attempt  to  understand  the  nature  of  lateral  waves  excited  on 
transition  layers,  the  reflected  field  from  four  profiles  is  investigated 
here.  These  profiles  are  :  linear,  parabolic,  Epstein  and  double  exponen¬ 
tial  ;  they  are  considered  in  Chapters  1  through  4,  respectively.  'The 
profiles  have  not  been  chosen  at  random,  but  instead  are  selected 
because  the  wave  functions  for  the  transition  can  be  related  to  well  known 
functions,  and  because  each  successive  transition  demonstrates  an  aspect 
of  lateral  waves  not  shown  by  the  previous  profiles.  Finally,  in  Chapter 
5,  the  conclusions  drawn  from  the  first  four  chapters  are  extended  to  an 
arbitrary  layer  when  possible  .  The  general  method  of  investigation  is: 
first,  to  formulate  the  integral  representation  for  the  reflected  field  ; 
second,  to  evaluate  the  formal  solution  asymptotically  for  thick  layers 

(k  L»l);  and  third,  to  perform  an  asymptotic  analysis  for  arbitrary 
o 

layer  thickness.  It  has  been  found  that  in  order  to  obtain  an  asymptotic 
estimate  of  the  reflected  field  for  arbitrary  layer  thickness,  it  is 
necessary  to  assume  that  the  observation  point  is  far  from  the  source. 
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CHAPTER  1. 

THE  LINEAR  TRANSITION  LAYER 
li  1  Introduction 

To  begin  oar  investigation  of  lateral  waves  vr«  rrlii  study  the 
reflected  field  from  a  linear  transition.  This  transition  is  composed  of 
a  ramp  dielectric  variation  between  two  homogeneous  dielectric  half  spaces. 
A  source,  placed  in  the  denser  half  space,  illuminates  the  transition  and 
causes  a  reflected  field.  This  reflected  field,  which  will  be  the  main  object 
of  our  investigation,  can  be  represented  by  a  continuous  sum  or  integral 
over  the  reflection  coefficient  of  the  layer  times  the  spectral  components 
excited  by  the  source.  The  formulation  and  some  of  the  properties  of  the 
reflection  coefficient  have  been  investigated  by  Hartree  .  However,  no 
attempt  has  been  made  by  him  to  investigate  the  source  problem. 

The  integral  representation  for  the  reflected  field  is  too  complex  to  be 
directly  integrated  and  therefore  it  must  be  asymptotically  approximated. 
Two  separate  asymptotic  evaluations  of  the  integral  representation  are 
carried  out :  first,  when  the  layer  thickness,  L,  is  large  compared  to  wave¬ 
length,  i.e. ,  k  L»l,  and  second,  when  the  observation  point  is  far  from 
o 

the  source.  In  the  first  case  we  will  obtain  an  asymptotic  approximation  to 
the  reflected  field  for  all  observation  points  for  large  ^L.  The  asymptotic 
approximation  obtained  by  the  second  procedure,  on  the  other  hand,  will  be 
valid  for  observation  points  which  are  far  from  the  source  compared  to  the 
layer  thickness.  In  both  evaluations  of  the  reflected  field,  special  emphasis 
will  be  placed  on  lateral  wave  contributions  and  their  interpretation. 

Before  proceeding  with  our  investigation,  some  of  the  pertinent  work 

(14) 

that  has  been  done  on  similar  problems  should  be  mentioned.  Orlov 
has  found  the  ray  trajectories  which  are  reflected  from  a  ramp  dielectric 
variation  as  shown  in  Fig.  1. 1  when  k  z_  »1.  Here  z_  is  the  width  of 

O  L  L 

th»*  tin  of  the  dielectric  layer  which  supports  propagating  waves. 

’  ..gat  modification,  these  results  can  be  used  to  give  the  structure 

of  the  reflected  ray  trajectories  for  the  linear  layer  when  k^L  is  large. 
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In  the  text  we  have  not  used  Orlov's  results  directly,  but  instead,  have 
rederived  them.  An  investigation  of  the  lateral  wave  which  is  excited  in  a 
medium  with  a  linear  velocity  transition,  has  bc-en  considered  by  Nakamura 
when  the  observation  point  is  far  from  the  souree.  Because  of  the  analogous 
behavior  of  acoustic  and  electromagnetic  waves  in  media  of  the  type  being 
considered  here,  Nakamura's  profile  is  equivalent  to  an  inverse  square 
dielectric  profile  as  shown  in  Fig.  1.2.  /  comparison  between  Nakamura's 
lateral  wave  and  our  results  will  be  made  in  an  attempt  to  determine  some 
of  the  invariants  of  transition  layers. 


1.2 


Formal  Solution 


An  electric  line  current  source  of  amplitude  J  is  placed  in  a  stratified 

dielectric  medium,  e(z),at  x  =  0  ,  z  =  z'  and  parallel  to  the  y  axis,  as  is 

shown  in  Fig.  1.  3.  Under  these  conditions  the  only  field  components 

excited  are  E  ,  H  and  H  .  The  electric  field  obeys  the  inhomogeneous 
Trie)  x 

wave  equation' 


[V2  +  k2  c(z)  ]  E  =  -  iau  J  6(x)  6(z  -  z')  . 
o  y  o 


(1.2.1) 


We  will  assume  that 

J  =  l/(iiuuo) 


(1.2.2) 


to  simplify  Eq.  (1.2.  1).  The  particular  dielectric  variation  to  be  used  is 
given  by 


e(z)  = 


A  z/L+  1 


,  z  >  0 
,  -L<  z  <  0 

,  z  <  -  L 


(1.2.3) 


and  is  shown  in  Fig.  1.  2.  The  dielectric  constant  of  the  homogeneous  region, 
z>0  ,  has  been  chosen  to  be  unity,  but  if  the  dielectric  constant  of  this  region 
is  not  unity,  Eq.  (1.  2. 1)  can  be  scaled  to  produce  an  equivalent  dielectric 


Orlov**  Dielectric  Variation 


Comparison  of  Nakamura's  and  Linear  Dielectric  Variation 


Fig.  1.3 
Source  Location 


Fig.  1.4 

op  Sheet  of  a  Four  Sheeted  Riemann  Surface 
p  -  plane  ' 
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constant  which  is  unity. 

Since  the  medium  ia  uniform  in  x,  we  will  represent  by  means  of 
ita  Fourier  transform  *(z,  p) ,  i.  e. , 

k  +•  ik  px 

Ey  =  IT  J  *(Z'P)  *  °  dp*  (1-2.4) 


By  using  the  above  representation  for  E  ,  Eq.  (1.  2. 1)  reduces  to  the  one- 
dimensional  Green's  function  problem 


+  k2  [c(z)-p2]l  *(s,p)  =  -  6(z-z') 


(1.2.5) 


with  the  boundary  conditions  that  t(z,p)  must  be  an  out-going  wave  as 

Z-*i". 


In  Appendix  A  the  formal  solution  for  the  Green's  function  t(z,  p)  is 
obtained  for  an  arbitrary  layer  variation  in  terms  of  a  pair  of  independent 
homogeneous  solutions  to  the  wave  equation  in  the  layer  region.  For  a 
layer  with  a  linear  variation  such  as  the  one  under  consideration  here, 
two  independent  solutions  to  the  wave  equation  in  the  layer  region  are 
Aj(§^  e*  where  A.  is  the  Airy  function^lnd  S2  =  (kQLM)2/3(Az/L+  1-p2). 
By  using  these  independent  solutions  in  Eq.  (A- 13)  and  in  Eq.  (A- 7)  we  obtain 
*  an  explicit  expression  for  ♦  (z,p) .  We  then  substitute  $  (z,  p)  into  Eq.  (1. 2.  4) 
and  obtain  an  integral  representation  for  the  field.  It  is 


E  =  E  ,  +  E 
y  yf  yr 


(1.2.6) 

(1.2.7) 


(1.2.8) 
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where 


?=-  £ 


(1.2.9) 


with 


t 

b 


)a;(^  .-*"»>+ 5j  .*  5t,/6a_(5j  «-in/3)j  )*<?*  .«"»>♦ 

'A'.k\  e‘in/3)±  ?2  e*  5n/6A.(?2  e‘in/3)J  |a'C? 


2  +in/3  _  iTT/6A  .2  in/3. 

2 C  ,±?2C  W  > 


5j  =  (T/A)1/3  Pl 
2  2 

p,  - 1*,-?2:* 


t  =  k  L  , 
o 


p2=  [i-p2]4  . 


(1.2.10) 


In  the  above  expression  for  the  reflection  coefficient,  T  ,  the  symbol 
A^(z)  =  dA.(z)/dz  .  The  branch  cuts  for  the  square  root  functions  and 
p^  must  be  specified  if  the  integrals  in  Eqs.  (1. 2.  7)  and  (1. 2.  8)  are  to  be 
completely  defined.  This  specification  has  been  made  in  Fig.  1.4  where 
the  top  sheet  of  a  four  sheeted  Riemann  surface  is  shown.  To  clarify  the 
designation  of  the  top  sheet,  the  sign  of  the  real  and  imaginary  parts  of 
and  p^  have  been  given  in  each  quadrant.  On  any  particular  sheet  of 
this  Riemann  surface  the  integral  is  a  single-valued  function,  however, 
the  integration  path  was  chosen  on  the  top  sheet  so  that  the  integral  would 
converge  properly  as  p-*±®  . 

The  expression  for  the  field  has  been  divided  into  two  parts  :  E  ^  and 

E  as  shown  in  Eq.  (1. 2.  6)  .  The  motivation  for  this  i*  that  E  ,  is  the 
yr  yf 

direct  field  from  the  source,  that  is,  it  is  the  field  which  would  exist  if 
the  entire  medium  were  homogeneous  while  E  is  the  reflected  field  from 


I 


tii*  layer.  A*  expected,  the  direct  field  integral,  Eq.  (1.2.  7)  integrates  to 


E  =  -i  H(1) 
yf  4  o 


(  kQ  7x2  +  <z- z')2  ) 


(1.2.11) 


where  H*  1  is  the  Hank  el  function  of  the  first  kind  and  zero  order.  We  see 
o 

that  ^  can  be  interpreted  as  a  cylindrical  wave  propagating  away  from 
the  Source.  The  integral  for  the  reflected  field  is  too  complicated  to  be 
integrated  directly,  so  approximate  techniques  must  be  used. 

Before  proceeding  with  the  approximate  evaluation  of  the  reflected 
field,  Eyr»  we  will  transform  the  integral  from  the  p  to  the  p^  plane. 

The  result  is 


i  ! 


p  r  ik  [p  (z  +  z')  +  px] 

1  O  £ 


(1.2.12) 


where  the  square  roots  p  =  -  p2  and  p^  =  +  p2  must  again  be  defined 

an  afour  sheeted  Riemann  surface.  The  upper  two  sheets  of  this  surface  are 
shown  in  Figs.  1.  5  and  1.  6  along  with  the  transformed  integration  path,  C. 

This  transformation  was  motivated  by  the  fact  that  the  interesting  contribu¬ 
tions  to  E^  come  from  the  region  near  the  branch  point  p  =  */c^  .  By  trans¬ 
forming  the  integral  to  the  p^  plane,  this  branch  point  and  the  branch  point 
at  p  =  ^  are  eliminated.  The  disappearance  of  these  two  branch  points  is 

accompanied  by  the  appearance  of  two  new  branch  points  at  p  °  .  For 

1  1  (13) 

a  more  detailed  explanation  of  the  transformation,  consult  Tamir  and  Felsen'  \ 


Evaluation  of  Reflected  Field  for  k  L»1 

o 


1.3.1.  Geometrical  Optics 


Before  attempting  a  rigorous  asymptotic  evaluation  of  the  reflected 

field  when  k  L»l,  we  shall  investigate  the  rays  which  are  emitted  from 
o 
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ths  source  when  the  layer  is  thick  compared  to  wavelength.  These  rays 
divide  into  three  basic  types  :  direct,  transmitted  and  reluming.  The  direct 
rays  are  emitted  from  the  source  in  the  angular  range  -n/2<3<0  where  6 
is  the  angle  of  emission  of  the  ray  from  the  source  ana  is  shown  in  Fig.  :.7. 
The  returning  and  transmitted  rays  are  emitted  from  the  source  in  the  angu¬ 
lar  ranges  0<9<6c  and  8^  <8<tt/2  respectively,  where  the  angle  6^  is  the 
critical  angle.  This  angle  is  defined  <s  the  angle  whose  corresponding  ray 
‘  (critical  ray)  has  its  turning  point  on  the  lower  interface .  The  three  basic 
types  of  rays  and  the  critical  ray  are  shown  in  Fig.  1.  7. 

An  examination  of  this  figure  shows  that  the  direct  and  transmitted  rays 
are  of  a  much  simpler  character  than  the  reflected  rays.  The  simplifying 
feature  of  these  ray  types  is  that  neighboring  ray  trajectories  do  not  cross 
one  another,  thus  making  the  formation  of  caustics  impossible.  No  such 
statement  can  be  made  about  the  returning  rays  which  shall  now  be  exam¬ 
ined  in  further  detail. 


The  returning  ray  trajectories  for  x>x^  can  be  found  by  integrating 
the  ray  equation 


=  *,+/ 


pd  t 


(.3.1) 


V«(T)-p 


where 


V  J, 


pd  t 


,  p  =  sin8 


z  VC(T)-p 


and 


e(zf)  =  p 


The  ray  parameter  p  is  the  same  as  the  integration  variable  in  Eq.  (1.  2.  8). 

We  will  not  investigate  the  returning  rays  in  the  region  before  they  turn 

(x<x  )  since  the  rays  do  not  cross  one  another  in  this  region.  Upon  inte- 
t 

grating  Eq.  (1. 3. 1)  ,  for  x>xt  we  obtain 
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Fig.  1.7 


Typical  Ray  Types  for  Linear  Layer 
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z  =  z '  + 


P2* 


-  4p 


z>  0 


2  »  - 


z  = 


(  2p  "  2p, 


P2)2  -  p2  fz<0 


where 


z  =  Az/L  ,  x  =  Ax/L  ,  z 7  =  AzVL 


(1.3.2) 


2  2 

and  c.  <  p  <1  .  It  is  interesting  to  note  that  p  =  e.  is  the  critically 

2  1 

reflected  ray.  When  the  ray  parameter  p  <  ,  the  rays  correspond  to 

transmitted  rays  and  are  no  longer  given  by  Eq.  (1.  3.  2)  . 


An  inspection  of  Eq.  (1.  3.  2)  shows  that  it  appears  in  a  normalized  form 

in  terms  Of  barred  coordinates,  i.  e.  ,  x  ,  z  and  z 1  .  As  a  function  of  these 

new  coordinates,  the  ray  trajectories  depend  only  on  the  source  coordinates 

¥'  .however,  the  equations  are  still  only  valid  for  Cj<p  <1  .  Since  we 

would  like  to  investigate  the  nature  of  the  rays  for  arbitrary  e  ,  i.  e.  , 

2  1 

0<Cj  <1  we  will  assume  0<p  <  1  .  This  corresponds  to  a  layer  with  =  0  . 

If  we  are  considering  a  layer  where  e  is  finite,  then  we  just  use  the  results 

2 

of  the  Cj  =0  case  and  eliminate  those  rays  with  0<p  <  . 


The  returning  rays  for  x>x^  cross  over  one  another  in  such  a  way  that 
they  form  a  caustic.  This  caustic  can  be  found  by  solving  the  constraint 

equation , 


0  =  x  -  2p2p' 


0  -  x2  +  z'x(2p3-p)  .  2px 

P2 


>2-' 
>  z 

4 

P, 


,  z  >  0 


,  z  <.  0 


(1.3.3) 


simultaneously  with  the  ray  equation  (1.  3.  2).  The  constraint  equation  can  be 
found  by  taking  the  partial  derivative  of  the  ray  equation  with  respect  to  p . 
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Th«  elimination  of  the  ray  parameter  p  between  Eqa.  (1. 3.  2)  and  (1.  3.  3) 
can  be  accomplished  only  when  z’>0.  The  result  is  quite  complicated  and 
will  not  be  preuented  here. 

As  an  alternative  method  of  finding  the  equation  of  the  caustic  directly 
we  shall  investigate  its  behavior  near  6  =  0  and  0  =  tt/2  .  With  this  information 
and  with  the  location  of  the  focal  points,  we  shall  have  a  fair  understanding 
of  the  caustic's  configuration. 

From  the  -ay  and  constraint  equation  we  find  that  the  point  on  the 
caustic,  corresponding  to  0  =  0  ,  is  located  at  x  =  0,  z  =  -  z^/(l  -z’,/4)  . 

The  slope  of  the  caustic  at  this  point  is  zero,  i.  e.,  dz/dx  =  0.  As  the 
angle  6  -*tt/2  ,  the  parameter  p^  becomes  small  and  the  equations  simplify 
which  allows  us  to  obtain  an  asymptote  to  the  caustic  in  this  region.  It  is 

z  =  -  ““  .  P2  «  °.  (1.3.4) 

x 

To  complete  the  description  of  the  caustic,  the  second  constraint  equation 
must  be  found.  This  equation,  together  with  the  ray  and  first  constraint 
equation,  will  give  the  location  of  the  foci.  The  second  constraint  equation, 
obtained  by  taking  the  partial  derivative  of  Eq.  (1.  3.  3)  with  respect  to  p  ,  is 


0  = 


4P,-1 


0  =  x  - 


.  —  ,2  3 

•*  Z  p 


P2[z'(4p2-  l)-2p2] 


z  £  0 


z  <  0  . 


(1.3.5) 


For  z£0 , 
It  is 


the  location  of  the  focus  can  be  found  explicitly  ina  simple  form. 


x  =  3^3/2  , 


z  = 


(1.3.6) 
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From  these  equations  we  see  that  there  is  one  focus  when  I'<  1/2  while 
no  foci  exist  for  7  >0  when  z*  >  1/2  .  For  7<  0  we  find  that  the  foci  occur 
at  ray  parameters 

2  .  *1  f 5±2  «A-67'  ] 

P2  *  2  L  8Y'  +  3  J  *  (1.3.7) 

When  7' >2/3  ,  the  ray  parameters  p^  are  complex  and  no  foci  occur. 

'  b  the  interval  l/2<7/<2/3  ,  two  real  p^  exist  and  two  foci  occur  for?<0  . 
When  0<7#<l/2  ,  two  real  p^  still  exist  but  one  of  the  ray  parameters  leads 
to  a  7>0  . 

To  summarize  our  investigation  of  the  foci,  we  have  found :  first,  for 
large  7* ,  no  foci  exist;  second,  as  7'  becomes  less  than  2/3,  two  foci 
appear  in  the  layer ;  and  third,  as  zf  becomes  less  than  1/2  ,  one  of  the 
two  foci  in  the  layer  region  moves  into  the  homogeneous  region.  It  is  also 
interesting  to  note  that  at  7'  =  2/3,  two  foci  are  located  at  the  same  spot. 

A  check  shows  that  the  third  constraint  equation  is  zero  at  this  point. 

A  graph  of  the  caustic  has  been  plotted  for  each  of  the  three  cases 
given  above.  These  graphs  appear  in  Figs.  1.8,  1.  9.  and  1.  10  along  with 
the  locus  of  turning  points  and  the  critical  ray.  The  portion  of  the  caustic 
near  p  =  0  is  not  shown  in  these  figures  since  >0  .  To  alleviate  this,  we 
have  presented  a  sketch  of  the  caustic  for  =0  and  z'  >2/3  in  Fig.  1.  11. 

The  behavior  of  the  portion  of  the  caustic,  corresponding  to  rays  with  small 
8  ,  remains  essentially  the  same  when  7 '  <2/3  . 


1.3.2  Asymptotic  Evaluation 

As  mentioned  in  the  previous  section,  the  integral  representation  for  the 
reflected  field,  appearing  in  Eq.  (1.  2.  8),  is  too  complicated  to  be  integrated 
directly.  To  effect  its  evaluation  we  will  expand  the  reflection  coefficient  in 
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Az 


a  geometric  series  and  then  make  use  of  asymptotic  techniques  to  integrate 
each  term  of  the  series  individually.  This  will  result  in  a  representation  for 
the  reflected  field,  E  ,  which  can  be  interpreted  in  terms  of  geometrical- 
optic  ray  contributions. 

The  reflection  coefficient  given  by  Eq.  (1.  2.  9)  can  be  written  as 


r 

oa 


rr  (r  .  -  r 

L  ob  oc 

i  -  rT  r 

L  o. 


) 


(1.3.8) 
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where 


oa 


,  2  in/3  _  1tt/6  ..2  irr/3 
A-  (?2  e  >  +  ?2e  Ai<?2*  * 

A  /  ..2  iTT/3.  I  iTV6.  ..2  iTT/6. 

Ai<V  )  -  52«  A.(?2«  ) 


rob  “ 


a/._2  -in/3.  i 5n/6  ..2  -in/3. 

Ai  (?2 e  >  -?2®  Ai(?2*  * 

.  /  #-2  +in/3  '  in/6  2  +in/3. 

A.(?2e  )-?2e  A.(?2.  ) 


A//F2  -in/3.  i  5tt/6  2  -irr/3. 

Ai  (?2  ®  )  +  ^2  6  Ai(?2  6  • 

^oc  A  /  .»2  +in/3  i  5tt/6  ._2  +in/3. 

Ai  (?2 e  )  +  ?2«  a-(?2  e  ) 


,  2  in/3  _  in/6  2  in/3 

.  At(V  >  +  ;ie  Ai(5ie  1 

L  ~  a//p2  -in/3  i5n/6  2  -in/3 

Ai<?l  e  )Mje  A.fSj  e  ) 


(1.3.9) 


(1.3.10) 


(1.3. 11) 


(1.3.  12) 


The  denominator  of  the  second  term  in  Eo.  (1.  3.  8)  can  be  expanded  in  a 
geometric  series  if  Ir^T  |<  1  •  ^is  condition  must  be  met  at  all  points 
along  the  integration  path  C.  An  evaluation  of  shows  that  its  magnitude 

is  less  than  one,  except  for  the  portion  of  the  path  Re  p^  =  0  ,  •/K<lxnpl<a>  . 
Along  this  portion  of  the  path  |r  1^  |  =  1  ,  however,  we  can  show  that  a  slight 
deformation  of  the  integration  path  C  to  the  right  of  the  Impj  axis  in  this 
region  makes  | |  <  1  .  The  use  of  this  series  expansion  in  the  reflection 
coefficient  gives  us 


00 

r  =  r  +  Y  r 

oa  L  n  n 
n  =  0 


(1.3.13) 


where 


r  =  r  (r  -  r  )  r®  rn+i 

n  oa  oa  oc  ob  L 


(1.3.  14) 
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We  now  use  the  series  expansion  forT  in  Eq.  (1.  2.  8)  and  interchange  the 
intergration  and  summation  signs.  We  obtain 


E  =  E  +  Y  E 
*r  ro  nr  0  n 


(1.3.  15) 


where 


ro 


Piro.  ikotp2,2+Z')  4  px] 


c  ppa 


dP, 


(1.3.  16) 


and 


i  r  pirn  ik0tP2(^')  +  P*] 

4TT  J  ""  C  dP> 


c  pp* 


(1.3.17) 


The  integrals  E^q  and  E^  will  now  be  evaluated  by  the  method  of  steepest 
descents  when  is  large.  Before  doing  this  however,  we  will  assume 
koL»l  and  replace  1^  and  by  their  asymptotic  approximations.  These 
approximations  are  obtained  by  using  the  asymptotic  expansions  for  the  Airy 
function  and  its  derivative.  The  expansion  for  A^(z)  is  given  by 

GO 

Ai(z)  ~  Z"1/4e’C  I  (_1)kckC’k  »  Urg*| <TT  (1.3.18) 

n=  0 

where 

c  =  2z3/2/3  ,  C  a  1  ,  Cj  *  5/72 

while  A^(z)  can  be  obtained  by  formally  differentiating  Eq.  (1.  3.  18)  .  The 

(171 

additional  coefficients  can  be  found  in  Abramawitz  When  the  argz 

is  close  to  ±rr  Eq.  (1.  3. 18)  is  no  longer  valid  and  another  asymptotic 

(18) 

expansion  including  this  sector  of  argz  must  be  used  *  .  The  expressions 

for  T  and  T  contain  the  Airy  functions  with  four  different  arguments, 
oa  n 


The  locus  of  possible  values,  when  the  argument  is  constrained  to 
have  an  angle  of  ±tt  ,  have  been  drawn  in  Fig.  1.12  for  the  four  '  iry 
functions.  This  figure  serves  to  make  the  region  of  validity  of  Eq.  (1.  3. 18) 
in  the  p^  plane  clear.  The  curves  shown  in  Fig.  1.12  are  the  same  on  both 
the  top  and  on  the  second  sheet  of  the  Riemann  surface  since  the  Airy 
functions  are  even  functions  of  p  . 


The  asymptotic  approximations  for  and  on  the  integration  path 


between  p^  =  and  p^  =  iJK  are 
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(1.3.19) 


and 


r  ~ 
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a  r  ”+1 

n  L 

#1  T  ^  ® 

<koL)  P2 


.n+1 


where  An=(-1)  (15A  C1/2)  . 

and  Pj  =  i t/K  ,  making  sure  that  p 
then  Eq. (1.  3.  20)  reduces  to 


i4(n+l)d/3  +  inTT/2 

e  (1.3.20) 


K  we  restrict  p^  to  lie  between  the  origin 
does  not  come  too  close  to  the  origin. 


r 

n 


A 


_ n 

"'oL,P2n 


i4(n+l)5*/3  +  i(n-2)n/6 

e 


(1.3.21) 


while  if  we  restrict  p^  to  the  real  axis  (0<Repj  <\Zi^)  and  again  keeping 
p^  away  from  the  origin,  Eq.  (1. 3.  20)  reduces  to 


n 


,  .2n+l  3n 

<koL)  P2 


3(n+l) 

pl 


i4(n+l)(^-?i  )/3  +  i(n-2)n/6 


r  ~ 

n 


e 


•  (1.3.22) 
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Fig.  1.12 

Region*  of  Validity  for  Asymptotic  Expansion*  of  Airy  Functions 


Fig.  1.14 

Deformed  Path  for  I  (Second  Sheet) 


Fig.  1.13 

Deformed  Path  for  If  (Top  Sheet) 
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1.22 


By  using  the  asymptotic  approximation  for  1^  in  the  integral  E^, 
the  saddle  point  equation  can  be  determined  by  taking  the  derivative  with 
respect  to  of  the  argument  of  the  exponential  term  in  the  integrand. 
The  result  is 


-  top  sheet 
-I-  second  sheet  „ 


(1.3.23) 


There  are  two  relevant  solutions  to  this  equation  on  the  top' sheet  which 
are  denoted  by  D  and  E  in  Fig.  1.13;  there  are  no  relevant  solutions  on 
the  second  sheet.  The  saddle  point  D  can  occur  anywhere  on  the  integra¬ 
tion  path  between  p^  =*/e|  and  p^  =  iJE  while  the  saddle  point  E  is  always 
located  at  p^  =  0 .  If  we  restrict  the  saddle  point  D  to  lie  on  the  imaginary 
axis  and  deform  the  integration  path  C  into  the  decay  regions,  (Figs.  1. 13 
and  1.  14),  the  field  E^is  asymptotically  approximated  by  the  two  saddle 
point  contributions  D  and  E.  The  contribution  of  E  is  exponentially  small 
however,  since  the  integrand  is  an  odd  function  of  p^  integrated  over  a 
symmetric  interval.  The  contribution  from  D  can  be  interpreted  as  a  ray 
reflected  from  the  z  =  0  interface  as  shown  in  Fig.  1.  19.  The  asymptotic 
order  of  contributions  of  D  and  E  are  shown  in  Table  1.1.  If  now  we 
restrict  the  saddle  point  D  to  the  real  axis,  the  saddle  point  E  is  not  inter¬ 
cepted  when  the  integration  path  C  is  deformed  and  the  total  contribution 
to  the  integral  comes  from  D.  The  interpretation  of  D  is  as  before.  We 
should  mention  at  this  point  that  the  assumption  has  been  made  that  no 
isolated  pole  singularities  exist  between  the  original  and  deformed  paths 
which  contribute  to  the  asymptotic  evaluation  of  E^ .  This  same  assumption 
will  apply  to  the  other  path  deformations  which  shall  occur  in  this  chapter. 
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Type 

Asymptotic  Order 

A 

oo.-<B+,/2>) 

n 

o 

B 

oo>-«2n+J/2>) 

n 

o 

C 

Oft-'-7'6', 

n 

o 

-3/2 

D 

0{k  ) 

o 

-ak 

E 

0(e  °),a>o 

Table  1.  1 

Asymptotic  Order  of  Geometrical 
Optic  Contributions 


The  asymptotic  evaluation  of  the  integral  Eq  can  be  treated  in  a 
similar  way  to  E  q  .  The  asymptotic  approximations  for  given  in 
Eqs.  (1.  3.  20) ,  (1.  3.  21)  and  (1.  3.  22)  are  substituted  ir.  the  integrand  of 
E  and  the  saddle  point  equations  are  then  obtained.  They  are  : 


/  4-  '  4(n+l)p  L. 

,  +£  +  - - - £-)=0,  Re p  =  0  ,  0<frnPl  <JE  (1.3.24) 

1  \  P  A  '  1  1 


x  ,  ♦<"+‘><P2-P1>L 
>1  +  D  +  4 


)  «  0  ,  0<Rep^<y€y,  Bnpj.Q 


(1.  3.  25) 


and 


<  1 


(1.3.  26) 
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where  the  upper  sign  in  Eqi.  (1.  3.  24)  and  (1.  3.  25)  refers  to  the  top 
sheet  and  the  lower  sign  to  the  bottom.  Equation  (1.  3.  26)  is  valid  on 
both  sheets. 

Upon  obtaining  the  saddle  point  locations  from  the  above  equations, 
we  deform  the  contour  C  through  these  saddle  points  and  into  the  decay 
regions.  Since  the  contribution  due  to  the  portion  of  the  integral  lying  in 
the  decay  region  is  exponentially  small,  the  integral  E^  will  be  asymptotic¬ 
ally  equal  to  the  sum  of  the  saddle  point  contributions.  The  contributions 
arising  from  saddle  points  obeying  Eq.  (1.  3.  24) ,  (1.  3.  25)  and  (1.  3.  26)  will 

be  known  as  type  A  ,  B  and  C  contributions  respectively.  The  subscript 
n  n  n 

n  indicates  the  integral  in  which  the  saddle  points  occur.  It  can  be  shown 
that  no  relevant  saddle  points  occur  on  the  second  sheet  for  any  n  and  there¬ 
fore  the  deformed  path  will  be  the  same  as  that  shown  in  Fig.  1.  14  . 

At  this  point  we  shall  focus  our  attention  on  the  relevant  saddle  points 

on  the  top  sheet,  their  contributions  to  the  reflected  field  and  the  geometrical 

interpretation  of  these  contributions.  We  shall  first  consider  the  integral 

E  since  E  contributes  the  dominant  terms  to  the  reflected  field  for  k  L 
o  o  o 

large  and  then,  at  a  later  time,  we  shall  consider  the  integrals  Eq,  n=  1, 2,  3, . . 

Saddle  point  contributions  of  type  A^  will  be  considered  first.  A 
comparison  of  Eq.  (1.  3.  24)  with  the  ray  equation  for  returning  rays  (*  > 0)  , 

Eq.  (1.  3.  2)  ,  shows  that  the  two  are  identical.  This  means  that  saddle 
point  contributions  in  the  interval  Re  p^  =  0  ,  0<hnp^  <  JE  correspond  to 
returning  ray  contributions.  At  this  point  we  can  use  our  knowledge  of  the 
returning  ray  trajectories  to  find  the  location  and  number  of  saddle  points 
in  the  interval  Re  p^=  0  ,  0<lmpj<v^.  Before  proceeding  however,  we 
shall  divide  the  portion  of  the  x-z  plane  with  z  >0  in  two  parts  separated  by 
the  critical  ray  as  shown  in  Fig.  1.19.  To  the  right  of  the  critical  ray 
(region  2)  there  is  either  one  or  there  are  three  returning  rays  passing 


-30- 


througn  an  observation  point.  If  the  observation  point  is  located  inside 
the  cusp  region  there  are  three  rays  while  outside,  there  is  only  one . 

To  the  left  of  the  critical  ray  (region  1)  there  are  either  two  saddle  points 
in  the  interval  Re p ^  =  0  ,  0 < Imp,  <  JK or  there  are  none.  The  two  saddle 
points  occur  when  the  cusp  region  lies  to  the  left  of  the  critical  ray,  and  the 
observation  point  is  located  inside  the  cusp  region.  The  saodle  points  are 
shown  for  the  four  cases,  considered  above  in  Figs.  1. 15  -  1. 18  ,  where 
they  are  denoted  by  A  (n=0)  .  The  contributions  to  the  reflected  fi«  d  at 
these  saddle  points  are  of  0(kQ  '  )  as  would  be  expected  since  they  repre¬ 
sent  returning  ray  contributions. 

Saddle  point  contributions  of  type  will  be  considered  next.  An 
examination  of  Eq.  (1.  3.  25)  shows  that  one  saddle  point  exists  in  th«*  interval 
0<  Re p^  <  •  hn  p^=  0  when  the  observation  point  is  in  region  1  hile  no 

saddle  points  exist  in  the  interval  when  the  observation  point  is  in  r  gion  2. 
The  saddle  point  contributing  to  region  1  can  be  interpreted  as  a  ray  reflected 
from  the  interface  at  z  =  -  L  .  A  ray  of  this  type  is  shown  in  Fig.  1.19.  Its 
asymptotic  contribution  to  the  reflected  field  is  of  which  is  lAc^ 

lower  than  the  returning  ray  contribution.  This  is  expected,  howev<  *  since 
the  ray  under  consideration  is  reflected  at  an  interface  when  the  die-*  ctric 
profile  has  a  discontinuous  first  derivative. 

Finally,  the  saddle  point  contribution  of  type  Cq  will  be  conside  d. 

An  examination  of  Eq.  (1.  3.  26)  shows  that  only  one  saddle  point  existt  and 
its  location  (p^=  0)  is  independent  of  the  location  of  the  observation  p  >int. 

An  analysis  of  the  decay  regions,  the  results  of  which  are  presented  in 
Fig.  1. 15  -  1.  18, has  shown  that  the  saddle  point  is  only  intercepted  when 
the  observation  point  is  located  in  region  2. 

We  will  now  attempt  to  asymptotically  evaluate  Eq  .  The  integral  Eq 
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CRITICALLY  REFLECTED  RAY 


f: j.  1.19 

Ray  Contributions  from  I  and  I 

r  o 


/ 

/ 


Fig.  1.20 

Ray  Contributions  from  I 

c* 
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in  the  vicinity  of  the  origin  doe*  not  appear  to  be  a  simple  steepest 
descent  integral  since  T  is  a  complicated  function  of  k  .  The  situation 

O  |  lA 

can  be  clarified  by  making  the  change  of  variables  =  (t/A)  p^  in  the 
integral.  The  portion  of  the  integral  in  the  vicinity  of  the  origin, 
becomes 


Vl 


ikoCp2(i  +  «')+px+|i  P\l 

e  d?, 


(1.3.27) 


where 


and  C'  is  the  transformed  path  C  in  the  vicinity  of  the  origin.  Now  T  is  no 
longer  a  function  of  k  but  p  and  p  are.  Because  of  the  simpler  dependence 

O  b 

of  p  and  p  upon  k  ,  it  becomes  possible  to  show  that  the  standard  steepest 
L  °(1 9) 

descent  techniques  can  be  applied  to  Eq.  (1. 3.  27)  .  Upon  the  application 
of  these  techniques  we  find 


ik  U/T(z  +  z')  +  +  A>JK L/3 ]  +  i  7TT/12 

A  e  °  1 

E  ~  _L _ 

L  (L  +  VTla/a)3/2 

p  1 


(1.3.28) 


where 


3*r 2(1/J)  e[/4L1/3 


L  =  x  -  >/t~TK  (z  +  z') 
p  1 


(1.3.29) 


and  T(z)  is  the  gamma  function  of  argument  z  .  This  result  can  be  inter¬ 
preted  as  the  contribution  from  a  ray  which  is  excited  by  the  critical  ray 
at  the  z  =  -  L  interface.  This  newly  formed  ray  travels  along  the  lower 

interface  for  a  distance  L  +4*/eT  L/^and  sheds  energy  into  the  upper  region. 

P  1 
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A  typical  ray  trajectory  is  shown  in  Fig.  1.  19.  The  lateral  wave  exicted 

on  an  abrupt  transition  has  a  similar  ray  interpretation  and  amplitude 

dependence  with  distance  along  the  interface.  The  two  lateral  waves. 

however,  have  different  excitation  coefficients.  The  excitation  coefficient, 

,  for  a  lateral  wave  on  a  thick  linear  transition,  depends  on  while 

the  excitation  coefficient  for  a  lateral  wave  on  an  abrupt  transition  has  a 
-3/2 

dependence  of  k^ 

Now  we  will  not  require  n  to  be  zero  and  we  will  investigate  the  asymp¬ 
totic  contributions  to  the  integral  .  To  do  this  we  must  locate  the  relevant 
saddle  points  and  then  find  their  contributions  to  the  reflected  field.  As  men¬ 
tioned  earlier,  the  general  saddle  point  equations  are  given  by  Eqs.  (1. 3.  24)  , 
(1,  3.  291  and  (1.  3.  26)  .  By  replacing  (n+1)  L  by  L  in  the  first  two  of  these 
saddle  point  equations  we  see  that  the  equations  are  the  same  as  the  n=  0  case 
with  L  replaced  by  X.  Since  the  n  =  0  case  was  done  for  arbitrary  L  ,  the 
location  of  saddle  points  of  type  Aq  and  Bn  has  already  been  investigated. 
However,  there  is  one  basic  change.  The  demarcation  between  region  #  1  and 
region  #  2  is  no  longer  the  critical  ray  reflected  from  a  layer  of  thickness  L 
but  is  rather  the  critical  ray  reflected  from  a  layer  of  thickness  of  L.  '’’his 
cirtical  ray  is  shown  for  n  =  2  in  Fig.  1.  20  .  The  third  saddle  point  equation 
given  in  Eq.  (1.  3.  26)  has  one  saddle  point  at  p^=  0  for  any  n  .  From  the 
above  argument  we  conclude  that  the  basic  structure  of  the  deformed  paths 
are  again  given  by  Figs.  1. 15  -  1.  18  . 


The  method  of  steepest  descents  can  then  be  used  to  evaluate  the 

asymptotic  contributions  at  each  saddle  point.  Those  contributions  of  type 

A^  correspond  to  returning  rays  reflected  n  times  from  the  z  =  0  interface 

while  those  of  type  B^  correspond  to  rays  reflected  n+1  times  from  the 

interface  at  z  -  -  L  .  The  C  are  lateral  wave  contributions.  These  are 

n 

excited  by  the  critical  ray  which  is  reflected  n  times  frem  the  z  =  0  interface. 
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The  ray  trajectories  for  these  contributions  are  shown  in  Fig.  1.  20  for 
n  =  2  .  The  order  of  each  type  of  contribution  is  given  in  Table  1.  We  see 
from  this  table  that  as  n  increases,  the  order  of  the  asymptotic  contribution 
decreases.  This  results  from  the  fact  that  the  order  of  a  ray  contribution  is 
reduced  by  0(l/kQ)  each  time  the  ray  is  reflected  from  the  z  =  0  or  z  =  -  L 
interface. 

To  summarize  briefly  we  have  found  that  the  reflected  field  from  a 

gradual  linear  layer  (k^L»  1)  can  be  decomposed  into  a  geometric  -  optic 

series.  The  dominant  terms  of  this  series  in  region  #  2  are  the  returning 

rays  which  are  of  0(k  ,  and  a  lateral  wave  contribution  which  is  of 

-7/6  ® 

0(kQ  )  .  This  lateral  wave  has  a  great  similarity  to  the  lateral  wave 
occurring  on  an  abrupt  transition  except  for  a  difference  in  its  excitation 
coefficient. 

1.4  Evaluation  of  Reflected  Field  fo?  Large  kQLp 

In  this  section  we  shall  explore  the  connection  between  the  lateral  wave 
contribution  observed  on  a  thick  linear  transition  (section  1. 3)  and  the 
lateral  wave  excited  on  an  abrupt  interface.  The  linear  layer  provides  an 
excellent  opportunity  to  do  this,  since  for  small  k^L  ,  the  layer  appears  to 
be  an  abrupt  transition,  while  for  large  k^L  the  layer  thickness  is  large 
compared  with  wavelength.  It  can  be  shown,  in  fact,  that  the  reflected  coe¬ 
fficient,  Eq.  (1.  2.  9)  ,  reduces  to  the  reflection  coefficient  for  an  abrupt 

transition  as  k  L-*0  . 

o 

Our  investigation  of  the  lateral  wave's  character  shall  be  carried  out 

by  performing  an  asymptotic  evaluation  of  the  reflected  field  for  observation 

points  which  are  far  from  the  source  compared  with  layer  thickness.  To  be 

more  specific  we  shall  assume  that  k  L  »1  and  L  »L  .  This  says  that  we 

o  p  p 
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shall  consider  those  observation  points  which  are  a  large  number  of 
wavelengths  and  a  large  number  of  layer  thicknesses  to  the  right  of  the 
reflected  ray,  («  +  *#)  . 

The  integral  representation  for  the  reflected  field  is  given  by 

[p2(s  +  s/)+  par] 

dPx  (1.4.1) 


y* 


-  —  r 

4TT  J 


ik 


PP, 


where 


with 


(1.4.2) 


The  representation  given  in  Eq.  (1. 4. 1)  is  the  same  as  that  presented  in 

Eq.  (1.  2. 12)  except  that  a  different  form  of  the  reflection  coefficient  has 

been  used.  The  alternate  form  of  this  reflection  coefficient  was  obtained  by 

2  2 

choosing  (^(s)  and  *2(*)  in  Appendix  A  as  BJ-?  )  and  A.(-5  )  respectively. 

These  two  functions  are  independent  of  one  another  and  thus  satisfy  the 

requirement  for  choosing  t^s)  and  •.,(*)  •  The  B.(-5  )  not  encountered 

previously,  is  the  Airy  function  of  the  second  kind.  The  reason  for  using  an 

alternative  form  of  reflection  coefficient  is  that  it  will  allow  us  to  put  our 

* 

results  in  a  form  which  is  more  suitable  for  computation. 


The  advantage  of  the  original  form  of  the  reflection  coefficient  was  that 
it  contained  Airy  functions  whose  asymptotic  expansions  were  of  a  simpler 
form  in  the  vicinity  of  the  integration  path  . 
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j£  we  now  use  the  assumption  that  k  L«k  L  then  the  reflectiSB 

o  o  p 

coefficient  in  Eq.  (1. 4. 1)  is  slowly  varying  compared  to  the  exponential 
in  the  integrand,  and  the  standard  steepest  descent  techniques  can  be  u^ed 
to  asymptotically  evaluate  the  integral  representation.  The  saddle  points 
are  found  by  taking  the  derivative  with  respect  to  p^  of  the  argument  of  the 
integrand's  exponential  and  setting  it  equal  to  zero.  The  result  is 


z  + 
P 


-  top  sheet 
+  second  sheet  t 


(1.4.3) 


This  Is  the  same  saddle  point  equation  encountered  in  the  asymptotic 
evaluation  of  E^q  .  The  relevant  saddle  points  which  occur  on  the  top 
sheet  of  the  Riemann  surface  are  located  at 


(1.4.4) 


where  we  have  assumed  that  x> \/e^7A  (z+z')  or  L  >0.  The  original 
integration  path  is  deformed  through  the  two  saddle  points  given  above 
and  into  the  decay  regions.  We  will  assume  that  there  are  no  singular¬ 
ities  of  importance  lying  between  the  two  paths.  The  decay  regions  and  the 
deformed  path  are  shown  in  Figs.  1.13  and  1. 14  . 


We  see  now  that  when  k  L  is  large,  the  reflected  field  is  composed 

of  two  contributions.  The  first,  of  these  occurring  at  p^  =  i  Jx  A  -(z+z*) 

is  a  reflected  ray  from  the  interface  at  z  =  0  .  Whon  l«k  L«k  L  this 
1  o  o  p 

reflected  ray  corresponds  to  a  returning  ray  that  has  its  turning  point 

close  to  the  z  =  0  interface.  When  k  L  is  small  the  contribution  is  the  same 

o 

as  a  reflected  ray  from  an  abrupt  interface.  The  contribution  from  the 

-1/2 

saddle  point  is  of  0(k  )  as  would  be  expected. 


The  second  contribution  to  the  reflected  field  comes  from  the  saddle 
joint  at  Pj  =  0  and  is  the  lateral  wave  contribution.  This  contribution  was 
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In  the  above  we  have  denoted  the  lateral  wave  contxiubtion  by  E  as  was 

1 -i 

dona  in  the  previous  section^Eq.  (1.  3.  27)  .  The  amplitude,  A(o)  and  phase 
f(o)  functions  have  been  plotted  by  computer  and  are  Bhown  in  Figs.  1. 21 
and  1. 22  .  Those  figures  also  show  the  approximate  lateral  wave  amplitude 
and  phase  when  the  layer  is  thick  compared  with  wavelength. 

The  physical  interpretation  of  the  lateral  wave  contribution  can  best 
be  gotten  from  examination  of  small  and  large  limits.  When  t  is  small  or 
when  the  layer  appears  to  be  abrupt,  we  find  A(c)-*1  and  t(cr)‘*0  .  By  using 

these  limits  in  Eq.  (1.  4.  5),  we  find  that  E  reduces  to  the  lateral  wave 

Li 

contribution  on  an  abrupt  interface.  On  the  other  hand,  when  l«k  L«k  L 

o  o  p 

then  E^  reduces  to 
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Fig.  1.22  Lateral  Wue  Phaie  vs.  Normalised  Transition  Thickness 
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ik  [^"(*+*/)  +  x  +  AjK  L/3]  +  17TT/12 

el  ■  ~W  <1-4-10' 

L 

p 

where  is  given  by  Eq.  (1.  3.  29).  When  this  result  is  compared  with 

Eq.  (1.  3.  28)  and  use  is  made  of  the  fact  that  k  L«k  L  ,  we  find  that  the 

o  op 

two  expressions  are  identical. 

As  was  stated  in  the  introduction,  Nakamura  has  obtained  a  uniform 

expression  for  the  lateral  wave  excited  on  an  inverse  square  dielectric 

profile  which  is  valid  under  the  same  conditons  assumed  in  this  section. 

A  comparison  between  our  resultyEq.  (1.  4.  5)^and  Nakamura's  shows  that 

a  lateral  wave  is  excited  for  all  layer  thicknesses  in  both  cases  and  both 

-3/2 

have  amplitude  dependences  of  L  .  However,  the  excitation  coeffi¬ 
cient  as  a  function  of  layer  thickness  appears  to  depend  on  the  detailed 
behavior  of  the  wave  functions  in  each  particular  medium,  and  nothing  in 
general  can  be  said.  When  the  layer  thickness  becomes  large  compared 

with  wavelength,  the  excitation  coefficient  simplifies  and  both  results  have 

-7/6 

a  wavenumber  dependence  which  is  0(k  )  . 
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CHAPTER  2. 

THE  PARABOLIC  TRANSITION  LAYER 


2. 1  Introduction 

In  Chapter  1  a  detailed  study  of  the  reflected  field  from  a  linear 
transition  layer  was  made.  There,  it  was  found  that  a  lateral  wave  was 
excited  on  the  layer  for  all  layer  thicknesses.  In  this  chapter,  we  will 
continue  our  study  of  lateral  waves  by  considering  the  reflected  field  from 
a  parabolic  transition.  Other  parameters  influencing  the  reflected  field, 
such  as  source  location  and  orientation,  will  remain  the  same  as  in  Chapter 
1.  The  dielectric  profile  to  be  considered  is  given  by 


*(*)  = 


A(z+L)2/L2  +  e. 


z  >0 
-L< z <0 
z<  -  L 


A  =  l  -  e. 


(2.1.  1) 


and  is  shown  in  Fig.  2.  1  .  The  basic  similarities  between  the  parabolic 
and  linear  transitions  are  their  thickness  L  ,  height  A  and  finite  slope  at 
z  =  0  .  Their  fundamental  difference,  on  the  other  hand,  is  the  slope  of  the 
profiles  at  z  =  -  L  .  The  linear  layer's  slope  at  that  point  is  finite  while  the 
parabolic  profile  has  zero  slope  there. 

Our  investigation  of  the  parabolic  transition  has  been  motivated  by  the 
behavior  of  the  critical  ray  trajectory  when  the  layer  is  thick  compared  to 

-w’ 

wavelength  .  This  critical  ray,  unlike  the  critical  ray  in  a  linear  transition, 
never  becomes  tangent  to  the  lower  interface  but  only  approaches  it  asymptot¬ 
ically  as  x-*00  .  In  the  previous  chapter  the  lateral  wave  contribution  has 
been  interpreted  as  a  wave  excited  at  the  point  that  the  critical  ray  becomes 
tangent  to  the  lower  interface.  If  this  supposition  is  true  for  the  parabolic 
transition,  no  lateral  wave  will  be  excited  when  the  parabolic  transition 
thickness  is  large  compared  to  wavelength.  In  an  attempt  to  clarify  the  above 


9 


TRANSMITTED  RAY 
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question  we  shell  asymptotically  approximate  the  formal  solution  first  for 

large  k  L  and  then  for  large  k  L  with  L  »L.  Here  L  is  the  same  as 
o  o  p  p  p 

that  defined  in  Chapter  1  .  When  kQL»l  or  the  layer  thickness  is  large 
compared  to  wavelength,  a  knowledge  of  reflected  ray  trajectories  will  be 
necessary  to  carry  out  an  asymptotic  evaluation  of  the  integral  representation. 
These  ray  trajectories  have  been  considered  in  some  detail  by  Orlov'  ' 
for  z>0  ;  his  results  will  be  used  in  the  text  . 

2.  2.  Formal  Solution 

The  integral  representation  for  the  field,  due  to  a  line  source  located 
above  a  parabolic  transition,  is  obtained  in  an  analogous  way  to  the  integral 
representation  found  in  Chapter  1  .  The  only  component  of  the  electric  field 
which  is  excited  is  E^  which  obeys  the  inhomogeneous  wave  equation 

+  k^e(z)]  Ey  =  -  iuuMQ  J  6(x)  6(z-z;)  (2.  2.  1) 

where  e(z)  is  given  by  Eq.  (2.  1.  1) . 

The  application  of  Fourier  integral  techniques  to  the  above  equation 

reduces  it  to  a  one-dimensional  Green's  function  problem.  The  formal 

solution  to  this  Green's  problem  has  been  presented  in  Appendix  A  in 

terms  of  two  independent  solutions  to  the  wave  equation  in  the  layer 

region.  For  the  parabolic  transition  which  is  being  considered  in  this 

chapter,  two  independent  solutions  to  the  wave  equation  in  the  layer  region 
2  *  2 

are  E(-a^  ,  ?)  and  its  conjugate,  E  (-a^  ,  §)  .  These  two  functions  which 

are  defined  by  Abramowitz^^have  an  order  -a^  and  an  argument  §  given  by 

a.  =  (j/Zt/K)^  p.  ,  ?  =  [2k  JZ/L^[z+L)  ,  T  =  k  L  (2.2.2) 

110  o 
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In  the  above  is  related  to  the  Fourier  transform  variable  p  as  it  was 
in  Chapter  1  ,  i.  e. ,  ■  P  )'  •  By  using  these  independent  solutions 

in  Eq.  (A-13)  and  in  Eq.  (A— 7)  we  can  obtain  an  explicit  expression  for 
the  one -dimensional  Green's  function  and,  from  this,  an  integral  repreot  ro¬ 
tation  for  the  field.  The  representation  for  E  is  given  by 


E  =  E  +  E 

y  yf  yr 


(2.2.  3) 


where  E^  is  given  in  Eq.  (1.  2.  7)  and  represents  the  direct  field  from  the 

source  in  the  region  z  >0  .  The  second  term  in  Eq.  (2.  2.  3) ,  E  ,  is  given 

yr 

by 

i  V°  f  ik  ^p  (z+z')  +  px] 

IT  -  -  -  A  C  Prv  tO  O  A 


4n  J  p., 

-cd  r2 


(2.2.4) 


where 

r  =  -  A  /  Ab  (2.2.  5) 

with 

jE'(-aj  ,  0)  +  iaiE(-aj  ,  0)j  j  E*  (-a*  ,  0)  +  ia^f-a*  ,  0)| 

|  E'(-aj  ,  c)±ia2E(-a^  ,  c  )|  | E*(-a^  ,  c  )  ±  ia2 E*(-a^  ,  c  )| 

a2  =  (t/27A)^  P2  ,  p2  =  (1-p2)^  ,  c  =  (2t  .  (2.2.6) 


2  2 

In  the  above  expressions  for  A^  the  symbol  E;(-a^  ,  5)  =  dE(-a^  ,  ?)/d§  . 

The  square  roots  are  defined  as  in  Chapter  1  by  Fig.  1.  4.  The  integral 
2 

for  the  reflected  field,  Eq.  (2.  2.  4)  is  transformed  to  the  p^  plane.  The  result 
is  given  in  Eq.  (1.  2.  12)  where  the  integration  path  C  and  the  square  roots  are 
defined  as  in  Chapter  1,  Fig.  1.  5  and  1.6. 
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2.  3  Evaluation  of  Reflected  Field  for  Large  Layer  Thichness  (k^L»  1 


2.  3.  1  Geometrical- Optica 


Before  performing  the  asymptotic  evaluation  of  the  reflected  field  for 

k  L»1  we  shall  make  use  of  the  method  of  geometrical-optics  to  predict 
o 

the  ray  contributions  to  the  reflected  field.  Then,  in  the  next  section,  we 
shall  asymptotically  evaluate  the  integral  representation  for  the  reflected 
field  and  compare  the  two  results.  In  this  way  we  will  see  the  limitation 
on  the  theory  of  geometrical- optics  when  applied  to  a  parabolic  layer  . 

The  rays  emitted  from  the  source  can  be  divided  into  three  types  : 
direct,  transmitted  and  returning.  They  are  defined  in  the  same  manner  as 
in  Chapter  1  .  There  the  reflected  and  transmitted  ray  types  have  been 
separated  by  a  critical  ray.  This  is  also  the  case  here,  however  the  critical 
ray  has  a  different  behavior  as  is  shown  in  Fig.  2.  1  . 

The  trajectory  for  this  ray  is  obtained  by  integrating  the  ray  equation 


with  the  appropriate  value  of  p  for  the  critical  ray  (p=  /F)  .  The  result  is 


x  =  JTJi  V  +  Lin  (^l)3 


(2.  3.  2) 


where  we  see  that  as  z-*-L,  x-®. 

The  difference  between  the  behavior  of  the  critical  ray  in  a  linear 
transition  and  its  behavior  in  '  parabolic  transition  can  be  understood  better 
by  applying  the  reciprocity  prinicple.  Consider  a  ray  that  is  progressing 
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along  the  z  =  -  L  interface  of  a  linear  transition.  This  ray  will  be 
forced  to  turn  away  from  the  interface  since  it  has  a  finite  radius  of 
curvature,  (de(z)/dz|z  _  ^  ^  0  )  .  The  situation  is  not  the  same  for  a 

ray  travelling  along  a  parabolic  interface.  There  de(z)/dz|  =0 

and  the  radius  of  curvature  is  infinite.  In  this  case  the  ray  continues  to 
propagate  along  the  interface.  Now,  applying  the  principle  of  reciprocity 
to  the  parabolic  layer,  we  see  that  if  a  ray  starting  in  the  interface,  cannot 
escape  then  a  ray  outside  the  interface  cannot  become  tangent  to  the  inter¬ 
face  unless  it  has  an  infinite  radius  of  curvature.  The  critical  ray  in  the 
parabolic  layer  only  obtains  an  infinite  radius  of  curvature  as  x-“. 

As  in  Chapter  1  we  shall  now  focus  our  attention  on  the  returning  rays. 
By  substituting  the  parabolic  dielectric  variation  into  the  ray  equation  for 
returning  rays,  Eq.  (1.  3.  1),  and  integrating  it  we  obtain 

x  =  p(z'  +  z)/p  +  cosh"1  7^7  *  z>  0 

2  lpl  *  "  (2.3.3) 


x  =  pz'/p,  + 


cosh 


JL 

lpil 


pL 

7r 


cosh 


-1 


^(z+L) 

|PllL 


z  <  0 


where  =  i  |  p ^  |  and  e,  <p  <  1  .  A  typical  ray  is  shown  in  Fig.  2.  1  . 
Since  these  rays  cross  over  one  another,  we  shall  require  the  constraint 
equation  for  the  caustic.  This  is  obtained  by  taking  the  derivative  with 
respect  to  p  of  Eq.  (2.  3.  3).  The  result  is 


To  be  more  accurate,  we  mast  say  that  the  ray  is  not  exactly  on  the  inter¬ 
face  but  just  a  small  amount  above  it.  This  removes  the  ambiguity  in 
d€(,)/dz|i=  I 
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(*  +  */)/ cosh-1  -j^-|  -  2Lp2/|p1|2p  f  z>  0 

i7p2  +  jt  [co’h‘1  j^|  +  co,h'1  1 


(2.3.4) 


f  —  ♦ 

Lp,  r. .  .  . 


1 2  L  P 


[A(z+LJ  -  |p 


,  z<0  . 

J2L2]iJ 


The  caustic  is  obtained  by  choosing  a  value  of  p  and  solving  for  z  ;  then  by 

using  the  ray  equation  to  find  x  .  We  see  that  the  first  equation,  z>0  ,  can 

0 

be  solved  explicitly  for  z  while  the  equation  for  z<0  cannot. 

It  will  also  be  of  interest  to  know  the  number  and  the  location  of  the  foci 
formed  by  the  returning  rays.  The  constraint  equation  for  the  foci  is  obtained 
by  taking  the  derivative  of  Eq.  (2.  3.4)  with  respect  to  p  .  Upon  doing  this,  we 
obtain 


To  obtain  the  location  of  the  foci  we  must  solve  Eq.  (2.  3.  4)  in  conjunction  with 
Eq.  (2.  3.  5) .  This  will  give  us  the  z  coordinates  of  the  foci ;  then  the  ray 
equation  can  be  used  to  obtain  the  x  coordinate. 
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As  has  been  mentioned  in  the  introduction,  Orlov  b*  -  studied  these 
returning  rays.  To  be  more  specific,  he  has  made  a  det  t  analytic 
study  of  the  caustic  and  foci  for  z£0  .  He  has  then  plotted  the  caust^u  for 
some  typical  values  of  parameters  when  both  z  iJO  and  z<0  .  We  have 
redrawn  two  of  his  graphs  which  are  representative  of  the  caustic's  general 
behavior.  These  are  shown  in  Fig  .  2.  2.  In  the  upper  graph  (a)  we  see  that 
five  cusps  appear  on  the  caustic  while,  in  the  lower  graph,  two  of  the  cusps 
disappear  when  2A  >(1. 076)  L.  We  note  that  the  bow  tie  configuration  which 
was  observed  in  Chapter  1  occurs  twice  in  *  and  once  in  (b).  As  L 
becomes  smaller  the  remaining  bow  tie  configuration  disappears  and  only 
one  focus  is  left.  This  is  shown  in  Fig.  2.  3  where  the  caustic  has  been 
plotted  rather  than  just  sketched  for  L=  ,  1. 

The  caustic  appearing  in  the  three  previous  graphs  share  two  character¬ 
istics  in  common.  We  see  that  at  large  distances  from  the  source,  the  lower 
branch  of  the  caustic  approaches  the  interface  while  the  upper  branch  appears 
to  have  corrdinates  that  become  unbounded.  The  asymptote  to  the  lower 
branch  of  the  caustic  can  be  obtained  from  Eq.  (2.  3.  4)  with  z<0  if  we  assume 

that  p  -*0  as  z  -«0  .  By  usir  this  procedure  we  find  the  asymptote  to  be 
£ 

x  —  z'  yr/(2A(-z)  )*  .  (2.3.6) 

In  addition  we  can  also  find  an  asymptote  to  the  upper  branch  of  the  caustic 
by  assuming  that  as  -*0  ,  z-*®  .  By  using  this  assumption  in  Eq.  (2.  3.  4) 
with  z>0  ,  we  find  that 

x  «  //«j/^(z  +  z')  +  Lin  (  I  (2.  3.  7) 

z  +  z'  -  2L6jA /  IpJ2 


with 


(2.3.3) 


Fig.  2.3 

Caustic  Formed  by  a  Parabolic  Layer 
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We  see  that  aa  the  rays  are  emitted  closer  to  the  critical  angle  (p1  -0)  , 
tha  caustic  formed  by  these  rays  tends  clost.r  to  x  =  (*  +  *')  .  We 

note  that  this  equation  represents  the  reflected  critical  ray  from  an  abrupt 
transition. 

.When  L.  is  small  the  asymptote  given  in  Eqs.  (2.  3.  6)  and  (2.  2.  7) 
da  sc  noe  the  caustic  except  for  the  region  near  the  focus.  The  approximate 
location  of  the  focus  can  be  obtained  by  assuming 

*  *  -  L  +  A Ir^  ,  ipJ-BL172  .  ^  «1  (2.3.9) 

where  A  and  B  are  constants  tc  e  determined.  When  these  assumptions 
are  Used  in  Eqs.  (2.  3. 4)  and  (2.  3.  5)  they  lead  to  an  asymptotic  approximation 
for  the  focus.  The  approximate  location  .of  the  focus  which  is  obtained  is 
given  by 

x  =  JTJK  ts'  +  la  {%'  fit  3  (2.3.10) 


and 

x  =  -  L  +  ljtlfz  L^2  .  (2.3,11) 

As  t-«$  the  rays  that  make  up  the  lower  and  upper  branches  of  the  caustic 
appear  to  have  the  same  ray  trajectories  respectively  as  the  reflected  and 
lateral  rays  associated  with  an  abrupt  transition.  We  also  note  that  the 
foCUi  tends  toward  the  point  at  which  the  critical  ray  is  reflected  from  an 
abrupt  transition.  The  comparison  made  above,  says  nothing  about  the  ray 
contributions  but  only  states  that  the  two  layers  being  compered  have  similar 
ray  trajectories. 
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2.3.2  Asymptotic  Evaluation 


We  will  now  asymptotically  evaluate  the  formal  solution  for  the 
reflected  field,  Eq.  {2.  2.4)  when  k^L»l  and  compare  our  results  to  the 
prediction  of  geometrical-optics  given  in  the  last  section.  To  effect  the 
evaluation  of  the  integral  in  Eq.  (2.  2.  4)  in  physical  terms  we  will  expand 
the  reflection  coefficient  in  a  geometric  series  and  then  asymptotically 
evaluate  each  term  individually.  The  manner  of  development  of  this  section 
will  parallel  section  1.  3.  2.  quite  closely. 


The  reflection  coefficient  given  in  Eq.  (2.  2.  5)  can  be  written  as 


r  =  r  + 

oa 


r  rT  (r  -  r  > 

oa  L  ob  oc 

i  -  r  r 

L  ob 


(2.  3. 12) 


where 


*'  2  *  2 
E  (-a2 ,  o)  +ia}  E  (-at  ,  o) 

E'  (-a*.  o)  +  iajEf-aj,  o) 


(2.3.13) 


*  2  *  2 
E  (-alfc)  +  ia2E  ^  ,  c) 

roa ~  ”  .  2  *  2 

E  (-ar  c)  -  ia2E  (-ar  c) 


(2.3.  14) 


rob= 


E'(-aj  ,  c)  -  ia2 E(-a^  ,  c) 

*'  2  *  2 
E  (-aj  ,  c)  -  ia2 E  (-aj  ,  c) 


(2.  3.  15) 


r  = 
*  oc 


E'(-*j  ,  c)  +  ia2  E(-a ^  ,  c) 

* 7  ?  *  2 
E  (-alfc)  +  ia2E  (-aj.c) 


(2.3.16) 
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If  al°ng  the  integration  path  C,  the  denominator  of  the  second 

term  In  the  reflection  coefficient,  Eq.  (2.  3.  12) ,  can  be  expanded  in  a  geo¬ 
metric  ieries.  A  calculation  of  jr_  T  |  shows  that  in  the  intervals 

ob 

0*  *epj  <•,  hnpj  =  0  and  Re  p^  *  0  ,  0  <Im  <  JE ,  the  term  is  less  than  one, 
however,  for  the  interval  Rep.  =0,  •JE'ixn  p.  <®  ,  we  find  |r  T  J  *  1.  In 
Chapter  1  a  similar  problem  had  arisen  and  it  was  circumvented  by  deform¬ 
ing  the  integration  path  C  a  slight  amount  to  the  right  of  the  Im  p,  axis.  We 
will  perform  a  similar  deformation  of  the  path  C  at  this  point,  however  it 

mult  be  assumed  that  k  L»1  in  order  to  show  that  |T_  F  ,  |  <1  alone  the 

o  L  ob'  6 

deformed  path. 


We  will  make  use  of  the  series  expansion  in  Eq.  (2.  3.12).  The 


result  is 


<!9 


(2.3.17) 


Where  1^  is  defined  in  Eq.  ( i .  3.  14)  and  where  it  is  assumed  that  the 

r.  ,  T  ,T  .  and  T  of  chapter  1  are  replaced  by  those  of  Chapter  c. 

Li  oa  ob  oc 

By  substituting  the  series  form  of  the  reflection  coefficient  given  above 
in  the  integral  representation  for  the  reflected  field,  Eq.  (2.  2.  4)  and  inter¬ 
changing  Orders  of  summation  and  integration,  we  find  that 

OD 

E  =  E  +  Y  E  (2.3.18) 

yr  ro  L .  n 

7  n=0 

where  E  and  E  are  defined  by  Eqs.  (1.  3.  16)  and  (1. 3.  17)  respectively. 


Before  evaluating  E^q  and  Eq  by  the  method  of  steepest  descents  we 

will  asymptotically  approximate  1^  and  to  reduce  the  complexity  of  the 

Integrands.  These  expressions,  1^  and  1^  ,  are  composed  of  parabolic 

cylinder  functions  whose  as\  nptotic  approximations  along  the  integration 

*  *' 

path  C  appear  in  Table  2. 1.  The  functions  E  and  E  do  not  appear  in 
the  table,  however,  their  asymptotic  approximations  can  be  obtained  by 
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taking  the  complex  conjugate  of  the  asymptotic  expansions  for  E  and  E 


The  asymptotic  approximations  for  E  and  E*  will  be  needed  close 
to  the  integration  path  C  when  subsequent  path  deformations  are  made. 
These  asymptotic  approximations  can  be  obtained  by  analytically  contin¬ 
uing  the  formulae  in  Table  2. 1  if  we  exclude  the  regions  near  p^  =  0,  »/K  . 
In  these  regions  we  can  use 


-TTa*/4  +  in/8  +  i«  /2 
E(-a1  ,  x)  =  JTo  1 


U(-iaJ 


itr/4. 

«  ) 


(2.3.19) 


E*(-aJ 


-TTa^/4 


in/8  +  i*2/2 


U(+ia, 


+in/4. 
xe  ) 


(2.  3.  20) 


where  U  is  another  form  of  parabolic  cylinder  function,  x  =  0  or  x»c  and 
=  I*(l/2  -  ia^)  ,  $2  =  0.  Olver^*^  jB(  obtained  the  asymptotic 

expansions  of  U  and  U* (dU  (a,x)  /dx  =  U,(a,x)  )  for  large  a^  ,  and  any  x  and 
arga^  .  His  expansions  for  U  and  U'  can  be  used  near  Pj  *0  and  Pj  , 
however,  we  cannot  came  too  close  to  p^  =0  since  a^  *  (t/2^)^p  and  the 
asymptotic  parameter  a^  will  tend  to  sero.  The  functions  U  and  U'  could 
have  been  used  over  the  whole  integration  path,  however,  their  asymptotic 
expansions  are  moi?  complex  and  unwieldy. 

By  using  the  asymptotic  expansions  found  in  Table  2. 1  to  calculate 
ahd  ?a  in  the  region  Rep^wO,  ^<Impj  <•  ,  and  then  by  substituting 
these  approximations  in  the  appropriate  integrals,  we  find  the  E_^  and  Eq 
have  no  asymptotic  contribution  (koL»l)from  this  section  of  their  integra¬ 
tion  path  on  the  top  or  second  sheet.  In  a  similar  manner  by  the  use  of  Eqs. 
(2.  3.  19)  and  (2.  3.  20)  along  with  Olver's  asymptotic  expansions,  we  can 
show  that  the  integration  path  in  the  region  of  the  branch  point,  p^  -JK , 
gives  no  asymptotic  contribution. 
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symptotic  forms  of  parabolic  cylinder  functions  on  the  integration  path 
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If  we  again  use  the  asymptotic  approximations  from  Table  2. 1, 
we  find  that 

r  ~  i  A/4t  p^  (2.3.21) 

when  p^  is  on  the  integration  path  between  p^  =<y/cj  and  p^  =  i»/E  .  We  now 
use  this  approximation  in  the  integral  for  and  apply  the  method  of 
steepest  descents.  We  will  not  go  into  the  details  since  they  exactly 
parallel  the  asymptotic  evaluation  of  E  in  Chapter  1.  It  will  suffice  to 
•ay  that  th»»re  are  two  relevant  saddle  points  denoted  by  D  and  E,  as  before. 
The  contribution  from  E  is  exponentially  small  while  the  contribution  from 
D  can  be  interpreted  as  a  ray  reflected  from  the  z  =  0  interface  as  shown  in 
Fig.  1.  19.  The  asymptotic  order  of  these  contributions  are  shown  in  Table 
2.  2  . 

We  new  asymptotically  approximate  E^  .  Fir  jt,  the  formulae  in  Table 
2.  1  will  be  used  to  simplify  I*  along  the  portion  of  the  integration  path 
between  p^  =  </e^  and  Pj  =  i  •/&■  When  0<Re  Pj<ve~j,Impj  =  0  we  find  that 


r  ~ 

n 


^2n+l 


,2n+l  3n  4n+4  3n+2 
4  p2  P1  T 


i[2(n+l)8j  -nTt/2] 
e 


(2.  3.  22) 


and  when  Re  p^  =  0  ,  0  <  Ln  p^  <JE  we  find  that 


r 

n 


,n/2 


n  3n  n 

8  p2  T 


i2(nil)  3 
e  4 


(2.  3.  23) 


where  8j  and 


are  defined  in  Table  2.  1  . 
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1  . 


Type 

Asymptotic  Order 

A 

n 

o 

B 

n 

0(k-(3n+5/2) 

o 

D 

fM 

fO 

1 

o 

-CL  k 

E 

0(e  °),ol>0 

Table  2.  2 

Asymptotic  Order  of  Geometrical>Optic  Contributions 


The  above  approximations  are  good  in  the  vicinity  of  the  integration  path 

2 

except  near  the  origin,  i.  e.  f(r/2is/5)p1«  1  . 

K  we  now  replace  the  1^  appearing  in  the  integrand  of  Eq.  (1.  3..  17)  , 
by  its  asymptotic  approximation  and  find  the  saddle  point  equations  as  in 
the  .previous  chapter,  we  obtain 


>1  ("Hr  +  f  +  2(°+l)L.inh'1-^)  = 


(2.  3.  24) 


and 


KRepj  <  ,  Impj  = 


^  -1''1  ipTf)  = 0  <—> 


Jt  |Pil 

Repj  =  0  .  0<lmp1<v^ 
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where  the  upper  sign  in  the  above  equations  refers  to  the  top  sheet 
while  the  lower  sign  refers  to  the  second  sheet. 

Upon  obtaining  the  saddle  point  locations  from  the  above  equations, 
we  deform  the  contour  C  through  the  saddle  points  and  into  the  decay 
regions.  Since  the  contribution  due  to  the  portion  of  the  integral  lying 
in  the  decay  region  is  exponentially  small,  the  integral  will  be  asymp¬ 
totically  equal  to  the  sum  of  the  saddle  point  contributions.  The  contribu¬ 
tions  arising  from  saddle  points  obeying  Eqs.  (2.  3.  25)  and  (2.  3.  24)  will  be 
known  as  type  and  contributions,  respectively.  It  can  be  shown  that 
no  relevant  saddle  points  occur  on  the  second  sheet  for  any  n  and  therefore 
the  deformed  path  will  be  the  same  as  that  shown  in  Figa.1.13  and  1. 14. 

We  shall  now  focus  our  attention  on  the  relevant  saddle  points  on  the 
top  sheet  and  their  contributions  to  the  reflected  field.  The  integral  Eq  will 
be  considered  first,  since  it  contributes  the  dominant  terms  to  the  reflected 
field  for  large  k^L  . 

We  will  now  investigate  the  saddle  points  of  type  .  A  comparison 
of  Eqs.  (2.  3.  3)  and  (2.  3.  25)  shows  us  that  the  two  are  identical.  This  means, 
as  in  Chapter  1,  that  each  saddle  point  corresponds  to  a  returning  ray.  If 
we  now  recall  our  discussion  of  returning  rays,  we  will  remember  that  there 
are  two  basic  regions  divided  by  a  caustic.  To  the  right  of  the  caustic  there 
are  two  returning  rays  ;  to  the  left  there  are  none.  For  certain  parameters, 
the  caustic  can  assume  a  bow  tie  configuration  for  z>0  as  shown  in  Fig.  2.  2(a). 
Inside  the  bow  tie  there  are  four  returning  rays  while  outside  the  bow  tie  there 
is  the  same  number  of  returning  rays  as  before.  With  the  above  information 
we  see  that  there  are  0,  2  or  4  solutions,  depending  upon  the  location  of  the 
observation  point,  to  Eq.  (2.  3.  25)  . 
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We  shall  next  examine  saddle  points  of  the  type  Bq  .  An  investigation 
of  Eq.  (2.  3.  24)  shows  that  there  is  one  saddle  point  for  any  given  observa¬ 
tion  point.  The  saddle  point  contribution  can  be  interpreted  as  a  ray  reflected 

-5/2 

from  the  lower  Interface.  The  amplitude  factor  of  this  ray  has  0(kQ  )  . 

As  stated  previously,  the  reflected  field,  E  ,  is  composed  of  a  sum  of 

o 

Aq  and  Bq  type  contributions.  The  deformed  path  on  the  top  sheet  is  shown 

in  Figs.  1. 1  5  and  1. 16  when  there  are  no  saddle  points  or  there  are  two 

saddle  points  of  the  A  type,  respectively.  When  there  are  four  saddle 

o 

points  of  type  Aq  the  path  is  similar  to  the  ones  shown,  but  it  is  not  shown 
here.  An  approximation  of  the  deformed  path  shows  that  it  can  be  deformed 
around  the  origin  for  any  given  observation  point  and,  therefore  no  contri¬ 
bution,  similar  to  type  Cq  ,  arises  as  we  observed  in  Chapter  1  .  This 
means  that  there  is  no  lateral  wave  contribution.  However,  there  are  two 
(assume  that  the  bow  tie  is  not  present)  rays  through  each  point  to  the  right 
of  the  caustic  while  in  a  similar  region,  in  Chapter  1,  there  is  only  one  ray 
through  each  point.  When  the  observation  point  is  located  far  from  the  source 
compared  to  layer  thickness  we  observe?  first,  that  the  returning  ray  in 
Chapter  1  has  a  similar  trajectory  to  one  of  the  returning  rays  of  Chapter  2; 
and  second,  that  the  trajectory  of  the  lateral  ray  in  Chapter  1  is  similar  to  the 
trajectory  of  the  other  returning  ray  observed  in  Chapter  2  . 

We  will  now  evaluate  explicitly  the  contribution  due  to  the  ray  whose 

trajectory  is  similar  to  a  lateral  ray.  The  evaluation  will  be  limited  to 

observation  points  which  are  located  far  from  the  source  compared  to 

layer  thickness  and  also  those  observation  points  where  L  »L  .  Here  L 

P  P 

is  defined  as  in  Chapter  1.  The  last  condition  has  the  following  physical 

explanation.  When  the  observation  point  is  located  far  from  the  source,  the 

upper  branch  of  the  caustic  is  approximately  described  by  L^  =  0  .  The 

condition,  L  »L,  says  that  the  observation  point  should  be  many  layer 
P 

thicknesses  from  this  portion  of  the  caustic. 


-61- 


By  using  the  above  approximations  on  the  saddle  point  equation, 
Eq.  (2.  3.  25)  ,  we  obtain  the  approximate  location  for  the  two  saddle 
points  in  question.  They  are  : 


Ip,J  *  zJE 


-JEi.  /2VTl 

p  * 


(2.  3. 26) 


The  notation  used  is  motivated  by  the  fact  that  the  contributions  from  saddle 
points  p^  and  p^  are  similar  to  reflected  and  lateral  rays  from  an  abrupt 
transition.  The  first  saddle  point  p ^  is  obtained  by  assuming  that  the  last 
term  in  Eq.  (2.  3.  25)  is  small  while  the  saddle  point  p  is  obtained  by 
assuming  p^  «0  .  From  Eq.  (2.  3.  25)  we  see  that  there  appears  to  be  a  saddle 
point  at  Pj  =  0  ,  however  the  asymptotic  approximations  used  to  simplify  the 
integrand  are  not  valid  there. 

Now  by  evaluating  the  integral  E>  at  the  saddle  point  p^  we  obtain 


1/4 

L  e  **  ik  +  r/)+  Jc.  x  +  JK 'Ll  +  i  3n/4 

2>c,TX  e  °  1 

1  o 


where  E^  is  the  contribution  from  Eq  which  is  due  to  the  saddle  point  . 

The  phase  of  the  above  contribution  corresponds  to  a  lateral  ray,  however 

the  amplitude  is  of  0(k  ^)  which  makes  this  contribution  a  classical  geo- 

o 

metrical-optic  result. 

An  examination  of  Eq.  (2.  3.  26)  reveals  one  basic  difficulty.  As  L  / L 

becomes  increasingly  large  |p1T  |  -  0  .  In  order  for  the  asymptotic  forms 

2 

of  Table  2.  1  to  be  valid  a^»l  or  using  Eq.  (2.  3.  26)  we  find 

-  jel  njr.  l 

p  1 


T  e 


»1  . 


(2.3.28) 
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TUa  inequality  essentially  says  that  as  1W  L  gets  very  large  compered  to 
T  ,  the  geometrical  -  optic  result ,  appearing  in  Eq.  (2.  3.  27)  ,  is  no  longer 
valid.  We  see  rays  emitted  from  the  source  close  to  the  critical  angle 
resemble  geometrical- optic  rays  with  lateral  ray  type  trajectories,  but, 
as  the  emission  angle  becomes  closer  to  the  critical  ray,  the  geometrical- 
optic  interpretation  can  no  longer  be  given. 


An  evaluation  of  E  for  arbitrary  n  will  show  E  to  be  the  dominant 

term  in  the  series.  We  will  not  go  into  detail  since  the  evaluation  exactly 

parallels  a  similar  development  in  Chapter  1  .  It  will  suffice  to  say  that 

the  contributions  of  types  A  and  B  will  be  encountered.  Those  of  type 

n  n 

A  are  returning  rays  reflected  from  the  z  =  0  interface  n  times.  On  the 

n 

other  hand,  the  Bq  contributions  are  those  rays  reflected  n  + 1  times  from 
the  a  =  -  L  interface.  The  asymptotic  order  of  those  contributions  is  shown 
in  Table  2. 2  . 


2. 4  Evaluation  of  Reflected  Field  for  Large  k  L 
_ 2_E 


In  this  section  we  will  investigate  the  reflected  field  from  the 

parabolic  transition  when  k^L,  »1  and  L^>>L.  From  this  investigation 

we  hope  to  learn  what  the  pertinent  field  contributions  are,  when  the  layer 

is  thick  compared  to  wavelength,  i.  e.:  k^L  » k^  L  » 1  .  In  addition,  we 

want  to  relate  these  contributions  to  the  field  constituents  that  occur  when 

k  L>>k  L  >>  1  and  when  k  L  » 1  ,  k  L«  1  . 
o  o  p  op  o 

The  integral  representation  for  the  reflected  field  is  given  by  Eq.  (1. 4. 1) 

where  the  reflection  coefficient  used  there  is  replaced  by 


r 


(2.4.  1) 
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whirt 


|w'(-ej,0}-  ia.W(-aj,0)| 

j-W'(-a*.  0)  -  ia^f-a^.O)! 

V  * 

bn 

*  # 

|w'(-aj,-c)+  ia2W(-aj,  -c)j  j 

1  r 

j-W'(-a*.  c)  +  ia2W(-aj  .  c)j 

W(>a^  ,  at?)  are  another  form  of  parabolic  cylinder  functions  which  are 

*  la  at  a  a 

defined  by  Abramowits^^  end  w'{-a^  ,?1=  dW  {-a ^  .  We  fee  that  the 

integral  representation  employed  here  is  the  same  as  that  used  in  the  last 
section,  except  with  a  different  reflection  coefficient.  This  reflection  coeff¬ 
icient  is  derived  via  Appendix  A.  by  choosing  ♦  (*)  =  W(-a7»  ~>  and 

2  i  a 

=  W(-aj  ,  ?)  where  ?  has  been  defined  in  the  beginning  of  Chapter  2. 

This  alternative  form  of  reflection  coefficient  makes  the  computations,  that 
will  follow,  easier  since  the  W{-*~  ,  are  real  functions  . 

Jf  w-  now  use  the  assumption  that  k  L«  k  L  then  the  reflection  coefficient 

o  o  p 

in  Eq.  (2.4.  1)  is  slowly  varying  compared  to  tl  .*  exponential  in  the  integrand, 
and  the  standard  steepest  descent  techniques  can  be  used  to  asymptotically 
evaulate  the  integral  representation.  We  will  not  go  into  detail  since  the 
development  parallels  Chapter  1,  Section  1.4. 

The  results  of  this  asymptotic  analysis  are  that  there  are  two  saddle 

point  contributions  which  make  up  the  reflected  field.  The  firs;  of  these  is 
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s imply  the  reflected  ray  contribution.  The  contribution  has  an  0(kQ  )  for 
all  layer  thicknesses  considered  here. 

The  second  contribution  to  the  reflected  field  comes  from  the  saddle 
point  occuriAg  at  p  (  =0  and  it  reduces  to  the  lateral  wave  contribution  on  an 
abrupt  interface  when  k^L-O  .  The  contribution  will  be  denoted  by  and 
it  is  given  by 
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el“ 


ua  ik :  C^SC* ♦  *')  +  v'cTjc]  +  ittVS’r)  +  i3n/4 

c|MA(^T)e  °  1 

JE?t(koL)3,Z 


(2.  4.  2) 


where 


Ja|  dT 
Tjdpi 


pr° 


rwgni/4)  [W#(Q,c)W(0.-c)  4  WW.clW'fO.  -c)] 


r(V4)  2  c2  2 
[(w'(o,-c)- w'(o,c))  +  J  (W(0.c)  +  W(0.-c)r3 


(2.  4.  3) 


♦(VS-T,  =  2  tan”1  U  -Xig.sLl 

flV®  1  Lc  W(0,  c)  +  W(0,  -c)  J 


(2.  4. 4) 


The  amplitude  A(^t)  and  tike  phase  f  (^/SV)  have  been  plotted  by  computer 
and  are  shown  in  Figs.  2.  4  and  2.5.  In  these  figures  the  thick  layer 
approximations  to  A(^V)and  t(^T)  are  also  shown  . 

The  physical  interpretation  of  E.  can  best  be  found  from  an  investigation 

of  the  small  and  large  t  limits .  When  t  is  small,  as  we  have  already  stated, 

E.  reduces  to  the  lateral  wave  on  an  abrupt  interface.  We  see  this  from  the 
Id 

fact  that  as  T-'O,  A(v^r)-1  and  t(<^T)-*0.  When  k  L»l,  or  the  layer  is 

o 

large  compared  wavelength,  reduces  to 


r2  (1/4)  « J/4  L1/2  i  ko[^S"(*  +  *')  +  JTX  t  +  L^T ]  +  iTT/2 

zl  e 

2(2n)  A  Lp  (kQLp) 


.  (2.4.5) 


H  we  compare  the  lateral  wave  contributions  excited  on  a  linear  and 

■* 

parabolic  layer  for  k^L  large,  as  given  in  Eqs.  (1.  4. 10)  and  (2.  4.  2) 


A(7Sr) 


Fig.  2.4 

Lateral  Ware  Amplitude  vs.  Normalised  Transition  Thickness  to  the  One  Half  Power 


♦Vat* 


Lateral  Wave  Phaae  v*.  Normalized  Transition  Thickness 


*espectively,  we  observe  that  both  have  a  similar  phase  dependence  and 
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an  amplitude  dependence  on  distance  of  0(L^  )  .  The  outstanding  differ¬ 

ence  between  the  two  contributions  is  that  the  lateral  wave,  excited  on  a 

linear  layer,  has  an  excitation  coefficient  of  Ofk  "^)  while  the  lateral 

°  -1 

wave,  found  on  a  parabolic  layer,  has  a  coefficient  of  0(k  )  .  We  see 

o 

for  high  frequency  that  the  parabolic  layer  wave  is  excited  more  strongly. 


The  lateral  wave  contribution  given  in  Eq.  (2. 4.  5)  also  sheds  some  li^ht 
on  another  area.  Id  the  previous  section,  a  geometrical- optic  ray  contribu¬ 
tion  was  found,  ft  has  a  ray  trajectory  similar  to  a  lateral  ray,  however, 
the  observation  points  move  away  from  the  source,  L  »L  ,  the  geo- 
metrfcal- optic  contribution  becomes  invalid.  We  see  now  that  this  contribu¬ 
tion  transforms  into  a  true  lateral  wave  contribution  when  k  L  »k  L»l. 

op  o 

1*110  transition  between  the  two  regions  of  validity  for  the  E  contribution 

L 

is  quite  complex  and,  as  yet,  transition  functions  have  not  been  found. 
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CHAPTER  3. 

THE  SYMMETRICAL  EPSTEIN  TRANSITION  LAYER 
3. 1  Introduction 

b  Chapter*  1  and  2  the  reflected  field  from  *  linear  and  parabolic 
transition  ha*  been  studied.  There  it  was  found  that  a  lateral  wave  was 
excited  for  all  layer  thicknesses,  and  that  its  character  was  chauged  *s 
tile  normalised  layer  thickness,  k^L  ,  was  varied.  For  both  of  the  above 
mentioned  cases  the  thick  layer  limit  was  particularly  interesting.  In  this 
limit  tike  lateral  wave  discussed  in  Chapter  1  appeared  to  be  excited  at  the 
point  that  the  critical  ray  was  tangent  to  the  lower  interface,  hi  Chapter  2, 
no  such  interpretation  was  possible,  but  instead,  the  lateral  rays  resulted 
from  returning  rays  emitted  close  to  the  critical  angle. 

In  an  examination  of  the  above  problems,  we  bad  noted  that  the  continuity 
of  the  layer  at  the  lower  interface  was  the  critical  factor  in  determining  the 
type  of  lateral  wave  that  could  be  excited.  A  logical  extension  of  the  studies 
undertaken  in  tike  first  two  chapters  could  consist  of  an  investigation  of  the 
reflected  field  from  a  layer  that  was  completely  continuous  at  the  lower  inter- 
face Wave  functions  for  the  layer  mentioned  above  could  not  be  found,  there¬ 
fore  we  rot  *rted  to  the  study  cf  a  layer  of  infinite  extent.  The  symmetrical 
Epstein  layer,  to  be  considered  in  this  chapter,  is  such  a  layer.  Its  dielectric 
profile  is  given  by 


e  (a) 


1 

2 

Cj  +  A  sech  z!  L 


s  >0 


<0 


A  “  1  -  e  .  (3.1.1) 


Tht.  layer  has  the  property  that  as  z  --  •  ,  e(s)  >  t,  and  all  derivatives  of 
c(s)  approach  aero.  Thus,  we  see  that  the  layer  is  completely  continuous 
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at  tfaa  lower  homogeneous  medium  which  occurs  at  minus  infinity.  The 
dielectric  layer  also  has  the  required  property  that  as  L-»6  the  layer 
approaches  an  abrupt  interface. 

In  this  chapter  we  will  first  investigate  the  structure  of  the  returning 
ray  trajectories  when  k^L»l.  Following  tide,  a  uniform  asymptotic 
approximation  for  the  reflected  field  will  be  derived  which  will  be  valid  for 
abritrary  values  of  k^L .  By  means  of  this  approximation  we  can  relate  the 
lateral  wave  on  a  sharp  interface  to  the  geometric- optic  rays  observed  when 
the  layer  is  thick  compared  with  wavelength. 

3. 2  Formal  Solution 

The  integral  representation  for  the  reflected  field  in  the  region  z  >0  is 
obtained  in  a  manner  similar  to  that  used  in  Chapters  1  and  2 .  The  source 
configuration  is  the  same  as  before  and  again  only  the  component  of  the 
electric  field  is  excited.  This  component  satisfies  the  wave  equation, 

Eq.  (1*  2. 1),  where  the  source  amplitude  has  been  assumed  to  be  normalized 
as  In  Eq.  (1.  2.  2) .  The  representation  zor  is  obtained  by  the  application 
at  Fourier  integral  techniques  to  Eq.  (1. 2. 1).  The  result  is  given  by  a  sum 
of  direct  and  reflected  fields 


E  *  E  ,  *  E 
y  yf  yr 


where 

E 


1  i*  Jkotp2U-z'|  +p*]  dp 
yf  “  ”  4TTi 


.00 


z  >0 


(3.2.2) 


and 
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yr  4tt 


+  ®  ^  ik  [p  (z  +  z  )  +  px]  _ 

r  o  2  dp  ,  z>0 

e  ~ 


L-  — 

T  i  .1  p_ 
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(3.2.  3) 


I 


with 


*k0P2  +  e'(o)/8(o) 


*0*2  ‘  0  C°l/©<0) 

p2-(:-p2)*  .  p,  -  («,  -  p2)4  . 


(3.2.4) 


(3.  2.5) 


Th«  function  0(z)  used  In  the  reflection  coefficient  T  is  defined  as  the  solution 

of  the  differential  equation 
2 

[—  +  kc  (A  *ech  (*/L)  +  Pj  )  1  8(«)  =  0  (3.2.6) 

da  ** 

which  satisfies  the  radiation  condition  as  z  —  -  “ .  The  derivative  of  9(z)  with 
respect  to  x  will  be  denoted  by  8'(k).  In  order  to  dofine  the  integrals 
appearing  in  (3.2.2)  and  (3.2.3)  completely,  the  square  roots  p  and  p 

*  Cm 

most  be  specified.  This  specification  is  shown  in  Fig.  1.  4  where  the  four 
branch  cuts  have  been  chosen.  The  integration  path  P  is  also  shown  in  the 
same  figure. 


We  will  bow  solve  (3.  2.  6)  exactly  Trarsforming  (3.  2.  6)  by 


?  *  tanh  (z/L) 


(3.2.7) 


reduces  it  to 


{u-S2)  —  '  +  [v<v+D  -  m2/(1-§2)]  }  ®=o 


(3.2.8) 


2  t 

n  =  +iT?i  ,  v=-  1/2  +  (1/4  +  t  ar  ,  T  =  kQL  .  (3.2.9) 


This  is  the  Associated  Legendre  equation  whose  two  independent  solutions  are 
(?)  ♦  Q  ^  (§)♦  If  we  represent  (§)  as  an  hypergeometric  function^23  and 
let  s^***  we  see  that  it  obeys  the  radiation  condition.  The  solution  to  (3. 2.  6) 


is  then 


e(x)  *  £  -  tanh(z/L)j 


(3.2. 10) 


The  initial  value  (z=0)  of  this  Associated  Legendre  function  and  its  derivative 

(24) 

are  given  in  terms  of  Gamma  function  .  Using  these  relations,  one  obtains 


0'(O)  2T(v/2  +  y/2  +  1)  T  (v/2  -  U / 2+ 1 )  ,  P  ,  ii4/;i 

8  (0)  "  '  LP (v/2+  u/2  +  l/2)r(v/2-u/2+l/2)tanLTT(v  +  W)/y 


8/(0)  _  2r(l/2  -  v/2  -  U/2)  r  (v/2  -  U/2  +  1) 

8  (0)  Ln-v/2-w/2)  r  (v/2  -  u/2  4-  1/2) 


(3.  2. 11) 


(3.2.12) 


In  order  to  determine  the  asymptotic  properties  of  the  reflection  co¬ 


efficient  the  location  of  the  poles  of  the  gamma  function  will  be  necessary. 


Consider  the  poles  of  T  (v/2  +  4/2  +  1),  for  example.  They  occur  when 

1  2  4 

-1/4+-J  (1/4+  1  LV  +  iTP]l/2+  1  =  -n  ,  n=«,  1,2,  ...  (3.2.13) 

is  satisfied.  Since 

T  (A)^  <  3/2  +  (1 /4  +  t2A)^  (3.2.14) 

for  all  r(A)^  ,  the  first  pole  lies  in  the  range  Rep^  *0  ,  (A)^  <Impj  .  Poles 
corresponding  to  n=  1  and  higher  have  larger  values  of  |p^  |  and  therefore 
occur  in  the  same  range  as  the  first  pole.  The  pole  locations  of  the  various 
gamma  functioca  appearing  in  (3.  2.  11)  and  (3.  2. 12)  are  shown  in  Table  3. 1  . 
Since  these  poles  occur  on  the  Rep  axis  for  |Rep|  >(e^)^  they  lie  directly  on 
the  integration  path.  If  a  small  amount  of  loss  is  added  to  the  dielectric 
medium,  the  poles  with  Re  p  >  0  shift  into  the  first  quadrant  of  the  p  plane 
while  the  remaining  poles  shift  into  the  third  quadrant.  W ith  this  information  , 
the  integration  path  can  be  deformed  around  the  poles  in  the  correct  manner  in 


the  limit. 6f  zero  loss. 


Gamma  Functions 

Pole  Locations 

r<v/2  +  u/2  +  l) 

Re  Pl  =  0 

F 

Im  Pj><A) 

r(v/2  +  m/2  +  1/2) 

Re  Pj  s  0 

Im  p  >(A)^ 

r(v/2  -  u/2  +  1) 

Re  Pl  =  0 

i 

Im  Pj<-(A)b 

r.(v/2  -  U/2  +  1/2) 

Re  Pl  =  0 

* 

ImPl<-(A) 

r(I/2-  v/2  -  U/2) 

Re  Pj  =  0 

4 

Im  Pj<(A) 

H-v/2  -  u/2) 

Re  Pl  =  0 

Im  Pj<(A)^ 

Table  3. 1 


role  Locations 


¥ 


L 
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Before  proceeding  we  shall  check  the  convergence  of  the  integral. 

To  do  this  we  asymptotically  evaluate  (3. 2. 10)  as  |p|  -•  .  We  find 

9  /(0)/9(0)  -  ko|pj|  (3.2.15) 

which  gives  a  reflection  coefficient 


_  !p2!  -  Ip,I 

r  = 


(3.2.16) 


|p2l  ♦  IpJ 

when  p  is  large.  Since  the  exponential  in  the  integral  decays  for  large  p,  the 
integral  converges. 

To  simplify  future  calculations,  the  branch  points  at  p  =  ±(e ^ are 


now  removed  by  means  of  the  transformation 

p  =  ,vpi2,i  •  p2  =  <i  +  p,2>4  • 


(3.2.17) 


This  transforms  the  integral  in  (3.  2.  3)  to 
E 


1  p  pir  ik<JpZ*z  +  I*  +  p*^ 

J  TT,  e  dpl  > 

C  2 


where  the  multivalued  functions  p|  p^  are  defined  on  a  four  sheeted  Riemann 
surface.  The  first  two  sheets  of  this  surface  are  shown  in  Figs.  1.  5  and  1.  6  . 
The  multivalued  character  of  the  mapping  (3.  2.  17)  leads  to  a  transformed  inte¬ 
gration  path,  C  ,  in  the  p^  plane.  Parts  of  this  path  appear  on  different  sheets 

since  any  two  points  symmetrically  located  on  path  p  lead  to  the  same  value  of 

2 

pi  • 

3f  3  Ray  Trajectories 

When  the  parameter  L  of  the  Epstein  layer  is  large  in  comparison  with 

a  wavelength,  the  medium  can  be  considered  slowly  varying,  and  the  methods 

of  geometrical  optics  can  be  used  to  find  the  ray  trajectories.  In  this  section, 

the  ray  family  will  be  found  by  using  the  ray  equation  for  a  stratified  medium 
(25) 

given  in  Kelso  .The  caustics  of  this  family  of  rays  will  also  be  found  and  the 
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bAitio?.  of  Om«  CMrtki  uL1*  0  will  bo  investigated  mad  rolotod 
to  (ha  half  ipoco  problam. 

>  • 

First,  tha  reflected  rayt  So  tha  homogaaaoaa  m odium,  m  >  0,  wiU  bo 
considered.  Tbooo  rayt  oro  straight  limes  (iron  by 

n/p  -(*  +  *')/  p2  -  Ln/Jpj  |  «  0  (3.3.1) 

where  p^  io  imoginory  for  rofloctod  rays,  i.  o. ,  p1  =  i|p  |.  The  same  ray 
equation  io  recovered  from  the  exact  solution, (3.2. ]  8)  ,  if  T  is  replaced  by  its 
asymptotic  approximation  (kQ  L  » I )  and  then  the  saddle  point  condition  is 
applied  to  the  resulting  integrand.  The  caustic  is  found  by  eliminating  fp^l 
between  (3.  3. 1)  and  the  derivative  of  (3.3.1)  with  respect  to  |pj  .  We  obtain 

x2/3  .  (y*)1'3  ♦  ,')2/3+  |W/4J'3.  (3.3.2, 

This  caustic  coincides  with  the  critically  reflected  ray  found  in  the  half  space 
problem  in  the  limit  as  L  «*  0. 

The  rays  in  the  medium  z  <  0  are  divided  into  three  types:  transmitted 
rays;  critical  ray;  and,  totally  reflected  rays.  The  transmitted  rays  are 
not  of  great  interest  to  us.  The  critical. ray  is  given  by 

x  =  £s/  -  L  sinh(z/L)J  ,  IpJ  =  0.  (3.3.3) 

This  ray  never  turns  since  z  is  a  monotonicaliy  decreasing  function  of  x.  The 
totally  reflected  rays  satisfy  the  equation 

coew  s  1  -  2|pj  i  2P2”2  *inh2(z/L)  (3.  3. 4) 

w  =  2  |pt  j(pL/_1(x-xo),  Xq=  p  z'/p2>  0<w<  2tt  . 


where 
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These  rays  sntsr  the  inhomogeneous  medium  at  x  =  x  ,  turn  at  x  =  x  +  x  . 

o  o  t 

(x  *  TrpL(2|p  I)'1)  and  return  to  the  linear  medium  at  x  =  x  +  2x  .  A 
v  1  /  o  t 

typical  ray  is  shown  in  Fig.  3.1 .  The  rays  have  been  plotted  for  some  typical 
parameters  in  Fig.  3.  2  . 

The  caustic  formed  by  these  returning  rays  is  found  by  eliminating 
|pj  |  between  the  derivative  of  {3.  3.4)  with  respect  to  |  |  ,  i.  e. , 


P3!- 


V 


Ip, Ip, 


tan  ( w/2 ) 


(3.3.5) 


and  the  ray  equation,  (3i3«4)  .  The  caustic  is  shown  in  Fig.  3.  2  .  It  consists  of 
two  branches:  one  that  crosces  the  z  =  0  line  and  connects  to  the  caustic 
in  the  homogeneous  space,  and  another  that  tends  asymptotically  toward  the 
z  =  0  line  as  x  *•  ®.  The  ray  and  constraint  equations  can  be  solved  exactly 
for  the  exit  point  of  the  first  branch  into  the  homogeneous  medium.  A 
calculation  shows  the  exit  point  lies  on  the  beginning  of  the  free  space  caustic. 
The  asymptotic  character  of  the  second  branch  can  be  found  by  assuming 
p  «1.  Using  this  approximation  in  the  ray  and  constraint  equations  gives 
the  asymptote 


-1 


(3.  3.  6' 


z/L  =  -  z 

This  asymptote  tends  toward  the  interface  as  L  -*  0. 

As  the  two  branches  progress  deeper  into  the  inhomogeneous  medium 
they  meet  in  a  cusp.  The  location  of  this  cusp  can  be  obtained  with  the 
additional  constraint  equation  found  by  taking  the  derivative  of  (3.  3.  5)  with 
respect  to  |pj  I  .  The  result  is 


Fig.  3.2 

Reflected  Rays  and  Caustics  From  an  Epstein  Laye 
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wher* 


..  ,  A  -lUr(v/2  +  U/2  +  1)  r(W2-u/2  +  1)  rv 

♦(Pj.(A)  )=  Un  jTp2r(v/2 -H/2+l/2r(v/2+u/2+l/2)tanL2(vfU)j(,{3-4*3) 

The  multivaluednes s  of  the  function  tan  ^  does  not  appear  in  the  overall  integrand 
2itan*1 

since  e  is  single  valued. 


The  uniform  asymptotic  evaluation  of  the  above  integral  is  performed 
by  first  finding  the  saddle  points  for  all  values  of  0<t  <  »  .  In  general,  these 

saddle  point  locations  will  depend  on  the  large  parameter  k  and  the  normalized 

o 

length  t ,  The  location  of  these  saddle  points  can  be  found:  first,  by  the 
method  of  successive  approximations  and  second,  by  using  an  asymptotic 
approximation  to  the  function  f  for  large  T.  The  first  method  is  a  generalized 
sharp  interface  technique  and  is  valid  when  the  length  L.  is  small  compared  to 

z'  and  k  is  large,  while  the  second  method  gives  the  geometric-optic  result 

d 

and  holds  when  T  is  large.  Since  the  two  methods  overlap  for  a  restricted 

range  of  t,  the  saddle  point  locations  can  be  found  for  any  t. 

The  saddle  point  equation  is  given  by 

p  [( z  +  z')  /p  -  x/p  ]  +  k  Sdf/dp, )  =  0  (3.  4.  4) 

l  c  o  i 

on  the  top  sheet  and 

Pj  [(z  +  z')  /p2  +  x/p  ]  +  kQ  \d<f/dpj)  =0  (3.  4.  5) 


on  the  second  sheet  First,  we  will  apply  the  method  of  succensive  approxi¬ 


mations  to  (3.  4.  4)  and  (3.4.5).  The  solutions  to  (3.4.  5)  neglecting  the 


1 


perturbing  term  kQ  (dijr/dpj)  are 


o 

(0) 

1L 


(0)  r 

P,,  =0  .  P,r=±il 


x  A  -  (z+z 


(3.4.6) 


l 
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where  the  saddle  points  p  and  p  {+)  correspond  to  the  lateral  and  reflected 

lL  1  x 

wares  found  in  the  half  space  problem.  The  superscript  is  used  to  denote  the 
order  of  approximation.  The  first  approximation  is  obtained  by  evaluating 
the  perturbing  term  in  (3. 4 . 4)  with  the  unperturbed  values  given  in  (3.  4.  6)  and 
then  assuming  p^  =  ^  +  6p^  with  6p^  «1.  The  first  order  saddle  points 


P(1L  =  f  Sdf/dpj)  ,  L=x~{Cl /*)*(*+ z')  (3.4.7) 

Pl=  ° 

Pli)apS,’{[(A,\V*/P23  +  X/l>3)(X(A)*+‘(Cl)*)]  <d*/dPi,}[p  =p<°>  {3‘4  8) 

rl  ~  1  r 

where  L  is  the  distance  the  lateral  wave  travels  along  the  interface  in  the 
P 

half  space  problem.  The  approximate  method  is  only  valid  when  «1 . 

For  this  to  be  true  it  is  sufficient  to  have 


L  ) 

o  p 


j  I  d*(Pl,T(A^) 


^ -  «1»  0<T  <  •,  Rep^=  0,  0<Impj<(A)  .  (3.4.9) 


w 

One  can  show  for  small  and  moderate  values  of  T  that  dt(p,(A)  t)  /  is 

1  /  apj 

bounded  except  at  p^=  0  and  for  large  values  of  t  that  df  ^dp^M-inT, 
With  the  aid  of  the  aboye  information,  ws  see  the  method  of  successive 


approximations  is  valid  when 

k  L  >>1 

o  p 


L/L  <<1,  0<T<<*>. 

/  p  - 


(3.4.  10) 


When  the  method  of  successive  approximation  is  applied  to  (3.  4.  5) ,  one  saddle 
point  is  found  on  the  negative  imaginary  axis  near  the  origin. 

The  second  method  of  finding  the  saddle  points  is  to  use  the  asymptotic 
approximations  for  the  Gamma  functions  in  (3.  4.4)  and  (3.  4.  5).  The  resulting 
saddle  point  equations  are 


- DEFORMED 

PATH 


Fig.  3.3 


Top  Sheet  of  Riemann  Surface-Decay  Regions 


Fig.  3.4 


Second  Sheet  of  Riemaan  Surface-Decay  Region* 
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/  z  +  z  —  x 

piV~T~  + 


f)*lr,U--0.  {; 


-  top  sheet 
second  sheet 


(3.4. 11) 


Here  again,  we  obtain  a  similar  arrangement  of  saddle  points  to  those  obtained 
by  the  method  of  successive  approximation.  Three  saddle  points  appear  on  the 

imaginary  axis  of  the  top  sheet.  Two  occur  for  Q<  Imp^<  (A)^  and  one  occurs 

i 

for  -(A)  <  Imp^  <  0.  On  the  bottom  sheet  we  find  only  one  saddle  point  on 

* 

the  imaginary  axis  for  -(A)  <  Imp^  <  0.  Solutions  to  (3.  4.  11)  give  a  good  approxi¬ 
mation  to  saddle  point  location  for  arbitrary  L/L  and  7»1.  The  two 

P 

approximate  methods  have  a  common  region  of  validity  when  T»1  and 


L/L  «i. 

P 

The  decay  regions  relative  to  the  saddle  points  have  been  investigated 
$ 

in  Appendix  C.  They  are  shown  in  Figs.  3.3  and  3.  4.  .  The  shaded  areas 
correspond  to  decay  regions  of  the  integrand,  fhe  integration  path,  C, 
will  now  be  deformed  into  the  path  C  '  as  shown  in  Figs.  3.  3  and  3.  4.  The 
deformed  path  has  been  chosen  to  go  through  the  two  saddle  points  which 
lie  on  the  integration  path,  and  everywhere  else  to  lie  in  a  decay  region. 

In  Appendix  B,  an  investigation  of  the  singularities  of  the  integrand  t  at  lie 
between  the  original  and  deformed  path  is  made.  There,  it  is  foun^  that 
those  singularities  give  rise  to  residue  contributions  which  are,  at  most, 
an  exponentially  small  order,  and  can  be  neglected.  As  a  result,  the 
integrals  over  the  two  paths  are  asymptotically  equivalent. 

Siiice  the  deformed  path,  C ' ,  lies  in  a  decay  region  except  at 
the  two  saddle  points,  the  integral 

The  saddle  points  and  decaying  regions  are  for  the  dominant  term  of  the 
asymptotic  expansion  of  the  integrand.  Other  terms  in  the  integrand's 
asymptotic  expansion  will  not  affect  our  results. 
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er“  4m  J 


1  -  P1  iko[P2(*  ’  *')  +  I“  +  £  *(pl ■  A*T>] 


C'  "» 


(3.4.  12) 


U  asymptotically  equal  to  tha  two  saddle  point  contributions.  The  method 

of  steepest  descent  can  be  used  to  evaluate  this  integral  if  we  generalise 

tiie  method  to  include  saddle  points  wUih  depend  on  the  large  parameter  k 

o 

The  dominant  term  of  the  asymptotic  expansion  of  the  integral 

»>i  lk„»<vV 


I  =  f  —  e  "  *  "  dp, 

J  PPj 


S.D.  o. 


(3.4.  13) 


where 

•(Pl-kJ*  P2(s  +  s/)  +  px+  *(p1,(A)*t)  (3.4.14) 

o 

and  the  path  of  integration  is  the  steepest  descent  path,  S.  D.  P.  is  given  by 


where 


I  -  i  k 


G(0)  e 


ik  *(p, 
o  Is 


ko> 


(3.4.  15) 


G(°)  =  £>U/<P,P„>][>"/  l*'(P„>l]i- 


dkin/4 


(3.4.  16) 


Here  we  denoted  the  second  derivative  of  #(p,k)  with  respect  to  p  as  $'(p.  ) 

o  1  Is 

and  have  assumed  there  is  only  one  saddle  point,  p^ ,  on  the  S.  D.  P.  The 
choice  of  the  +  or  -  sign  depends  on  the  direction  of  the  path  and  sign  of 
|*(p^).  The  above  evaluation  is  good  only  if  no  singularities  of  the  integrand 
approach  the  saddle  point  ask^  becomes  large.  For  the  saddle  points  treated 
in  this  paper  this  Is  the  case  as  is  shown  in  Appendix  B. 

The  asymptotic  evaluation  of  (-3.4. 12)  i».  now  performed.  For  L  going  to 
aero,  we  can  identify  the  contribution  of  the  lower  saddle  point  as  the  lateral 
wave,  while  that  of  the  upper  saddle  point  is  the  reflected  wave.  This 
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identiflcation  will  be  retained  for  arbitrary  values  of  L,  i.  e. ; 


where 


el* 


E  ~  E,  +  E 

yr  L  r 


ilr  fn  (eIc'Ud  *  +k  * 1  *tv  .  in/4 

oL*ai-*  *1*  o  '-la.-  •  j 

plLe 


(3.4.  17) 


2(2">PLP2L 


P2L  PL  ° 


(3.4.18) 


and 


E  = 
r 


‘k  fp2r<I  +  ''l  +  prx  V*(»,r^T)>  1  ivli 

pl,«  L _ 

2«2">*prp2r  [t' 'I  -'.W 


2r  rr 


(3.  4.  19) 

with  ♦ "  being  the  second  derivative  of  ♦  with  respect  to  p^.  The  positive  sign 

has  been  chosen  for  p,  in  (  3.  4.  6)  and 

lr 


-  J1) 

P1L  P1L 


P  =P(U 
Flr  Flr 


1  2  \ 
*  P2L  =  +  P1L)  ’  PL  =  (€1  “  P1L) 


2  i 


2  i 


2r 


=  ^  Plr  )  .  Pr  =  («j  '  Plr  )  • 


Since  *'(p  ,  (Aj  T)/k  ~  O  (k  )  for  all  t,  -;t  can  be  neglected  in  (3.4.18)  and 

1  r  o  o 

(3.  4.19)  if  only  the  dominant  term  of  the  asymptotic  approximation  is  desired. 

It  is  interesting  to  see  how  these  integrals  reduce  when  t  is  in  the  range 
where  the  saddle  points  can  be  obtained  by  the  method  of  successive  approxi¬ 
mations,  i.  e. ,  (3.  4.  7)  and  (3.  4.  8)  .  The  contribution  from  the  reflected  wave 
part  of  the  integral  is 


i 
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,  lko[P2r  <*+*' )  +  px  -Kk  1  *  (pj  r,  (A)*t  )  |f  in/4 
I  © 


E  - 


P2r  Pr 


(3.  4.  20) 


plr  *  *(x  2a-<*+*')2  . 


This  is  the  same  contribution  one  would  obtain  from  a  sharply  bounded  half 

space  problem  except  for  the  additional  phase  term.  When  the  method  of 

successive  approximations  is  applied  to  E  ,  we  find 

L* 

II*  r(A)*T  1  +  <  V  **+ k„' I*t<a>*T ]j+i3"/4 

EL - 1  (3.4.21) 

<2">MkoLJ 


where 


1  J  [  T(A)fr  (v 


/  2+ 1  /  2) 


tan(n  v/2){  (3.4.22) 


a[(A)*t]=  (A/4)*(d*/dPl] 


nT2Asec2(nv/2)r2(v/2-H)i^(v/2-H/2) 
Pj«0  4tan2(nv/2)l'4(v/2+l)+T2Ar4(v/2+l/2) 


*  (3.4.23) 


Equation  (3.4.  21)  in  the  limit  of  small  t  reduces  to 


ikJ^(s+z/)+  +  i  3rr/4 


(2n)  A  (k  L  ) 
o  p 


»  T  «1  .  (3.  4.  24) 


This  is  the  same  as  a  lateral  wave  on  a  dielectric  half  space.  On  the  other 
hand,  if  we  assume  t  is  large  enough,  the  asymptotic  approximation  for  the 


Lateral  Wave  Amplitude  vs.  Normalised  Transition  Thickness 


V  hickness 
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Gamma  functions  can  be  used  and  we  obtain  the  geometric  optical  result. 
It  is 


/4L  ikji^(z+z')  +  (s^x  H  TtL  (A'^  iT7,'4 

Er - 1 - 1  c  1 

2(2A)*L  (k  L  )* 


.  (3.4.25) 


pop 


The  lateral  wave  amplitude,  A^(a)^t  j  and  phase,  f£(A)^T  J,  have  been 
plotted  in  Figs.  3.5  and  3.6  respectively.  The  straight  line  (shown  dotted) 
appearing  in  both  graphs  is  the  approximate  solution  when  the  transition  is 
slowly  varying  with  respect  to  wavelength.  The  dash-dot  line  is  the  approximate 
solution  when  the  transition  region  is  small  compared  with  wavelength.  In  this 
case  the  reflection  coefficient  has  been  approximated  by  the  first  few  terms 
of  the  Taylor  series  expansion  of  T  about  T  =  0,  i.e, 

.2 


f  .1*2!  ♦  2SS.T.2 


Vi 


plp2  (VPl)Y  *  L  p2  2  ,  ,  ,w 

2  T  ,  t  «l  (3.4.26) 

(P2+  Pi) 


where 


v  =  2Aln  2  +  0(Pj  )  . 

The  lateral  wave  contribution  to  each  term  in  (3.  4.  26)  was  found  and  then  these 
three  lateral  wave  contributions  were  r-»*t  in  the  form  of  (3.4.21)  to  give  an 
approximate  lateral  wave  amplitude  a£(A)^t  J  and  phase  (£(A)^t  j  . 


3.5 


CONCLUSION 


For  purposes  of  investigating  the  high  frequency  reflected  field  from 
the  symmetrical  Epstein  layer  the  homogeneous  half-space,  z  ^0,  divides 
naturally  into  two  regions:  first,  the  region  between  the  z  axis  and  the 
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caustic  (Fig.  3.2)  where  the  reflected  field  is  exponentially  small  and  , 
second,  the  region  to  the  right  of  the  caustic  where  there  are  two  major 
contirubtions  to  the  field  at  each  point.  The  first  of  these  contributions  is 
the  reflected  ray  contribution  ,  (3.4.  20)  ,  whose  amplitude  variation  with 
frequency  is  kQ  .  The  frequency  dependence  does  not  change  as  t  is  varied. 
The  second  field  contribution  to  the  right  of  the  caustic  (3.4.  21)  has  an  ampli¬ 
tude  which  is  dependent  on  t  .  For  small  t  the  contribution  approximates  a 

-3/2 

lateral  wave  on  a  sh?*o  interface  whose  amplitude  dependence  is  (k  L  )  . 

o  p 

The  caustic  in  Fig.  3.  2  provides  the  natural  boundary  for  this  contribution 
since  for  small  t  it  coincides  with  the  critically  reflected  ray  in  th»j  sharp 
interface  probl-ar-  When  t  becomes  large,  this  second  field  contribution 
changes  chaj?.cter  from  a  lateral  ray  to  a  reflected  ray  contribution.  The 
amplitude  dependence  now  has  a  (k^L^)^  variation.  Approximate  forms  of  the 
seco.:d  contribution  are  given  in  (3. 4. 24)  and  (3.  4.  25)  for  small  and  large  t 
respectively.  An  examination  of  these  two  expressions  shows  that  the  lateral 


distance  dependence  L 


remains  unchanged  as  t  is  varied. 


The  two  field  contributions  to  the  right  of  the  caustic  correspond  to  the 
rays  through  each  point  in  this  region.  In  the  homogeneous  region  the  ray? 
that  make  up  the  caustic  are  the  rays  that  correspond  to  the  contribution  E.  , 
given  in  (3.4.  21)  .  The  rays  that  form  the  caustic  which  approaches  the  z  *  0 
line  asymptotically  correspond  to  the  contribution  ,  or  (3.  4.  20)  .  Identi¬ 
fication  of  contributions  corresponding  to  rays  that  pass  near  the  focus  is  more 
difficult  and  must  be  considered  when  L  is  small. 

In  Chapters  1  and  2  we  found  that  a  lateral  wave  was  excited  on  the  layer 

for  all  values  of  t  .  This  lateral  wave  changed  its  amplitude  dependence  on 

frequency  from  (kL  )  ^  to  (k  L  )  ^  for  the  linear  layer  and  from  (k  L  ) 

®  P  ®  P  ®  P 

to  (k  L  )~  for  the  parabolic  layer  as  t  went  from  small  to  large  values.  This 
o  p 

contribution,  however,  always  remained  a  diffraction  effect,  not  predictable  by 
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dassical  geometric  -  optics.  In  the  Epstein  transition  treated  in  this  chapter 
the  lateral  wave  is  excited  for  small  r  but  as  r  becomes  large  this  contribu¬ 
tion  becomes  a  reflected  wave,  which  can  be  predictable  solely  on  the  basis  of 
geometrical  optics.  This  shows  that  the  added  continuity  of  the  dielectric  trans¬ 
ition  has  a  marked  effect  on  the  behavior  of  the  lateral  wave.  The  transition 
chosen  in  this  chapter  has  the  drawback  that  it  has  an  infinite  width.  This  has 
the  effect  that  a  ray  which  emerges  from  the  inhomogeneous  medium,  having 
traveled  a  large  distance  along  the  interface,  penetrates  deeply  into  the  strati¬ 
fied  medium.  It  would  be  of  interest  to  investigate  transitions  with  a  finite  width 
but  more  continuity  than  the  layers  treated  in  Chapters  1  and  2  . 
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CHAPTER  4. 

THE  DOUBLE  EXPONENTIAL  TRANSITION  LAYER 


4. 1  Introduction 

In  tlu  last  chapter  the  reflected  fields  from  half  of  a  symmetrical 
Epstein  layer  were  investigated  for  arbitrary  layer  thickness.  When  the 
layer  thickness  was  large  in  comparison  to  wavelength,  the  reflected  field 
was  composed  of  geometrical -optic  rays  forming  a  two  branched  caustic 
which  met  in  a  cusp.  As  the  layer  thickness  decreased,  a  portion  of  these 
reflected  waves  changed  character  and  formed  a  lateral  wave. 

As  the  next  step  in  our  investigation  of  the  reflected  field  from 

transition  layers,  it  would  be  ideal  for  us  to  study  a  dielectric  profile  that 

is  completely  continuous,  i.  e. ,  one  having  all  derivatives  of  e(z)  continuous 

(4) 

for  «•<*<*.  The  Epstein  v  transition  layer  is  an  example  of  such  a  profile 
but  the  asymptotic  properties  of  the  wave  functions  have  not  been  investi¬ 
gated  in  enough  detail  to  make  the  problem  manageable  . 

Since  this  is  not  possible  at  present,  we  will  study  the  double  exponential 
medium  instead.  This  transition  is  composed  of  an  exponential  function  which 
approaches  as  a-*-*  for  z<0  and  an  exponential  function  which  approaches 
1  as  s-***  for  z>0  .  At  z  =  0  the  profile  and  its  derivative  are  continuous.  The 
double  exponential  profile  i3  similar  to  the  Epstein  profile  since  they  are  both 
composed  of  inhomogeneous  medium  for  all  z.  The  two  are  different  since 
the  double  exponential  medium  is  not  completely  continuous  for  all  z.  In  fact, 
the  double  exponential  medium  has  a  discontinuous  second  derivative  of  e(z); 
this  was  the  case  for  the  symmetrical  Epstein  layer  considered  in  Chapter  3. 

In  this  chapter  we  shall:  first,  evaluate  the  reflected  fields  when  the 
transition  is  thick  compared  to  wavelength,  and  second,  evaluate  the  reflected 
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fields  for  arbitrary  transition  thickness  where  the  source  and  observation 
points  are  many  transition  thicknesses  from  the  layer.  We  will  determine 
if  the  lateral  wave  in  this  medium  changes  its  character  as  observed  in 
Chapter  3  and  note  what  effects  occur  from  the  placement  of  the  source  in 
an  inhomogeneous  medium. 


4.  2  Statement  of  Problem  and  Formal  Integral  Representations 
Consider  a  dielectric  profile  given  by 


1  l-[4fc]  ,'2z,L 

z >0 

€(l)  °  j  ,  +2z/L 

(  +[A/2]e 

z<0 

A  =  1-€1 

(4.2.  1) 

where  an  electric  line  current  source  is  located  at  x  =  0  ,  z  *  z' ,  z'>0  and 
parallel  to  the  y  axis.  It  is  assumed,  as  before,  that  the  e  time 

dependence  is  understood.  Because  of  the  source  location  only  the  y 
component  of  the  electric  field,  E  ,  is  excited  and  it  obeys  the  equation 


[ 


k^'  e(z)  1  E  =  -  6  (x)  6(z-z/) 
o  J  y 


(4.  2.  2) 


The  source  magnitude  has  been  adjusted  to  make  the  coefficient  of  the  delta 
function  -1  and  is  given  by  Eq.  (1.  2.  2).  Since  the  medium  is  uniform  in  the 
x  direction,  it  becomes  convenient  to  introduce  the  transformation 

k  + »  ik  p 

E=  J  *(Z»P)  e  °  dp  (4.2.3) 

^  a  CO 

into  Eq.  (4.  2.  2) .  The  resulting  equation  for  ty(z#  p)  is  given  by 


[■^2  +  k^(M*)-p2)]  t(*#p)  *"  6(z-z') 


(4.  2.  4) 


where  \)(z,p)  must  obey  the  radiation  condition  as  z  .  The  above 
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equation  is  a  one  dimensional  Green's  function  problem.  The  Green's 

(27) 

function  %(z,  p)  is  given  by  Friedman  as 

^(z.p)  =  W(*2,*1)  (4.2.5) 

where  ?2  and  ^  are  homogeneous  solutions  of  Eq.  (4.  2.4)  which  satisfy 
the  radiation  condition  at  plus  and  minus  infinity  uespectively.  The  notation 
z<  means  that  z^  is  equal  to  the  lesser  of  z  or  z'  while  means  that  z>  is 
greater  of  z  or  z  .  The  denominator  of  Eq.  (4.  2.  5)  is  the  Wronskian  of  * 
and  (fj  which  is  defined  by 

1 

W(*2,ti)  =  *?(z)  d^(z)/dz-  ^(z)  dt2(z)/dz  ,  (4.2.6) 


We  will  now  restrict  ourselves  to  finding  f(z,p)  in  the  region  z>0. 

Upon  introducing  the  dielectric  profile  €(z)  for  z>0  into  Eq.  (4.  2.4)  and 
l?  8) 

using  Abromowitz1  '  ,we  find  that  the  homogeneous  solution  iji  (z)  is  given  by 

£> 

♦2(z)  =  Jv  (-iXe"z/Ij  )  ,  z>  0  (4.2.7) 


where  Jv(z)  is  the  Bessel  function  of  order  v  and  argument  z  .  The  symbols 
v2  and  X  are  defined  by 


(4.2.8) 


The  solution  ty^(z)  is  given  by 


V*)  = 


AJ  (Xez/L) 
V1 


,  z  <  0 


(4.2.  9) 


J  (-iXe"z/L)  +  rj  (-iXe'z/L)  ,  z>0 
'2  2 


where  A  is  constant  independent  of  z 


.  Vj*-  iTy.1  -  pz 


and 


-95- 


J  (-iX)  J'  {>.)  -  ij  (X)  j'  (-iX) 
2  1  1  2 
J  (-iX)  J'  (X)-  ij  (X)  J'  (-iX) 

2  1  1  2 


(4.  2.  10) 


The  prime  over  the  Bessel  function  indicates  differentiation  with  respect 
to  the  argument.  The  solution  (XeZ^)  satisfies  the  radiation  condition 

1  j  j  /  7 

as  z  oo  while  the  functions  {-iXe”z/  )  and  J_v  (-iXe"z/  )  are  two 
independent  homogeneous  solutions  of  Eq.  (4.  2.  4)  for  z  >0  .  A  linear  comb¬ 
ination  of  these  solutions  is  used  for  ^(z)  in  such  a  way  as  to  make  ^(z) 
and  d\Jij(z)/dz  continuous  at  z  =  0.  The  Wronskian  W(^i  ,  )  is  found  to  be 


\^2»  =  (2/ttL  )  sin  tt 


(4.  2.  11) 


•4_(29) 


by  using  Abramowitz'  1 .  By  substituting  ^(z)  ,  ^(z)  and  W *n 
Eq.  (4.  2.  5)  and  then  using  the  resulting  expression  for  i|f(z,  p)  in  Eq.  (4.  2.  3), 
the  formal  integral  representation  for  E^  is  obtained.  It  is 


ik  px 


Ey  =  4  [  ^sin(7TV2)^  -iXexp(-z</L)  Jv  j^-iX  expt-z^/Ljje  °  dp 


(4.2.  12) 


ik  px 


t  '■  r  r  -1  r  1 

+  4  J  iinNlnj"  Jv  [-awcpt-^/MJ  J  [-aexpl-z^Dje  dp  _ 


For  these  integrals  to  be  completely  defined,  the  square  roots 


must  be  properly  defined  on  a  four  sheeted  Riemann  surface  which  is  shown 
in  Fig.  1.4  .  The  multivalued  character  of  the  Bessel  function  need  not 
concern  us  since  the  integration  variable  is  not  contained  in  the  argument 
of  any  of  the  Bessel  functions  in  Eq.  (4.  2.  12)  . 


r 
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T  -0  , 

Eq.  (4 

.2.  12)  reduces  to 

1 

+  ® 
r 

ik«  (P? \z~z'  1  +  Px) 

e  dp 

4TTi 

"  00 

P2 

1 

+  ® 

1* 

ik  [p  (z+z')+px] 

~  O  £ 

1  e 

-  P2_P1 

4TTi 

L 

dP , 

f2 

'  p^+Pj 

(4.  2.  13) 


We  see  from  the  last  equation  that  as  t  -0,  or  as  the  layer  becomes  thin 
compard  to  wave  length,  the  formal  integral  representation,  given  in  Eq. 
(4.2. 12),  reduces  to  the  half  space  integral  representation. 


As  in  Chapter  1,  the  branch  points  at  p  =  are  now  removed  by 

means  of  the  transformation 

/ - T 

P  =  v  ej-Pl 


=  Jb  +  p? 


This  transforms  the  integral  in  Eq.  (4.  2.  12)  to 


T  (•  rl  r  r  -  ikQPX 

Ey  S  *  4  J  p sin(TTv  )  J- v  [-iXe*p(-.>/L)je  dPj 

c  z  2  z 

n  F  ’!»  (4.  2.  14) 

T  r  P1  r  r  11 -  n  lk«px 

-  4  J  psin(TTV0 )  JVo[-^exp(-z</L)JJvJ-iXexp(-z>/L)j 

c 


r  r 


e  dp  . 
1 


where  the  multivalued  functions  p.j^are  defined  on  a  four  sheeted  Riemann 
surface.  The  first  two  sheets  of  this  surface  are  shown  in  Fig.  1„  5  and  1.6. 

For  reasons  that  will  become  obvious  at  a  later  point,  it  ‘\ould  be 
advantageous  to  obtain  a  different  form  of  Eq.  (4.  2.  14).  This  can  be  done  by 
using  the  relation 

J_v(z)  =  eWi  Jv(z)  -  i sin(TTv)  H^2)(z)  (4.  2.  15) 


A. 
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in  Eq.  (4.  2. 14).  The  resulting  integral  representation  is 

Ey  =  ^rj  “T  H^^-iXenpt-.^L)]  Jv  [-iXexp^/L)].  °  dPj 

q  2  2 

(4.2.16) 

Pi  lk  px 

+  ~4~J  Jv  ['4XeatP(-t</L']  Jv  ["lXeXp<"*>/L)]  *  °  dpl 

g  2  2 

where  H^C-iXJj;  (X)  -Uv  (X)  (-iX) 

**  2  1  12 

r  “  Jv  (X)  -i^  (X)jJ  (-iX)  '  (4* 2‘ 17) 

2  1  i  2 


4.  3  Asymptotic  Evaluation  of  Integral  Representation  for  Large  keL 

4.3.1.  General  Gonslderatiopa 

Since  the  integral  representations  for  S  which  have  been  presented  in 
the  last  section  are  too  complicated  to  be  integrated  directly,  approximate 
procedures  must  be  used  to  simplify  the  integrals  further.  If  we  assume 
that  the  medium  is  slowly  varying  with  respect  to  a  wave-length  and  that  the 
observation  point  is  not  near  the  source,  we  can  asymptotically  approximate 
the  integral  representation.  To  be  more  specific  :  if  the  medium  is  slowly 
varying,  i.  e. ,  kQL»  1,  the  uniform  asymptotic  approximations  to  the  Bessel 
functions  can  be  used  to  simplify  the  integrand  of  Eqs.  (4.  2. 14)  and  (4.  2.  16). 
By  using  this  simplified  form  of  the  integrand,  the  integrals  in  Eqs.  (4.  2. 14) 
and  (4.  2. 16)  can  by  asymptotically  approximated  by  the  method  of  steepest 
descents  for  large  kQr  where  r  =  */x^  +  . 

Before  proceeding  with  the  asymptotic  evaluation  of  the  exact  solution, 
the  methods  of  geometric -optics  will  be  applied  to  .the  problem  in  the  next 
section,  4.  3.  2.  This  will  allow  us  to  represent  the  field  as  a  sum  of  ray 
contributions.  The  ray  family  will  be  investigated  in  detail  with  its  assoc" 
iated  caustics  and  foci.  A  comparison  of  the  field  contributions,  obtained 


from  the  asymptotic  evaluation,  will  be  performed  in  section  4.  3.  3. ,  and 
the  method  of  geometric -optics  will  clearly  distinguish  the  geometric -optic 
from  the  diffraction  effects. 


As  stated  in  the  last  section,  the  methods  of  geometric -optics  will  now 
be  applied  to  the  double  exponential  medium.  We  will  assume  L » 1  sin.-< 

the  method  requires  the  medium  to  be  slowly  varying  with  respect  to  wave¬ 
length.  By  using  this  method,  the  field  can  be  characterized  by  a  family  o i 
rays  emanating  from  the  source.  In  this  section  the  structure  of  this  ray 
family  will  be  studied  including  its  associated  caustics  and  foci.  Also  a 
detailed  investigation  of  this  ray  family  will  be  made  a»  L-*0  while  keeping 

k  L  large.  This  limiting  ray  family  is  important  since  a  comparison,  with 
o 

the  ray  family  for  a  sharp  interface  problem,  gives  important  clues  to  the 
changing  nature  of  the  lateral  wave  a  L  increases. 


The  ray  trajectories  can  be  obtained  by  an  integration  of  the  ray 
equation  i.  e. , 


(4.  3.  1 ) 


where  the  plus  and  minus  signs  are  used  for  rays  with  positive  and 
negative  slope  respectively.  The  symbol  p,  used  above,  is  the  ray  para¬ 
meter  and  is  related  to  the  initial  ray  angle  9  by  p  =  J e(z')  sin8  ;  the  angle 

o  o 

9  is  measured  between  the  tangent  to  the  ray  emerging  from  the  source  and 
o 

the  line  z  =  z' ,  x>0,  and  is  considered  positive  in  the  clockwise  direction. 


The  rays  divide  naturally  into  three  types  :  direct,  reflected  and  trans¬ 
mitted  rays.  The  direct  rays  are  restricted  to  three  regions  of  the  x  -  z 
plane  as  shown  by  the  Roman  numerals  in  Fig.  4. 1.  Region  I  is  bounded  by 
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a  portion  of  the  z  axis  and  the  glancing  ray,  i.  e. ,  the  ray  emitted  from  the 
source  at  0.  The  name,  ''glancing  ray",  stems  from  the  fact  that  the  ray 
just  glances  along  the  interface  in  the  limiting  case  of  an  abrupt  interface 
(L-*0).  The  glancing  ray's  trajectory  is  given  by 


x/L 


(1-&/2  e~2z// 
./ATT  e‘Z  /L 


* 


cosh 


-If  (z-z')f  LA 

\B  J 


(4.  3.  2) 


Region  II  is  bounded  by  the  locus  of  turning  points  and  the  critical  ray.  The 

equations,  describing  the  locus  of  turning  points,  are  given  in  Eqs.  (4.  3.  4) 

and  (4.  3.  8)  and  will  be  discussed  in  detail  at  a  later  point.  The  critical  ray 

is  emitted  from  the  source  at  the  angle  6  where  6  is  the  critical  angle  given 
_j  ^  _  c  c 

by  6  =  sin  Jc  /e(  z)  .  This  ray  is  found  by  integrating  Eq.  (4.  3.  1)  with 
cl 

p  =  */"e^  •  I*8  ray  trajectory  is 


x/L  = 


£cosh  1  {^-JTe7  -  cosh  1 ‘^^eZ^ J  ,  z>0 


(4.  3.  3) 


“  £cosh  *  ^  v^eZ  +  «/Ze  Z^  -  JT  -cosh  *  J  ,  z<0  # 


The  direct  rays  are  not  the  only  rays  existing  in  regions  I  and  II.  They  are, 
however,  the  only  rays  emitted  directly  by  the  source  into  regions  I  and  II  . 
Region  III  is  bounded  by  a  portion  of  the  x  axis,  the  zaxis  and  the  critical  ray. 

Direct  rays  emitted  in  the  angular  sector  0<9<6c  are  called  returning 
rays  after  they  pass  through  the  locus  of  turning  points.  These  returning 
rays  have  z  coordinates  that  tend  toward  plus  infinity  as  the  observation 
point  becomes  far  from  the  source.  Direct  rays  which  have  been  emitted 
in  the  angular  sector  9^<  9  <  tt / 2  and  have  passed  through  the  interface 
(z  =  0)  are  called  transmitted  rays.  These  rays  have  z  coordinates  that  tend 
toward  minus  infinity  when  the  observation  point  is  far  from  the  source. 


Because  of  the  monotonic  character  of  the  dielectric  profile,  the  direct 
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rays  never  cross  one  another;  the  sarr?  is  true  for  transmitted  rays.  As  a 
result,  these  two  ray  types  can  never  generate  caustics  or  foci,  therefore 
making  further  analysis  of  these  rays  unnecessary.  The  returning  rays,  on 
the  other  hand,  do  cross  one  another  and  must  be  studied  further. 

The  returning  rays  can  be  subdivided  into  two  ray  types:  those  rays  with 
turning  points  in  the  upper  half  space  and  those  rays  with  turning  points  in  the 
lower  half  space.  The  two  ray  types  are  divided  by  the  tangent  ray  which  has 
an  initial  angle  donoted  by  9^  as  shown  in  Fig.  4.  1. 


b.  Turning  Points 

Before  considering  the  caustics  formed  by  the  returning  rays,  an 
investigation  of  the  locus  of  turning  points  will  be  useful.  The  locus  points, 
zt>0,  are  found  by  integrating  Eq.  (4.  3.  1)  and  noting  that  p  =  J efz^)  .  The 
resultant  equation  is 

r  -2z  /Li|  z  /L  r  (z'-z  )/Li 

x^/L  =  ./2/A  U-(A/2)e  J  e  cosh  Le  J  ,  z^>0  .  (4.3.4) 

Some  interesting  leatures  of  this  curve,  as  illustrated  in  Fig.  4.  1,  are: 
first,  the  locus  intersects  the  source  point  and  has  a  zero  slope  at  this  point; 
second,  the  curve  has  one  maximum  (dx^/dz^  =  0}  when 

[  -2z7l_.  -1  z'/L.  .. 

J  1-e  cosh  e  ->  l-A/2  (4.3.5) 

and  is  monotonically  decreasing  when  inequality  (4.  3.  5)  is  not  fulfilled. 

The  location  of  the  maximum  will  be  denoted  by  (x_  z.  .)  . 

tM  tM 

If  {z'-z  )/L«  1  or  (z*-z  )/L»l  ,  Eq.  (4.  3.  4)  can  be  simplified.  By 

t  V 

use  of  the  former  approximation,  the  locus  of  the  turning  point  reduces  to 
x^  2  J (z#-z^)  L / A  eZ  ^ 


(4.  3.6) 
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This  describes  the  section  of  the  locus  of  turning  points  near  the  z  =  z' 
line.  When  (zJ -z  )/L»l  ,  on  the  other  hand,  we  obtain  the  approximate 
form  of  the  locus  of  turning  points  in  the  region  z7>z^>o  but  not  including 
those  points  close  to  z  .  This  approximation  is 

-2z  /L  £  z  /L 

=  J2JK  |^1- (A/2)  e  J  (z'-z^)  e  .  (4.3.7) 


If  we  now  assume  that  I-  is  small,  we  see  from  inequality  (4.  3.  5)  that 
the  locus  of  turning  points  has  a  maximum.  As  L-*0  the  x^  coordinate  of 
this  maximum  teuds  toward  infinity.  The  approximate  forms  of  the  locus 
above  and  below  the  maximum  are  described  by  Eqs.  (4.  3.  6)  and  (4.  3.  7) 
respectively.  As  L-0  we  see  from  Eq.  (4.  3.  6)  that  the  portion  of  the  locus 

above  z  tends  toward  the  straight  line  z-z'  with  x>0,  while  we  see  from 

CM 

Eq.  (4.  3.  7)  that  the  p>ortion  of  the  locus  below  z^  tends  toward  the  straight 
line  z*0  with  J"2jh  (l-A/2)^  z'<x<®  .  A  sketch  of  the  locus  for  small  L  is 
shown  in  Fig.  4.  2  . 


The  turning  points  when  z^<0,  are  obtained  in  a  similar  manner  as  those 
for  zt>0,  i.  e.  ,  by  integrating  Eq.  (4.  3.  1 )  but  now  z ^  must  be  considered  as 
negative.  The  resultant  equation  is 

x^/L  =  (p/p2>  {cosh  1^<y2^P2eZ  ^LJ-cosh  J 


(p/  |p2  |)  cosh"1  £  |pj  |  JVh\ 


(4.  3.  8) 


with 


+z  /L 

lp1 I  =  /A/2  e 


(4.  3.  9) 


As  would  be  expected  from  the  continuity  of  the  dielectric  constant 
and  its  derivative  at  z  =0,  the  locus  of  turning  points  and  its  derivative 
are  also  continuous  there.  As  we  approach  the  critical  angle,  i.e.  ,  |p^  |-*0, 
we  find  from  Eq.  (4.  3.  9)  that  z  -*•  00  .  Using  this  limit  in  Eq.  (4.  3.  8)  showr 
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that  the  locus  of  turning  points  approaches  the  critical  ray,  Eq.(4,  3.  3), 
as  z  “•  The  locus  in  the  region  z^<0  has  one  extremum  when  inequality 
(4.  3.  5)  is  satisfied,  and  it  has  no  extremum  when  it  is  not  satisfied. 

For  small  L  inequality  (4.  3.  5)  is  satisfied,  and  the  locus  of  turning 
points  has  one  minimum  which  is  given  asymptotically  by 

=  JTjEz  +  OfL^3)  ,  ~  C(LlnL1/3)  .  (4.3J0) 

tm  1  tm 

In  the  limit  of  L-»0,  we  see  x. _  -•  J c- .  /A  z '  and  z, _ -0  .  If  we  assume 

tm  1  tm 

|z^|/L»l,  the  portion  of  the  locus,  that  approaches  the  critical  ray,  can 
be  approximated  by 

UJ/L 

x  s»  J <c  /Az7  +  ti  Jt  /2A  Le  ,  L«  1.  (4.  3.  1 1 ) 

V  *  * 

If  we  let  |ztl  ~  Lln(8/L)  ,  then  x^  becomes 

~  Jl~jEz'  +  tt  ,/ e  1  / 2A  (4.3.12) 

If  0  is  varied,  all  values  of  x  in  the  interval  ^e./Az'  <x  <®  are  obtained 
with  the  exception  of  the  region  close  to  the  point  x  =  J e^/A  z'  .  In  the  limit 
of  small  L,  this  section  of  the  locus  of  turning  points  tends  towards  tie 
interface.  All  z  of  order  larger  than  Lln(8/L)  lead  to  x  which  approach 

W  V 

infinity  as  L-*0  . 

The  only  section  of  the  curve  not  investigated  as  yet,  is  the  section 
between  the  minimum  and  the  zt  =  0  crossing  point.  Since  |p^  |  is  non  zero 
from  Eq.  (4.  3.  9)  we  see  zt~*0  as  L~*0,  and  as  a  result,  this  section  of 
the  curve  also  approaches  the  interface  as  L-»0.  To  summarize,  all 
portions  of  the  locus  of  turning  points,  which  lie  in  the  finite  x  -  z  plane  as 
L-»0,  approach  the  interface  and  the  z-z  line. 
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c.  Returning  Rays,  z^>  0 

As  mentioned  previously,  those  rays,  that  are  turned  by  the  medium, 
cross  one  another  after  they  have  turned.  This  leads,  as  shall  be  seen,  to 
the  formation  of  a  caustic.  In  this  section  we  shall  investigate  the  section 
of  the  caustic  formed  by  those  returning  rays  with  turning  points,  z^  >0  . 

The  ray  equation  for  returning  rays  is  found  by  choosing  the  minus  sign 
in  Eq.  (4.  3.  1)  and  integrating  from  z'to  z^  and  then,  choosing  the  plus  sign 
and  integrating  from  z^  to  z  .  The  resultant  ray  equation  is 

x/L  =  (p/p2)  jcosh  1  j^p2  JTJK  eZ  ^Lj  +  cosh  eZ^LJ  j  (4.  3.  13) 

To  find  the  caustic,  the  constraint  equation  will  be  needed.  This  equation, 
found  by  taking  the  derivative  of  the  above  equation  with  respect  to  p,  is 

cosh  '^Py^/A  eZ  ^Lj  r  cosh  eZ^LJ  =  p2p2£l/p2(z')  +  l/p2(z)j 

(4.  3.  14) 

where 

p2 (z)  =  [p2  "  (A/2)  e  2z/LJ^  ,  =  lim  p2(z)  ,  (4.3.15) 

z  — 00 

Since  both  Eq.  (4.  3.  13)  and  Eq.  (4.  3.  14)  are  transcendental  in  p  ,  it  appears 

Ct 

impossible  to  eliminate  the  parameter  directly  between  the  two  equations  and 

obtain  the  caustic  directly.  However,  It  is  possible  to  obtain  an  asymptote  to 

i  z7l 

the  caustic  in  the  region  where  (z-z,)/L»l  and  p2(z  )  e  «1  .  By  using 
these  relations  to  simplify  the  constraint  equation,  a  relation  is  obtained 
between  z  and  p2(z')  •  When  this  relation  is  substituted  in  the  ray  equation, 
the  asymptote  obtained  is 

z«  x£(2/A)  e2z/L -l]^  +  z'  -  Lln2  ,  (4.3.16) 
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This  asymptote  is  a  straight  line  which  is  valid  when  z  >z*  . 

Another  general  feature  of  the  caustic  is  its  intersection  with  the 
line  z  =  z' .  It  is  shown  in  Appendix  D  that  the  ray  passing  through  an 
extremum  of  the  turning  point  locus  also  passes  through  the  caustic  at 
z  =  z# .  Thus,  if  the  rays  forming  a  caustic  have  a  locus  of  turning  points 
with  an  extremum,  the  caustic  must  cross  the  line  z-z'.  The  condition  that 
the  locus  of  turning  points  has  an  extremum  for  z  >0j  as  expressed  in 

*  I 

Eq.  (4.  3.  5)  ,  is 

f.  -2z'/L  -1  /  z7L\  ,  . 

•J  1-e  cosh  \e  )>l-A/2  (4.3.17) 

and,  if  satisfied,  the  section  of  the  caustic  for  z^>0  crosses  the  z  -z 
line  once. 

Because  of  the  complex  analyt.c  character  of  the  ray  and  constraint 
equations,  a  GE235  computer  was  used  to  plot  the  caustic  for  several 
different  transition  thicknesses,  L.  The  basic  procedure  for  doing  this  was 
to  pick  a  particular  value  of  p  ,  and  then  solve  the  constraint  equation  for 
z.  Since  this  equation  is  transcendental  in  z,  a  Newton  Raphson  method 
was  used  to  find  z.  The  complete  caustic  has  been  plotted  for  L=  10,  1, .  1 
and  is  shown  in  Figs.  4.  3,  4.  4  and  4.  5  respectively. 

For  L>=  iO,  1  the  caustic  consists  of  two  branches  meeting  in  a  cusp. 

The  cusp,  having  the  ray  with  z  =  0  passing  through  its  tip,  is  not  actually 

*  * 

a  regular  cusp  since  the  local  condition  is  not  satisfied  at  its  tip.  The 
cusped  nature  of  the  caustic  is  instead  formed  because  of  the  discontinuous 
character  of  the  dielectric  distribution  at  z  =  0.  The  lower  branch  of  this 
caustic,  having  z^  >0,  is  the  caustic  described  in  this  section. 


♦ 


The  focal  condition  states  that  the  second  partial  derivative  of  the  ray 
equation  with  respect  to  the  ray  parameter  must  be  zero  at  the  focal 
point. 


CAUSTIC 

LOCUS  OF  TURNING  POINTS 
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In  Fig.  4.  5  the  structure  of  the  caustic  changes  because  of  the 
appearance  of  two  new  cusps  :  one  below  the  x  axis  and  one  above.  Both 
of  these  new  cusps  satisfy  the  focal  condition  at  their  tip.  An  enlargement 
of  Fig.  4.  5  in  the  area  of  the  cusps  is  shown  in  Fig.  4.  6  .  We  note  f  he 
regular  cusp  below  the  x  axis  is  shown  clearly  in  this  figure,  while  the  two 
cusps  appearing  above  the  x  axis  occur  too  closo  to  one  another  to  be  distin¬ 
guishable.  To  rectify  this  situation  a  detail  of  Fig.  4.  6  has  been  plotted  in 
the  region  of  these  two  cusps.  This  detail  is  shown  in  Fig.  4.  7.  The  coord¬ 
inates  used  in  this  figure  are  RT  and  DL.  The  quantity  DL  is  the  distance 
measured  along  a  straight  line  passing  through  the  tip  of  the  irregular  cusp  ; 
DL  is  equal  to  zero  at  the  point  the  straight  line  passes  through  the  cusp's  tip. 

The  slope  of  this  line  is  the  same  as  the  slope  of  the  caustic  at  the  irregular 

* 

cusp's  tip  .  The  quantity  RT  is  measured  perpendicular  to  DL  as  shown  in 
Fig.  4.  7  .  The  coordinates  of  the  irregular  cusp's  tip,  (x^.,  z^.) ,  and  the  angle  , 
9,  that  the  DL  axis  makes  with  the  x  axis  are  given  at  the  bottom  of  Fig.  4.  7  . 

An  examination  of  Fig.  4.  7  shows  that  the  caustic  forms  a  bow  tie 
configuration  similar  to  that  observed  in  the  linear  layer  (Chapter  I).  The 
portion  of  the  caustic  tending  toward  the  lower  part  of  the  graph  is  composed 
of  rays  with  turning  points  z^<0  ,  and,  if  extended,  it  would  connect  to  the 
cusp  shown  in  Fig.  4.  6.  On  the  other  hand,  the  portion  of  the  caustic  that 
appears  to  coincide  with  the  DL  axis  corresponds  to  rays  with  turning  points 
z^>0,  if  extended,  it  would  produce  the  lower  branch  of  the  caustic  shown 
in  Fig.  4.  5  .  We  also  note  that  the  irregular  cusp  occurs  at  DL=  0  =  RT  as 
would  be  expected  from  the  definition  of  the  coordinates  . 


* 


The  slope  of  the  caustic  exits  at  the  cusp's  tip  and  is  unique. 


Detail  of  Caustic  for  Double  Exponential  Medium, 
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As  L  becomes  increasingly  small,  various  approximations  can  be 
made  which  enable  us  to  simplify  Eqs.  (4.  3.  13)  and  (4.  3.  14)  that  represent 
the  lower  branch  of  the  caustic.  First,  we  note  from  section  4.  3.  lb  that 
for  z^>0  and  L  small,  a  maximum  of  the  locus  of  turning  points  always 
exists  and  as  a  result,  the  caustic  crosses  the  line  z  =  z'  .  An  asymptotic 
approximation  to  the  caustic  can  be  obtained  below  this  line.  If  we  assume 
(z/-z)/L»l  and  (z7-*  )/L»l  ,  the  constraint  equation,  4.3.14,  reduces  to 


-2z  /L_ 


z  -  z  «*  L3  ;  1  -  (Ul)  e  4  j  /2(z'  -  z  )2, 
t  L  J  t 

By  using  the  above  approximation  in  the  ray  equation,  we  find 

2z  /L  i 

x  =  ^(2/A)  e  1  -lj*  {z'-z^  , 


(4.  3.  18) 


(4.  3.  19) 


An  examination  of  Eq.  (4.  3.  18)  as  L-0  shows  us  that  z-z^.  Equation 

4.  3.  19  divides  naturally  into  two  cases  :  first,  those  z  attL  which  lead  to 
1  ,  1 
(21  h-  1 )  z  <x^<®  while  z^  with  a  weaker  dependence  on  L  have  x^  coord¬ 
inates  which  tend  toward  infinity  as  L-0  .  With  the  above  information  and 
with  the  previous  knowledge  that  z  -0  for  finite  x  ,  we  conclude  that  this 
portion  of  the  caustic  coincides  with  the  z  =  0  line  for  (2/A  -1)E  z/<x<®. 

A  comparison  of  Figs.  4.  3,  4.  4  and  4.  5  will  show  how  this  limit  is  approached 
as  L- 0  . 

The  original  asymptote,  Eq.  (4.  3.  16),  still  holds  when  L  is  small. 

From  Eq.  (4.  3.  16)  we  see  that  any  z  >  z  has  a  corresponding  x  coordinate 
which  tends  towards  infinity  as  L-0  .  To  summarize  in  the  limit  of  L-*0  , 
a  portion  of  the  caustic  lies  along  the  x  axis,  while  other  points  on  the  caustic, 
corresponding  to  non- zero  z  values,  have  x  coordinates  that  lie  at  infinity. 

The  focal  condit  on  for  this  section  of  the  caustic  is 

3  Ll/p2<z/)  +  1/P2^Z0  +  p2[1/P2^,)  +  =  ° 


(4.  3.  20) 
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where  )  is  defined  in  Eq.  (4.  3.  15).  As  can  be  seen  this  equation  is 
never  satisfied  for  real  ray  parameters  and  thus  no  foci  exist  on  this 
section  of  the  caustic.  In  particular  this  equation  is  not  satisfied  at  z^  =  0 
which  corresponds  to  the  ray  that  passes  through  the  tip  of  the  irregular 
cusp. 


d.  Returning  Rays,  z^<0,  z>0 

The  rays  with  turning  points,  <0,  penetrate  :nto  the  medium  z>0  . 
Their  ray  equation,  as  in  previous  cases,  is  obtained  by  an  integration  of 
Eq.  (4.  3.  1).  The  ray  equation  is 


x/L  =  (p/p2)j-2  cosh  ,/2/A  p2  j  +  cosh  1  A  p^  eZ  ^  L] 


(4.  3. 21) 


cosh'1|^/27Ap2  eZ/L]  }+(2p/|Pl|)  cos  I/2M  ,  z^<  0  ,  z  >0  . 


These  rays  form  another  section  of  the  caustic  found  in  the  previous  section 
4.  3.  c.  The  constraint  equation  for  this  caustic, 

0  =  2  cosh  p2J  -  cosh  *£,/2/A  p^  eZ  ^LJ  -  cosh  p^  eZ^j 

+  (ZclP2/IP||3)co8"1['®^  Ip[I]  +(4P2P2/ IP]|2)[p2  W-3-22) 

+  p2^j[p2-!/W)e'2z7L]‘i  +  [p 2  -  (A/2)e'2*/Lj' 

was  obtained  by  taking  the  derivative  of  the  ray  equation  with  respect  to  p^  . 
Since  elimination  of  the  parameter,  p^  ,  between  Eqs.  (4.  3.  21)  and  (4.  3.  22) 
seems  impossible,  the  caustic  was  plotted  on  the  computer  by  a  similar 
procedure  to  that  described  in  the  previous  section. 
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The  upper  branch  of  the  caustic  shown  in  Figs.  4.  3  and  4.4  is  the 
section  of  the  caustic  whose  rays  have  z  <0  ,  i.  e.  ,  the  part  of  the  caustic 
described  in  this  section.  In  Figs.  4.  5  and  4.  6  ,  the  part  of  the  caustic, 
resulting  from  rays  with  zt  <0,  z>0,  is  more  complicated.  It  consists  of 
the  upper  branch,  starting  at  the  x  axis,  as  shown  in  Fig.  4.  5  and  the  section 
between  the  tip  of  the  irregular  cusp  and  the  x  axis,  as  shown  in  Figs.  4.  5 
and  4.  6.  Since  the  section  of  the  caustic  between  the  irregular  cusp  and  the 
x  axis  is  very  small,  an  enlargement  has  been  made  in  Fig.  4.  6  . 

As  can  be  seen  from  an  examination  of  Figs.  4.  3  -  4.  6  ,  the  general 
structure  of  the  caustic  changes  as  L  decreases  ;  however,  the  irregular 
cusp  always  remains.  The  analytic  behavior  of  this  irregular  cusp  can  be 
understood  better  by  noting;  first,  the  caustic  and  its  first  derivative,  dx/dz, 
are  continuous  at  the  tip  of  the  cusp,  z^  =  0  ;  second,  as  z^-0  from  the 
negative  side,  z  increases  for  all  values  of  A  and  L.  This  behavior  gives 
the  caustic  in  this  region  a  cusp-like  appearance. 

As  the  transition  thickness  gets  smaller,  the  irregular  cusp  approaches 
the  x  axis,  and  two  ad«.  itional  regular  cusps  appear.  One  of  these  cusps  is 
formed  from  rays  with  z^<0  ,  z  >0  which  are  being  considered  in  this  section. 
The  exact  location  of  the  focus  is  difficult  to  find  because  of  the  complicated 
analytical  character  of  the  focal  condition.  As  L  becomes  small,  the  focal 
condition  simplifies  and  an  approximate  location  of  the  focus  can  be  found. 


The  general  focal  condition  is 

|5.  -If 


0  =  (6*xa/  Ipx  |  )  cos  |p1|v^7Aj  + 


2  |  p.  I2p2  +  3e  (p£  -A/2) 


ipl  I4  (P^  "  A/2)' 


(4.  3.23) 


'3[I/p2(z')  +  •  P2[‘/P2<"')  + 
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wher*  p  ( z)  is  defined  in  Eq.  (4.  3.  15),  If  we  assume  L  is  small  and  z  -0  , 

2  t 

the  constraint  equation,  (4.  3.  22),  reduces  approximately  to 

z7L=(l-A/2)  [p2  -  (A/2)  e"2z/LJ  ^  (4.3.24) 

We  now  let  p  =  JEfZ  +  6  where  6«1  since  z  is  small  and  use  this  with 

b  t 

Eq.  (4.  3.  24)  in  the  focal  condition  Eq.  (4.  3.  23).  We  obtain 

6  =  32/a72(1-A/2)6  (z'f6  L6  (4.3.25) 

We  now  see  that  our  assumption  of  small  z^  was  valid  when  L  is  assumed 
to  be  small .  By  using  this  value  of  6  in  Eq.  (4.  3.  24)  we  find  the  approximate 
focus  location  is  given  by 

z  =  2(l-A/2)2(z,)"Z  L3  .  (4.  3.26) 

When  L  is  small,  an  approximate  equation  for  the  upper  branch  of  the 
caustic  can  be  obtained  by  using  z/L»l  and  z7L»l  in  Eqs.  (4.  3.  21)  and 
(4.  3.  22).  The  equation  obtained  is 

x/L  =  (z+z')/L  +  tt^3  €j  1/6  A"1/Z(z+z')1/3  L273  .  (4.  3.  27) 

The  use  of  Eq.  (4.  3.  27)  in  conjunction  with  Eq.  (4.  3.  23)  shows  that  there 
are  no  focal  points  on  the  upper  branch  of  the  caustic  for  small  L>. 

If  we  now  consider  the  limiting  case  as  L-*t  we  find  :  first,  the  section 
of  the  caustic  between  the  irregular  cusp  and  the  x  axis  shrinks  to  zero  as 
L-»0  ;  second,  the  upper  branch  of  the  caustic  approaches  the  line 
X  =  ^(zrz')  which  is  the  critically  reflected  ray  from  a  sharply  bounded 
interface  of  height  A  . 
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e.  Returning  Rays  z^<0,  z  <  0 

The  returning  rays  with  z<0  are  given  via  Eq.  (4.  3.  1)  by 

x/L=  (p/p2)  i  -cosh  1  !\/2/A  p2J  +  cosh  1  ./2/A  p2  eZ  ^Lj| 

(4.  3.  28) 

+  (p/!p1l)^cos’1[|p1  |  ^/Aj  +  cos"1^|p1  i  v£7Xe"z/Lj  j 

The  constraint  equation  is  obtained  by  taking  the  derivative  with  respect  to 
P2  of  Eq.  (4.  3.  28).  The  result  is 


0  =  cosh  1|_^7Ap2j  -  cosh  A£/2/A  p2  eZ  /-L,j  +  2p^p2  (p2  -  A/2)S  |pJ  | 


-l  r 


z#/L- 


Z  ,  2  .  |  I  -2 


+  (e1P2/lp1  I3)  -rcos_1[^  |  p  x  |  JtJI\  +  cos  ^  |Pj  |  JlTh  e"z/L  (4.  3.  29) 


,2  r  2 

+  P  P2  !  P 


2  -  (A/2) e  2Z  /L1  p2p3  |pj  |  2  [  (A/2)e2z/L  -  | Pj  I 


2]-* 

r  J  * 


A  computer  calculation  of  this  sector  of  the  caustic  was  made.  An  examina¬ 
tion  of  these  results  shown  in  Figs.  4.  3.  -  4.  6,  seems  to  indicate  that  if  a 
focal  point  does  not  exist  below  the  x  axis,  then  there  exists  no  caustic  below 
the  x  axis.  Applying  the  focal  condition  to  Eq.  (4.  3.  28) ,  gives  us  the  second 
constraint  equation.  This  equation,  in  conjunction  with  Eqs.  (4.  3.  28)  and 
(4.  3.  29) ,  enables  us  to  find  the  focal  points  .  Nothing  in  general  can  be  said 
because  of  the  complication  of  the  focal  condition. 

When  L  is  small,  Eqs.  (4.  3.  28) ,  (4.  3.  29)  and  the  focal  condition  simplify 
substantially,  enabling  us  to  locate  the  foci  if  they  exist.  We  find  that  one 
focus  exists  which  tends  towards  the  interface  as  L-*0,  i.  e.  ,  z  -»0,  x-*/e. ^/a  z' 
as  L-»0.  The  two  branches  of  this  cusp  tend  towards  the  interface  also,  as 
L-0  . 
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f.  Summary 

The  investigation  of  the  rays  emitted  from  a  point  source  in  a  double 
exponential  layer  has  led  us  to  conclude  :  first,  all  rays,  except  the 
returning  rays,  have  a  simple  character  ;  second,  the  returning  rays  form 
a  caustic  with  an  irregular  cusp  resulting  from  the  discontinuity  in  the  die¬ 
lectric  ;  third,  as  L  becomes  small,  two  additional  regular  cusps  appear 
one  below  the  x  axis  and  one  above  ;  fourth,  as  L-»0  ,  one  portion  of  the 
caustic  tends  toward  the  interface  while  the  other  portion  tends  toward  the 
critically  reflected  ray;  fifth,  the  focus  below  the  x  axis  tends  toward  the 
reflection  point  of  the  critically  reflected  ray;  and  sixth,  the  two  focii  above 
the  x  axis  tend  toward  the  turning  point  of  the  tangent  ray  as  L-0  . 


4.  3.  3  Asymptotic  Evaluation  of  the  Formal  Solution 

a.  Introduction 

The  asymptotic  evaluation  of  the  formal  solution  will  now  be  carried  out. 

As  has  been  mentioned  previously,  the  integrands  of  Eq.  (4.  2.  14)  or  Eq.  (4.  2.  16) 
will  be  asymptotically  evaluated  for  large  k^L  and  then  the  resulting  integrals 
will  be  asymptotically  approximated  for  large  k^r . 

b.  Asymptotic  Approximations  for  Large  k^L 


The  integrands  of  the  integrals  appearing  in  Eqs.  (4.  2.  14)  and  (4.  2.  16) 
contain  the  following  Bessel  functions  :  Jv  (-iXe  ,  j  v  (-iXe  )  t 

(-iXe  ^^ )  t  (- iX.  e  ,  Jv  (X)  and  the  derivatives  of  these  functions 

with  respect  to  their  arguments.  A  significant  simplification  of  the  integrands 
will  occur  if  the  above  functions  are  replaced  by  their  asymptotic  approxima¬ 
tions  for  k^L  large.  These  asymptotic  approximations  need  only  be  known  in 
a  strip-like  region  centered  around  the  integration  path,  C. 
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(3°)  (1) 

Olver  has  derived  uniform  asymptotic  expansions  for  J  (vz),  H'  '  (vz) 
(7.\  v  v 

and  '(vz)  and  for  their  derivatives  when  v  is  large  and  |argv|  <tt/2  . 
These  expansions  are  valid  in  the  whole  z  plane  when  a  branch  cu  intro¬ 
duced  along  a  curve  in  the  second  quadrant,  as  described  in  Append  £, 

Since  the  asymptotic  expansions  will  be  required  when  argv  =  -  ^  ,  Olver's 
results  have  been  extended  in  Appendix  E  to  include  this  end  point.  The 
extension  shows  that  the  same  expansions,  which  are  valid  when  |  argv|  <tt/2, 
are  also  valid  when  argv=  -  tt/2  .  These  uniform  asymptotic  expansions  are 


Jv  Ai(ivi»t 


2/3, 


V  “  1  " W2 
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where 


2  ,V2  v2  /  ..  ,  . 

1  >11  i  ‘  P 


«b/L 


(4.  3.30) 


(4.3.31) 


(4.  3.32) 


(4.  3.  33) 


(4.  3.  34) 


The  branch  cuts  which  define  the  correspondence  between  the  5  and 
planes  are  defined  in  Appendix  E.  By  using  the  definition  for  given  in 
Eq.  (4.  2.  8),  and  the  fact  that  the  above  asymptotic  approximations  are  valid 
for  Rov.,  >0,  we  find  that  these  approximations  are  valid  in  the  first  and 
third  quadrants  of  the  p^  plane  including  the  real  and  imaginary  axes. 
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In  Fig.  4.  8,  the  regions  of  validity  are  shown  by  shaded  regions  composed 
of  slanted  lines.  No  branch  cuts  due  to  the  square  root  p  =  ^”^+|p  are 
shown  in  this  figure  since  these  particular  Bessel  functions  are  independent 
of  p. 


The  Bessel  functions  ,  given  in  Eqs.  (4.  3,  30}  -  (4.  3.  33) ,  have  turning 
points  which  are  located  at  w^  ~  ±  l  .  The  turning  point  corresponding  to 

w  =  1  appears  in  two  places  on  the  top  sheet  of  the  p  plane  as  shown  in 

2 

Fig.  4.8.  These  two  points  are  located  at  p^  =  ±v.(z)  where  v^(z)=  (i-v^(z)) 
v  (z)  =  ./A/2  exp  {-z/X i)  .  The  two  points  corresponding  to  the  turning  point 
located  atw'-l,  are  located  on  the  bottom  sheet. 

-z/I 

",'o  ohi-ain  ih;.  asymptotic  approximation  of  J_V-. (-iX  e  ')  in  the  same 
region  of  validity  as  described  above,  we  will  use  the  connection  formula 

0  -z/L._  _  r  V2TTi  (l)  ..  -z/L.  ,  "V1  „(2h  -z/L.'' 

JP  {-i*e  |  =  le  ¥T  (-iXe  )rs  H  X-iXe  }  ' !  2  (4.3.33) 

2  V2  VZ 


where  the  asymptotic  approximations  given  by  Eq.  (4.  3.  31)  are  used  for 
A j  --/L 

H^y(-iX  e  ^  )  .  Since  it  can  be  shown  that  in  a  region  close  to  the  integra¬ 

tion  path,  C,  the  second  term  of  Eq.  (4.  3.  35)  is  exponentially  small.  The 

-z/L 

asymptotic  approximation  to  J_v?(-iXe  )  is  given  by 


_  ,  -z/L .  !v2|n+iTT/3  r4|v2!  ? 
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? -  ■ 
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(4.  3.  3fc) 


By  a  similar  procedure,  the  asymptotic  approximation  for  the  derivative 
of  JVz  can  be  obtained.  It  is 
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The  asymptotic  approximation  for  J..  (X)  and  its  derivative  can  be 

1 

found  by  using  Olver's  results  as  were  used  before. 

The  resultant  asymptotic  approximations  are 
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(4.  3.  38) 


{4.  3.  39) 


where  §  is  definec  as  in  Eq.  (4.  3.  34)  with  being  replaced  by  = 

-j£tZt ipj  .  The  region  of  validity  for  these  asymptotic  approximations  is 
the  upper  half  of  the  Pj  plane  including  the  Repj  axis  as  shown  by  the 
horizontal  lines  in  Fig.  4.  8 .  The  turning  point  of  Jv^  corresponding  to  w^  =  1  , 
is  located  on  the  Im  p,  axis  in  the  interval  0 <  Im  p^  <JK  while  the  turning 
point  corresponding  to  w  =  -  1  is  in  the  conjugate  position. 


The  asymptotic  approximation  to  the  Bessel  functions  Jv  J_v  and 
(a)  2  2 

Hv^  will  be  needed  in  a  narrow  strip  just  to  the  left  of  the  integration  path. 

C,  where  0<lmp  <JK.  By  making  use  of  the  formula  relating  different  types 

1  (31) 

of  Bessel  functions  given  by  Watson  we  can  show  that  the  asymptotic 
expansions,  which  are  valid  on  the  Imp^  axis,  are  also  valid  in  a  small  strip 
to  the  right  of  it.  These  asymptotic  expansions  are  not  valid  near  the  branch 
points,  Pj  =  ±i  JK  since  the  large  asymptotic  parameter,  tends  to  zero 
there . 


By  applying  these  same  ideas,  the  asymptoti;  approximations  for  Jv^  , 
J.v,  »  and  Jv  are  analytically  continued  to  a  narrow  strip  below  the  positive 
real  p^  axis  .  As  in  the  case  of  the  Bessel  functions  of  order  ±v^  ,  the  asymp¬ 
totic  approximations  of  Bessel  functions  of  order  v^,  are  not  valid  in  a  small 
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circular  region  about  the  origin  of  tne  plane. 

The  approximate  forms  of  the  Bessel  functions  previously  mentioned, 
can  be  simplified  further  by  using  the  asymptotic  formula  for  the  Airy 
functions  .  This  is  only  possible  when  the  integration  va  riable  p^  is  not 
loo  close  to  a  turning  point  of  the  particular  Bessel  function  under  consider¬ 
ation.  This  has  been  done  in  Tables  4.  1A  and  4.  IB  for  values  of  p^  along 
the  integration  pathC  .  These  formulae,  however,  can  be  analytically 
continued,  so  that  they  will  be  valid  in  a  narrow  strip,  centered  around  the 
integration  path.  The  asymptotic  approximations  presented  in  Table  4. 1A 

and  4.  IB  are  valid  in  the  neighborhood  of  p  =  0  and  p  =  0 .  Although  Olver's 

1  2  (32) 

uniform  expansions  break  down  at  these  points,  the  Debye  approximations  , 

which  utilize  the  large  parameter  \ ,  can  be  used  to  extend  the  range  of 
validity  of  the  approximations  appearing  in  the  tables  to  include  the  neighbor¬ 
hood  of  pj  =  0  and  =  G  . 

c.  Approximation  of  Formal  Solution 

The  asymptotic  evaluation  of  Eq.  (4.  2.  14)  or  its  alternate,  Eq.  (4.  2.  16)  , 
will  now  be  performed  with  the  a  of  the  asymptotic  approximations  obtained 
in  the  last  section.  Before  proceeding  with  this  evaluation,  we  would  like  to 
note  some  important  characteristics  of  the  integrands. 

The  integrands  as  expressed  in  section  4.2,  are  completely  defined. 

There  is  no  need  to  consider  the  multivalued  nature  of  the  cylinder  functions 
contained  in  these  integrands,  since  the  cylinder  functions  are  only  multi¬ 
valued  with  respect  to  their  arguments  which  are  a  constant  along  the  path  of 
integration.  These  cylinder  functions  have  turning  points  located  on  the  Im  p^ 
axis  athnp^  ±Vj(z<)  ,  ±  v^(z^)  ,  iv^O)  where  Vj(z< )  >Vj(z^)>  v^O)  .  The 
asymptotic  formulae,  used  to  approximate  the  cylinder  functions,  change 
characteristics  upon  the  passage  of  the  integration  path  through  a  turning  point 


0<  Re  p 


Table  3.  IB  -  Asymptotic  approximations  along  Integration  path 
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and,  as  a  result,  the  turning  points  separate  different  regions  of  importance 
along  the  integration  path. 

The  asymptotic  approximation  of  the  formal  solution  given  in  Eq.  (4.  2.  14) 
along  interval  Rep^  =0,  Vj(z^)  <  Impj  <  ”  Is  carried  out  on  the  top  sheet  by 
introducing  the  approximate  forms  of  the  cylinder  functions  found  in  Tables 
4.  IA  and  4.  IB,  into  the  two  integrands  and  then  integrating  these  integrands. 

Because  of  the  exponentially  increasing  nature  of  Jv  and  J_v  in  this 

L  2 

region,  both  integrals  in  Eq.  (4.  2.  14)  are  large  for  large  k^L.  Since  the 
difference  of  the  two  integrals  becomes  indeterminate  as  k  L-®,  the  repre¬ 
sentation,  given  by  Eq.{4.  2.  14), cannot  be  used  over  this  portion  of  the  integra¬ 
tion  path.  To  alleviate  this  difficulty,  we  use  the  alternative  representation 
given  in  Eq.  (4.  2.  1 6) .  Proceeding  as  before,  we  find  that  this  portion  of  the 
integration  path  gives  an  exponentially  small  contribution  tc  the  field  for  all 
observation  points  with  z>0  .  The  integration  path  was  defomed  a  small 
amount  around  the  branch  cut  in  order  to  facilitate  the  evaluation  of  the 
integral.  The  same  portion  of  the  integration  path,  on  the  second  sheet,  gives 
a  similar  result  since  the  asymptotic  approximations  used  in  the  evaluation 
remain  the  same  (they  are  independent  of  p )., 

The  original  representation  for  the  field,  i.e.  ,  Eq.  (4.  2.  14),  will  be  used 
to  evaluate  the  contribution  from  the  remaining  portion  of  the  integration  path. 
The  only  motivation  for  introducing  the  alternate  representation,  4.  2.  16,  was 
to  alleviate  the  difficulty  encountered  above  .  Before  proceeding  further  we 
shall  divide  the  field  E^  ,  as  represented  in  Eq.  (4.  2.  14) ,  in  a  sum  of  its 
individual  integrals,  i.  e.  , 

E  =  E  +  E_  (4.  3.40) 

y  Dir  Ref 

where 

EDir  =  f‘D(‘<’ V,l)dPl  '  ERef  =  J£R(Z;<  ’  S  ■  >dP|  (4’3-4‘> 
c  c 
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The  choice  of  symbols  E  and  E  was  motivated  by  the  fact  that  the 

L)~T  K6i 

evaluation  of  E-  . ^  yields  the  direct  ray  contribution,  while  the  evaluation  of 
E^^.  yields  the  returning  ray  contributions. 

The  asymptotic  approximations  for  Jv  and  J_v  ,  listed  in  Table  4. 1A, 

2  2 

are  now  used  to  approximate  fD(z<  .  z>.  p^  along  the  section  of  the  integration 
path  from  p^  ~  JL  on  the  top  sheet  to  =  JK  on  the  second  sheet .  The  inte~ 
grand  is  approximately  given  by 
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(4.  3.44) 


where 


LD(p1>  =  *2<Vp2M2(z>'p2)  ±p* 


upper  sheet 
-  lower  sheet  » 


(4.  3.45) 


For  large  kQ,  the  main  contributions  to  E^.^  will  be  in  the  vicinity  of  the 
saddle  points  of  the  integrand,  i.  e. ,  where  dL^(p,  )/dp^  =  0  .  By  performing 
the  indicated  differentiation  on  Eq.  (4.  3.  45) ,  the  saddle  point  locations  are 
given  by 

+  upper  sheet 


lower  sheet 


(4.  3.46) 
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The  above  saddle  point  equation  is  the  equation  for  the  direct  rays  when 
the  plus  sign  is  chosen.  A  sketch  of  the  different  ray  regions  has  been 
drawn  in  Fig.  4.  9.  The  basic  regions  defined  in  this  figure  are  valid  when 
L  is  not  too  small.  We  have  drawn  in  the  portion  of  regions  B  and  C, 
appearing  below  the  z-z'  line,  but  these  may  not  appear  for  some  parameter 
ranges.  When  L  is  small,  as  we  have  seen  in  the  section  on  rays,  an  add¬ 
itional  c  isp  is  introduced  into  the  region  just  above  the  z  =  0  line.  The 
sketch  of  Fig.  4.  9  cannot  be  used  to  find  the  ray  regions  in  the  immediate 
vicinity  of  this  cusp,  but  the  other  ray  regions,  with  minor  modifications, 
remain  essentially  the  same  for  small  L. 

By  using  the  nomenclature  of  Fig.  4.  9,  these  rays  are  located  in  regions 
A,B,  and  C  of  the  x-z  plane.  Since  the  rays  do  not  cross,  only  one  solution 
of  Eq.  (4.  3.  46)  exists  for  each  point  in  the  above  mentioned  regions.  The 
saddle  point  equation,  4,  3.  46,  has  no  solutions  on  the  lower  sheet  for  any 
location  in  the  x-z  plane. 

The  method  of  steepest  descent  can  now  be  used  to  evaluate  the  integral 
[Ep-rl  asymptotically.  Without  gomg  into  detail,  we  find  that  the  integral 
is  approximated  by  the  contribution  due  to  the  direct  ray  saddle  poiut,  when 
the  observation  point  is  located  in  regions  A,  B,  or  C.  On  the  other  hand,  the 
integral  is  exponentially  small  when  the  observation  point  is  located  in  regions 
other  than  A,  B,  or  C. 

The  reflected  field  integrand,  f  (z  ,  z  ,  p. )  will  now  be  asymptotically 
approximated  along  that  portion  of  the  integrand  path  from  p^  =  \JK on  the  top 
sheet  to  p^  =  i  Jh  on  the  bottom  sheet.  The  uniform  asymptotic  approximations 
found  in  section  b  can  be  used  for  this  purpose.  However,  an  additional  simp¬ 
lification  of  the  integrand  can  be  obtained  by  using  the  approximations  given 
in  Tables  1.  1A  and  4.  IB,  but  portions  of  the  integration  path  lying  close  to 
turning  points  of  cylinder  functions  must  be  omitted.  We  will  use  the  later 
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asymptotic  approximations,  and  then,  deform  the  integration  path  around 
the  turning  points. 

Before  doing  this,  it  will  be  advantageous  to  divide  the  portjon  of  the 
integration  path  under  consideration  into  six  intervals  :  three  on  the  top 
sheet  and  three  on  the  bottom  sheet.  The  intervals  arc  :  l)Repj  =  0  . 

JhJl  <  Impj<  v^z^)  ,  2)Repj  =  0.  0  <Imp^  <J A/2  and  3)  C<ReOj<<\/e^  * 
Impj  =  0  .  The  intervals  on  the  top  and  bottom  sheets  are  the  same.  If  we 
recall  that  v^(z<)  and  */V2  are  the  locations  of  turning  points.  Then  the 
approximations  to  be  used  fcr  f  (z  ,z  ,  p. )  will  not  hold  near  these  points, 

By  using  the  asymptotic  approximations  appearing  in  Tables  4.  1A  and 
4.  IB,  the  integrand  for  the  first  interval  on  the  top  sheet  is  given  approxi¬ 
mately  by 


ik  L 
o  R  > 


fR(Z<’Z>,Pl)  r  2  2  . 1/4 f  2  2.  ,1/4 

4np[p2-v2(Z>)]  [pz-v2(Z<)] 


M  <Imp1  <  v1(z>) 


(4.  3.47) 


where 


JR>=  *2+(^>.P2)  +  »2+(i<,p2)  +  px 


(4.  3.48) 


with  $2+  being  in  Table  4.  IB.  The  equation  for  saddle  points  of  this  integrand 
is  given  by  dL^>/dp^  =  0  .  By  performing  the  differentiation  on  Eq.  (4.  3.  48), 
we  find  that  the  resultant  equation  is  the  equation  for  the  returning  ray  with 
z^>0  given  in  Eq.  (4.  3.  13).  The  returning  rays  are  located  in  regions  D,  E 
and  F  of  Fig.  4.  9  . 

Proceeding  in  a  similar  manner,  the  integrand  for  the  second  interval 
on  the  top  sheet  is  given  by 

p  eik°LR<(pl)dp 

£R<WP1>  "  '  ,  2  2,1  v1l/4r  2  2, 1  .,1/4  '  0<^P1<v® 

4iTpLp?  -  v2(z<)J  lp2-v2(z>)J 


(4.  3.49) 


-no¬ 


where 

LR<  ‘  *2+^Z<’  P2^  +  *2+^Z>'  P2^  ’  2*2+^°’  p2^  +  +  px  »  (4.  3.  50) 

The  saddle  points  for  the  above  integrand  are  given  by  di^^^/dp^  =  0.  If  we 
perform  the  indicated  differentiation  upon  Eq.  (4.  3.  50),  we  find  i+  to  be  the 
same  as  the  equation  for  the  returning  rays  with  z^<0,  z  >0,  Eq.  (4.  3.  21). 
Therefore  saddle  points  appearing  in  this  region  correspond  to  returning 
rays  with  z^<0.  These  returning  rays  are  located  in  regions  B,C,D,  E, 
and  G  of  Fig.  4.  9  . 

To  show  the  correspondence  between  the  number  of  saddle  points  in  a 

particular  interval  and  their  associated  regions  in  the  x-z  plane,  we  have 

prepared  Table  4.2.  In  the  table,  only  saddle  points  of  E  have  been 

Rcl 

considered,  while  saddle  points  of  E^.^  have  been  omitted  because  of  their 
simpler  form. 

By  using  the  asymptotic  formulae  to  approximate  the  integrand  in  the 
third  interval,  we  obtain 


Region  in  x-z  pla^e 

Number  of  Saddle  Points 

0<ImpjVA/2 

^<inp1  <JK 

A 

0 

0 

B 

2 

0 

C 

1 

1 

D 

1 

2 

E 

2 

1 

F 

0 

1 

G 

1 

0 

Table  4.2 

Correspondence  between  Saddle  Points  and  Ray  Regions  of 

toteg'a!  ERef_ 


-  1  3  i  - 


fR(l<'VPl>  ' 


ik„Lr<Pl> 

r  2  2  ,l/4r  2  2  ,1/4*  0<R'P1 

4-p  [p2- v2{z<)  j  P2-V2^z>)j 


Pj  A(t,  Pj  )  e 


(4.  3.51) 


where 

Lr(pi)  =  *2+*Z<,p2*+*2+*2>,p2*_2*2+*°,P2*  ’  A(t, Pj)  -  Of-r-2)  •  (4.3.52) 


We  can  show  that  the  function  A(‘rl  p  )  in  the  above  integrand  is  of  order 

-2  1 

.  This  means  that  any  saddle  points,  occurring  in  this  interval,  will 

rot  give  a  dominant  asymptotic  contribution  to  jthe.-efore,  these  contrib¬ 

utions  can  be  neglected. 

It  is,  however,  interesting  to  learn  about  the  physical  significance  of 

these  saddle  points  if  they  occur.  The  saddle  points  are  given  by  **  ^^1^  =  0 

dpi 

An  investigation  of  this  saddle  point  equation  shows  that  it  represents 
the  equation  for  direct  rays  which  are  reflected  from  the  discontinuity  at 
z  =  0.  Only  those  direct  rays,  to  the  left  of  the  critically  reflected  ray,  shown 
in  Fig.  4.  9,  are  obtained  from  integral  4.  3.  51.  It  is  strongly  suspected  that 
direct  rays,  reflected  from  the  interface  and  lying  between  the  critically 
reflected  ray  and  the  tangent  ray,  would  be  found  from  a  higher  order  asymp¬ 
totic  analysis  of  integral  4.  3.  49  ;  however,  this  has  not  been  done. 

A  similar  asymptotic  analysis  of  the  integrand  for  the  three  intervals 
on  the  bottom  sheet  shows  that  no  saddle  points  occur  along  this  portion  of 
this  integration  path. 


Since  we  now  know  the  location  of  all  the  saddle  points  occurring  on  the 
integration  path,  the  field  E  ,  can  be  evaluated  by  using  the  method  of  steep- 

R0X 

est  descents.  We  deform  the  original  integration  path  into  the  adjacent  decay 
region  and  through  the  appropriate  saddle  points.  The  resultant  field  will 
then  be  given  asymptotically  as  a  sum  of  a  relevant  saddle  point  contributions 


^rrovronwro 


Sketch  of  Ray  Region# 
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plus  any  contributions  that  come  from  singularities  between  the  original 
and  deformed  paths.  At  this  point  we  will  restrict  our  attention  to  cases 
where  the  observation  point  lies  in  regions  B  and  C .  The  evaluation  of  the 
field  in  these  two  regions  will  illustrate  the  main  features  of  the  reflected 
field  and  will  be  useful  for  comparison  purposes  in  a  later  section. 

Two  saddle  points  occur  on  the  integration  path  when  the  observation 

point  is  located  in  either  regions  B  or  C.  For  region  B  both  these  saddle 

points  lie  in  the  second  interval,  i.  e.  ,  Rep^O,  0<lmp1  <  JEjZ  and  they 

are  shown  in  Fig.  4.  10  where  the  lower  saddle  point  is  denoted  by  and 

the  upper  one  by  p  .  The  nomenclature  has  been  motivated  by  the  fact 

that  for  small  L  the  rays  corresponding  to  and  <  look  like  the 

lateral  and  reflected  rays  on  an  abrupt  interface  respectively  .  For  region 

C  one  saddle  point,  p^  ,  lies  in  the  second  interval  while  the  second  saddle 

point,  p  ,  is  located  in  the  first  interval,  i.  e.  <v  (t  ), 

I  1  > 

Rep,  =  0.  The  naming  of  saddle  points  in  this  region  bar  a  similar  motivation 
to  those  saddle  points  associated  with  region  B. 

The  decay  regions  which  are  adjacent  to  the  integration  path  on  the  top 
sheet  have  been  investigated  by  expanding  L^Cp^).  L^<(p1)  and  L^tp^) 
about  an  arbitrary  point  in  integration  intervals  one,  two  and  three  respect¬ 
ively.  If  Pj  is  close  to  the  integration  path  the  appropriate  phase  function  can 
be  approximated  by  its  first  one  or  two  terms.  Use  of  the  approximation  in 
the  integrands  of  Eqs.  (4.  3.47)  ,  (4.  3.  49)  and  (4.  3.  Si)  yields  the  decay  regions 
for  Band  C,  shown  in  Figs.  4. 10,  4,  1 1  and  4.  12. 

The  original  integration  path  is  then  deformed  and  evaluated  along  this 
deformed  path.  The  portions  of  the  integral  passing  through  saddle  points 
are  evaluated  by  the  method  of  steepest  descents,  while  portions  of  the  inte¬ 
gration  path  in  decay  regions  are  exponentially  small. 
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The  singularities  that  can  occur  between  the  original  and  deformed 
paths  are  of  two  types  :  branch  cuts  and  simple  poles.  There  are  no  branch 
cuts  intercepted  by  the  deformed  path  as  can  be  seen  from  Figs.  4.  10,  4.  11 
and  4.  12.  The  asymptotic  form  of  the  reflection  coefficient  has  no  poles 
located  between  the  two  paths  except  possibly  in  the  turning  point  region, 

Pj  =  JKTl  .  Because  of  the  complicated  nature  of  the  reflection  coefficient 
in  this  region,  the  singularities  have  not  been  investigated  here.  We  will 
assume  that  singularities,  contributing  other  than  exponentially  small 
contributions,  are  not  present. 

Upon  applying  the  method  of  steepest  descent  when  che  observation 
point  is  located  in  region  B,  we  obtain 

(4.  3.  53) 


(4.  3.  54) 

and  doing  the  same  in  region  C,  we  obtain 


where 


ERef|B  "‘WW 
ikoLR<(pl’  +  i"/'* 


P1  "P1R< 


MP,>  * 


pl« 


1  1  z^ptp^-v^u,.)]1'4  CP2-»2<*>>i1/4r'‘0iii<<p1>liI/z 


Ref 


-  iii<P,> 


"VP1L 


wl 


(4.  3.  55) 


'P1  =  PIR> 


where 


W  4 


pi« 


ikoLR>(pl)4i,’/4 


2^  pfp^-v^)]1'4^  -  ''2<y3‘'’tkoL'R><P1)i 


1/4, 


,1/2 


(4.  3.  56) 


To  find  the  total  field  for  regions  B  and  C,  we  shall  have  to  add  the 
contribution  from  the  first  integral  in  Eq.  (4.  2.  14)  due  to  the  direct  ray. 
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The  asymptotic  formulae  obtained  in  Eqs.  {4.  3.  54)  and  (4.  3.  55.1  break 
down  as  the  observation  approaches  the  tangent  ray.  Here  a  saddle  point 
approaches  the  turning  point  at  =  JETz  and  the  simple  asymptotic  analysis 
used,  can  no  longer  be  applied. 

d.  Conclusions 


We  have  taken  the  exact  solution,  Eq.  (4.  2.  14),  and  approximated  it  for 
large  k^L,  i.  e. ,  we  have  assumed  that  the  medium  is  slowly  varying  with 
respect  to  wavelength.  The  asymptotic  approximations  developed  by  Olver 
have  enabled  us  to  shew  that  the  first  integral  in  Eq.  (3.  2. 14)  corresponds 
to  direct  ray  contributions  while  the  second  integral  corresponds  to  return¬ 
ing  ray  contributions.  In  addition,  contributions  due  to  direct  rays,  reflected 

2 

from  the  interface,  were  found  in  the  second  integral  but  were  of  order  1/k 

o 

lower  than  the  dominant  contributions.  The  above  finding  again  confirms  the 

method  of  geometric -optics  since  all  of  the  above  results  can  be  predicted 

without  resorting  to  the  asymptotic  approximation  of  the  exact  solution ;  this 

was  demonstrated  in  section  4.  3.  2.  The  reflected  waves  for  the  discontin- 

(33) 

uity  were  not  mentioned  in  section  4.  3.  2,  but  Chester  and  Keller  have  shown 

that,  for  a  dielectric  profile  having  a  discontinuous  second  derivative,  the 

2 

reflected  waves  will  be  1/k  lower  than  the  incident  wave. 

o 

4.4  Asymptotic  Evaluation  of  the  Field  for  Large  kQLp 

a.  Introduction 


In  the  last  two  sections,  4.  3.  2  and  4.  3.  3,  the  high  frequency  field  has 
been  investigated  when  the  transition  profile  is  slowly  varying  with  respect 
to  frequency.  We  shall  now  investigate  the  high  frequency  field  for  all  trans¬ 
ition  thicknesses,  k^L,  with  special  attention  placed  upon  the  changing  charac¬ 
ter  of  the  lateral  wave,  as  the  transition  thickness  increases.  The  method 
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used  will  essentially  follow  Chapter  3  ,  i.  e. ,  we  will  write  the  reflection 
coefficient  in  an  exponential  form  and  then  use  the  method  of  successive 
approximations  to  find  the  saddle  points. 


b.  Restriction  of  Field  Evaluation  to  Homogeneous  Regions 

An  examination  of  the  field  representation  given  in  Eqs.  (4.  3. 40)  ,  (4.3.41) , 
(4.  3. 42)  and  (4.  3. 43)  shows  that  both  the  integrals  for  the  direct  field,  E^^ 
and  the  reflected  field  Ej^  are  of  a  different  generic  type  than  the  integrals 
that  were  uniformly  evaluated  in  Chapter  3.  The  basic  difference  is  the 
presence  of  the  Bessel  functions  Jv^-i^e"1^)  and  J_v  (-iXe”*^*)  .  In  order 
to  apply  the  methods  of  Chapter  3,  the  Bessel  functions  must  be  approximated 
and  this  will,  in  turn,  limit  the  region  in  which  the  uniform  asymptotic  approx¬ 
imation  is  valid. 


Since  the  integral  evaluated  in  Chapter  3  resulted  from  an  exact  field 
representation  in  a  homogeneous  region,  it  would  seem  likely  that  placing 
the  observation  and  source  points  far  from  the  interface  [z  =  0)  might  put 
the  integrals  for  EDir  and  Epef  in  *■  **iore  tractable  form.  For  large  z  and 
z'  the  arguments  of  the  Bessel  functions  become  small,  and  the  Bessel 
functions  can  be  represented  approximately  by  the  first  term  of  their  series 
representation. 


sented  as 


To  put  this  on  a  more  rigorous  basis,  the  Bessel  function  can  be  repre- 
(34) 


I  (-iXe'2/L)  .(^ 
v  '  '  \  2  / 

Ct 


..  V  ik  p,  z 
-iX  ^  2  ok2 
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t -iX  \^m  -2m  z/L 
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®  /- iA.  \ 

y  oil 

A  ..  I  t 


m=0  m!  Hv^+m+l) 


(4.4.1) 


The  rate  of  convergence  of  the  series  can  be  measured  by  the  absolute 

value  of  the  ratio  of  adjacent  terms  :n  thv  series.  The  ratio  of  the  (m  +  1) 
th 

to  the  m  term  is 


% 
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.  2  -2z/L 

A  t  e _ 

8m  [(m  +  2  +  Tp..)^  +  f 

L  VL\'  K2rJ 


(4.  4.  2) 


where  p^  and  p^.  are  the  real  and  imagniary  parts  of  p^  respectively. 

An  examination  of  the  above  ratios  shows  that  a  sufficient  condition  for 
the  Bessel  function  to  be  approximately  equal  to  the  first  term  in  the  series 
representation  is  R  «  1.  A  slightly  more  restrictive  condition  than  Rq«  1 


is 


Are 


-2z/L 


«  1 


(4.  4.3) 


This  more  restrictive  condition  can  be  expressed  in  a  more  physically 
understandable  form,  i.  e. , 

(T/pz)  «'•  <4-4-4' 

As  z  becomes  large,  e(z)  tends  toward  unity,  and  the  above  condition 
requires  that  e/(z)  must  be  small  compared  with  p^/t  • 

This  more  accurately  defines  what  is  meant  by  assuming  that  the  source 
and  observation  points  lie  in  the  "homogeneous"  region.  We  note  that  rays 
close  to  the  glancing  angle  will  have  to  be  excluded  from  our  treatment  since 
these  rays  have  p^  —  0. 

It  must  be  mentioned  at  this  point  that  one  representation  for  the  field 
will  be  sufficient  for  an  asymptotic  calculation  along  the  whole  integration 
path.  If  we  recall  in  section  4.  3.  3,  two  different  representations  are 
necessary  to  calculate  the  field.  The  representation  given  by  Eq.  (4.  2. 14) 
does  not  give  an  uideterminate  asymptotic  expression  for  the  field  if  the 
source  and  observation  points  are  located  in  the  "homogeneous"  region. 

As  a  result  we  will  not  need  the  alternative  representation  in  Eq.  4.  2.  16 
in  this  section. 
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!  I 


i  ! 


By  using  the  first  term  of  the  series  representation  for  the  Besf-1 
function  in  the  integral  for  E^  ,  Eq.  (4.  3.41),  we  have 


E  * 
Dir 


p  +ik  Lp_ |z  -  z'|  +pxj 

I  O  L 


4"i  -  PP- 

C  L 


dPi 


(4.  4.  5) 


We  can  recognize  the  above  expression  as  the  incident  field  emitted  from 
a  source  located  in  a  homogeneous  medium  ha  vine;  a  dielectric  constant  of 
unity.  The  asymptotic  evaluation  of  the  above  integral  leads  to  an  interpre¬ 
tation  of  the  field  as  consisting  of  straight  line  rays.  It  can  be  shown  that 
the  "homogeneous"  approximation,  in  general,  limits  us  to  the  region  where 
the  rays  are  approximately  straight  lines. 

By  substituting  the  same  approximation  for  the  Bessel  function  in  the 
expression  for  Ej^  »  Eq.  (4.  3.  42) ,  we  obtain 
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1  Pl^  ikQ  [p2(z  +  z '  +  px] 

4TTi  pp  C  **P1 


C  PP* 
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r  = 


r(i-v2)r 

r<1  +  Vpz 


(4.  4.  6) 


(4.  4.7) 


In  the  above  expression  T  is  the  plane  wave  reflection  coefficient  as  z-*®. 
Hers  again  the  rays  in  the  "homogen ious"  region  will  be  straight  lines  but 

A 

will  have  a  phase  delay  due  to  T  . 


c.  Asymptotic  Evaluation  of  E 


Ref 


Following  Chapter  3,  we  rewrite  Eq.  (4.  4.  6)  as 


E  .  « 

Ref  4tt 
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(4.  4.  8) 
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(4.  4.  9) 

(4.  4.  10) 


The  multivaluedness  of  the  logarithm  does  not  appear  in  the  integrand 
since  the  function  e^n  ^  is  single  valued. 

The  uniform  asymptotic  evaluation  of  the  above  integral  will  be  per¬ 
formed  by  the  method  of  steepest  descents.  As  was  the  case  in  Chapter  3, 
the  saddle  points  will  be  found  by  the  method  of  successive  approximations 
and  will  have  a  location  in  the  p^  plane  that  varies  with  the  large  parameter 

k  . 
o 

The  equation  that  describes  the  location  of  the  saddle  points  is 


p  r(z  +  z')/p  Tx/p]  +  (l/k  )d*/dp  =  0  ,  “top  sheet  #  (4.4.U) 

L  °  «  +second  sheet 

We  shall  solve  this  equation  in  the  region  to  the  right  of  the  critically 
reflected  ray.  When  t  is  small  or  moderate  in  value,  the  position  of  the 
critical  ray  is  closely  approximated  by  the  critically  reflected  ray  from 
an  equivalent  abrupt  interface  problem.  When  t  is  large,  the  position  of 
this  ray  is  shown  in  Fig.  4.  9.  For  observation  points  in  the  "homogeneous" 
region,  the  critically  reflected  ray  in  Fig.  4.  9  and  the  upper  branch  of  the 
caustic  almost  coincide.  Therefore,  using  the  nomenclature  of  Fig.  4.  9 
for  large  T,  we  are  finding  the  saddle  points  in  regions  B  and  C. 

By  applying  the  method  of  successive  approximations  to  Eq.  (4.  4.  11) 
on  the  upper  sheet  when  k^  is  large,  the  first  corrected  location  for  the 
saddle  points  is  given  by 
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k  L  dp 
o  P  Ki!Pl  =  o 


,  L  -■  x  - VejM(z  t  z) 


(4.  4.  12) 


Pl(r)=  C  '  ^L1/j'/?l'oLp(<Z  +  Z,,/p2  +*/p3)(VSx+y^(i  +  z')[  |p^p<0) 


(4.  4. 13) 


where  the  superscript  is  the  order  of  approximation  with  p^  =•  ijbx.c"-(z  +  z#)t 
and  dili/dp^  =  i  1  dt/dp^  |  for  Re  p^  =  0,  0<lmp^<v^.  For  this  approximation  to 


be  valiJ  we  must  have 


di!r(p,  .  v'S’t)  / 

- — - /  L  «1  ,  0<t<»,  RePl  =0  ,  0<  ImPl  <JE .  (4.4.14) 


We  can  show  for  small  and  for  moderate  values  of  T  that  d^/dp^  is  bounded 
except  at  o^  =  0  and  for  large  values  of  t  that  dili/dp^  ~  Tf(p1)  where  f(p^) 
is  a  bounded  function  of  p  except  at  p^  =  0.  This  leads  us  to  the  conclusions 
that  the  method  of  successive  approximations  is  valid  when 

k  L  >•  1  ,  L/L  «1  .  (4.4.15) 

o  p  p 

Since  the  critically  reflected  ray  occurs  at  =  0,  these  conditions 
essentially  state  that  the  observation  point  should  be  located  a  number  of 
transition  lengths  L  from  the  critically  reflected  ray,  for  the  method  of 
successive  approximations  to  be  valid.  No  other  saddle  points  occur  in  the 
remaining  portions  of  the  integration  path  on  the  top  sheet  or  along  any 
portion  of  the  integration  path  on  the  lower  sheet. 

The  integration  path  will  now  be  deformed  in  the  adjacent  decay  regions 
and  through  the  two  saddle  points  just  mentioned  ;  it  can  be  shown  that  the 
decay  regions  for  small  or  for  moderate  t  are  similar  to  those  described  in 


4 
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Chap  ter  3  and  also  to  those  shown  in  Figs.  4.  10,  4.  1 1  and  4.  12.  For 
large  t,  the  decay  regions  have  been  calculated  in  the  last  section,  4.  3.  3.  , 
and  are  shown  in  the  three  figures  mentioned  above.  No  singularities  of  the 
branch  type  occur  between  the  original  and  deformed  paths.  The  nature  of 
pole  singularities  have  not  been  investigated  in  detail  except  in  the  case  of 
large  t  where  a  partial  investigation  shows  that  no  singularities  which  have 
a  substantial  effect  on  the  field  are  present.  We  will  assume  that  any  singu¬ 
larities  which  might  occur  between  the  two  paths  are  unimportant  and  their 
residues  can  be  neglected. 

The  integral  is  now  evaluated  along  the  deformed  path.  Portions  of  the 
path  in  the  decay  regions  have  exponentially  small  contributions  while  the 
dominant  contributions  come  from  the  portions  of  the  path  near  the  two 
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and 


M  =  rJi^T(-a)J^().) -iJo()-)J^  _{-!>-)]  TCi+i^)  .  (4.4.21) 


The  notation  E  and  E_  has  been  adopted  since,  in  the  limit  of  small  r  , 
r  L 

these  terms  correspond  respectively  to  the  reflected  and  lateral  wave 
contributions  on  an  abrupt  interface. 

The  contribution  E^  has  the  form  of  a  reflected  wave  for  an  abrupt 
interface  except  for  the  presence  of  the  phase  term,  v (p^ ,  </Kt)  .  This 
phase  term  increases  as  the  normalized  transition  thickness,  kQL  , 
increases.  For  large  T  ,  the  asymptotic  approximations  for  the  Bessel 
functions  can  be  used  in  i  (p^  ,  fK t).  We  obtain 


(plr  >^T)  ~ 


2  $i_(  Ip^  I )  -  2$z+(0,  p2r)  -  2  Lp2rln(2  ?2r  Jz7K)  -  2  Lp^ 

0<  |p,  |  <  s/572 


(4.  4.  22a) 


2  Lp2rln(2p2r‘v/^‘^_  2  Lp2r  *  <  lplrl<v}(z>) 


(4.  4.  22b) 


where  and  $2+  have  been  defined  in  Table  4.  IB  and  p2  =  Jk  -  | |  ^  . 

The  use  of  these  approximations  in  Eq.  (4.4.  17)  gives  us  a  portion  of  the 
reflected  wave  when  the  transition  is  slowly  varying  with  respect  to  frequency. 
This  portion  of  the  reflected  field  corresponds  to  contributions  from  rayo 
that  make  up  the  section  of  the  caustic  close  to  the  interface  as  shown  in 
Fig.  4.  5. 

The  contribution  E  has  the  form 
L 

E^  ~  A(  s/A  t)  e'^^  ^  |  conventional  lateral  wave  contribution^ 

(4.  4.23) 
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Here  A(>/&"t)  and  iU/K t)  form  a  transition  function  which  modifies  the  form 
of  the  lateral  wave.  For  small  t,E.  is  approximately  represented  by  the 

li 

conventional  lateral  wave  contribution  since  A(v^t)**1  and  t(^T)s»0.  As 
increases,  A{JK r)  and  t)  increase  in  such  a  way  as  to  modify  the  lateral 
wave  teem  and  change  it  into  a  geometric -optics  contribution.  For  large  r, 
the  asymptotic  forms  of  the  Bessel  functions  can  be  used  to  approximate 
A{JS t)  and  .  We  obtain 


...  ik  \jE(z  +  z')  +  JF.  x  +  ZjKtZ/JT-  1  -ln( 
>eJ/4Le  °l  1 

zJzSl  (k  L  f 

pop 


,  T»  1  . 


(4.  4.  24) 


-3/2 

We  note  that  the  dependence  for  small  t  in  Eq.  (4.  4.  18)  has  now 

been  replaced  by  a  k  ^  dependence.  This  signifies  the  change  in  the 
contribution  Ej^  from  a  diffraction  effect  to  a  geometric -optic  effect.  The 
amplitude  A(a/S’t)  and  the  phase  t)  have  been  plotted  in  Figs.  4.  13 
and  4.  14  respectively.  The  dash  line  in  these  figures  represents  the  slowly 
varying  approximation,  i.  e.  ,  when  A(V^t)  and  are  asymptotically 

approximated  for  large  t.. 


The  approximate  forms  of  and  E^  for  large  r  can  also  be  obtained 
from  the  results  of  section  4.  3.  There  the  field  in  region  B  is  given  by 
Eq.  (4.  3.  53)  while  the  field  for  region  C  is  given  in  Eq.  (4.  3.  55).  If  we 
simplify  the  expressions  by  placing  the  source  and  observation  points  in 
the  "homogeneous"  region  and  assume  L/L^»l  ,  these  reduce  to  the  large 


T  results  of  this  section.  To  be  more  explicit,  the  term  il^(p^)|p^  _  p 


reduces  to  Eq.  (4.  4.  24)  ;  the  term  I.  (p  ) 


reduces  to  (4.4.17) 


1L 


W 


lpl  =  P1R> 


rrripi  =  piR< 

i)  an 

reduce  to  Eq.  (4.  4.  17)  where  in  Eq.  (4.  4.  17)  is  approxi- 


with  ♦(p^  *  t)  approximated  by  Eq.  (4.4.  22a)  and  finally  the  term 


mated  by  Eq.  (4.  4.  22b). 


Fig.  4.  13 

Lateral  Wave  Amplitude  vs.  Normalized  Transition  Thickness 


Lateral  Wave  Phase  vs.  Normalized  Transition  Thickness 
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d.  Conclusions 


Our  basic  purposes  for  treating  the  double  exponential  medium  were 
twofold:  first,  to  see  if  the  lateral  wave  behavior  observed  on  the  Epstein 
half  space  carried  over  to  profiles  where  the  source  was  located  in  the 
inhomogeneous  medium,  and  second,  to  see  the  effect  of  relocating  the 
second  order  discontinuity  in  the  dielectric  profile.  We  found  that  by 
restricting  the  source  and  observation  point  locations  to  the  "homogeneous" 
region,  results,  similar  to  those  in  Chapter  3  were  found.  The  lateral  wave 
changed  to  a  geometric -optics  wave  as  the  transition  went  from  abrupt  to 
slowly  varying.  It  was  also  observed  here  and  in  Chapter  3  that  the  depend¬ 
ence  on  the  lateral  distance  L.^  remained  .  independently  of  the  trans¬ 
ition  length.  The  discontinuity  in  the  dielectric  profile  did  not  affect  the 
basic  nature  of  the  lateral  wave  expression  but  it  did  introduce  a  difficulty 
in  the  reflected  field,  E^  .  For  large  transition  length,  the  field  close  to 
the  tangent  ray  could  not  be  treated  by  the  classical  methods  of  geometric- 
optics. 
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CHAPTER  5. 

AN  ARBITRARY  TRANSITION  LAYER 


5. 1  Introduction 


To  conclude  our  investigation  of  lateral  waves  which  are  excited  on 
transition  layers,  we  shall  discuss  and,  in  some  cases,  analytically 
formulate  the  dominant  asymptotic  contributions  to  the  reflected  field  from 
an  arbitrary  monotonically  increasing  layer.  The  layers  to  be  considered 
are  divided  into  two  basic  classes  :  finite  and  infinite.  The  finite  layers  cif 
width  L  will  consist  of  a  dielectric  medium  varying  continuously  between  two 
half  spaces,  one  being  composed  of  vacuum  (z>0)  and  the  other  (z<-  L)  having 
a  dielectric  constant  of  «j<  1  .  The  dielectric  variation  of  the  layer  should  be 
strictly  increasing,  i.  e:  de(z)/az>0  except  at  the  lower  interface  where 
(n-1)  derivatives  of  e(z)  may  be  zero.  To  be  more  explicit,  the  dielectric 
profile  c(z)  has  a  series  expansion  at  the  lower  interface  which  is  given  by 


e(z)  = 


..♦■a-dsi 

1  j  n 

dz 


(z+L)n 
n  ! 


{5.1.1) 


The  infinite  layer,  on  the  other  hand,  will  again  consist  of  a  vacuum 
half  space  for  z>0  and  a  strictly  increasing  dielectric  profile  (z<0)  which 
approaches  asymptotically  as  z-»-®  .  In  addition,  we  will  require  that 
all  derivatives  of  e(z)  approach  zero  z-*-  ®  and  that  e(z)=  e(z/L)  where  L 
is  a  parameter  which  is  proportioned  to  the  average  thickness  of  the  layer. 
Both  classes  of  layers  have  this  property:  as  L  becomes  small,  the  medium 
approaches  an  abrupt  transition. 

The  line  source,  as  in  previous  chapters,  will  be  located  parallel  to  the 
z  =  0  interface  in  the  vacuum  medium  .  Sources  placed  in  the  inhomogeneous 
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portion  of  the  layer  will  not  be  considered,  since  the  results  of  Chapter  4 
indicate  that  such  a  source  location  has  little  effect  on  the  basic  character 
of  the  lateral  wave  excited. 

The  structure  of  this  chapter  will  be  similar  to  previous  chapters.  First, 
the  reflected  field  from  a  layer  that  is  thick  compare!  with  wavelength  will  be 
discussed  and  then,  an  asymptotic  evaluation  of  the  reflected  field  representa¬ 
tion  will  be  carried  out  when  k  L  »1  . 

o  p 

5.  2  Reflected  Field  from  an  Arbitrary  Layer  for  k>L»l 

The  reflected  field  from  a  layer  that  is  thici:  compared  with  wave¬ 
length  can  be  obtained  asymptotically  by  the  classical  methods  of  geometrical- 
optics  with  the  exception  of  those  regions  where  diffraction  effects  are  import¬ 
ant.  In  this  section  we  would  like  to  locate  those  diffraction  regions  for  an 
arbitrary  layer  and  try  to  predict,  by  using  the  results  of  previous  chapters, 
what  type  of  diffraction  effects  we  can  expect. 

Before  doing  this,  however,  the  ray  trajectories,  as  predicted  by  classical 
geometrical- optics,  will  be  investigated.  The  rays  emitted  from  the  source 
toward  the  layer  divide  into  two  ray  types  :  transmitted  and  returning ;  the  two 
ray  types  are  separated  by  a  critical  ray.  The  returning  ray  equation,  as 

,  ~  i  (20) 

given  by  Orlov  ,  is 


X  =  x  +  x.z/z'  ,  z>0  (5.2.1) 

o  i  — 

where  x^ ,  the  coordinate  at  which  the  ray  enters  the  medium,  can  be  related 
to  the  ray  parameter  p  used  in  previous  chapters.  Their  relationship  is 
x^p^z/p  •  The  coordinate  xq  is  the  point  at  which  the  returning  ray  leaves 
the  layer,  as  shown  in  Fig.  5.  1,  and  is  related  to  x^  by 

x  =  x.  +  2  x, 
o  1 


(5.  2.2) 


Typical  Returning  Ray  Trajectories 
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Here  x^  is  the  x  coordinate  of  the  turning  point  which  is  given  by 


Ct=-  J 


pdT 


Jz{ T)-  p2 


(5.2.3) 


in  terms  of  the  ray  parameter  p  . 

The  rays  described  by  Eq.  (5.  2.  1)  usually  cross  one  another  in  such  a 
way  as  to  form  a  caustic  and  foci.  The  constraint  equation  of  this  caustic 
can  be  found  by  taking  the  derivative  of  Eq.  (5.  2.  1)  with  respect  to  x.  .  The 
result  is 

dx  /dx.  =  z/z'  .  (5.2.4) 

o  1 


We  see,  frcm  this  equation,  that  the  caustic  goes  into  the  layer  at  points 
given  by  dx^/dx.^  0  .  An  additional  point  or  points  on  the  caustic  can  be 
found  by  using  the  result  of  Appendix  D  where  it  has  been  shown  that  the 
maximum  of  the  locus  of  turning  points  (dx^/dz^=0)  corresponds  to  points 
at  which  the  caustic  crosses  the  line  z  =  z  , 

The  locations  of  possible  foci  are  given  by  the  second  constraint  equation 

d^x  /dx^  =  0  (5.2.5) 

o  1 

when  used  in  conjunction  with  Eqs.  (5.  2.  4)  and  (5.  2.  1)  .  The  constraint 
equation  was  found  by  taking  the  derivative  of  Eq,  (5.  2.  4)  with  respect  to  x.  . 

Now  that  we  have  outlined  the  means  by  which  we  can  obtain  the  structure 
of  the  reflected  ray  family,  as  predicted  by  classical  geometrical- optics,  we 
will  discuss  diffraction  effects.  To  simplify  the  following  discussion,  it  will 
be  ?3sumed  that  the  profile  of  the  arbitrary  layer  is  completely  continuous  for 
-L<z<0  ,  i.  e:  all  derivatives  of  e(z)  are  continuous  in  the  layer  region. 
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This  assumption  eliminates  diffraction  regions  occurring  in  the  vicinity 
of  rays  which  have  their  turning  points  located  at  discontinuities  of 
dne(s)/dzn  . 

Diffraction  regions  are  those  regions  in  which  the  classical  theory  of 
geometrical-optics  gives  a  poor  approximation  to  the  field.  Typical  regions 
of  this  type  are  caustics,  foci  and  ray  regions  composed  of  rays  emitted  close 
to  the  critical  angle.  It  is  the  latter  of  these  diffraction  effects  that  we  wish  to 
discuss.  The  class  of  layers,  under  discussion,  will  be  divided  into  three 
parts  :  first,  layers  with  a  discontinuity  in  slope,  n=  1  ;  second,  layers  with 
discontinuities  in  higher  derivatives,  n>2  ;  and  third,  infinite  layers.  Each 
sub-class  of  layers  will  be  discussed  individually. 

The  layers  with  n  =  1  have  been  singled  out  since  these  layers  have 

critical  rays  whose  turning  points  have  finite  coordinates.  The  results  of 

Chapter  1  indicate  that  a  lateral  ray  should  be  exicted  at  this  turning  point 

and  it  should  travel  along  the  z  =  -  L  interface  while  shedding  energy  into  the 

1/3  3/2  it 

reflected  field.  The  general  form  of  this  wave  would  be  A(  (kL)  / (k  L  )  )e  *. 

o  os 

Here  f  is  the  phase  from  the  source  to  the  observation  point  along  the  lateral 

ray ;  L  is  the  distance  that  the  ray  travels  along  the  z  =  -  L  interface  and  A 
8 

is  the  excitation  coefficient  which  will  depend  on  the  slope  of  the  layer  at 
z  =  ~  L  .  We  note  that  the  k  dependence  is  smaller  than  the  returning  ray 
contribution,  however,  the  lateral  wave  contribution  might  be  important  if 
there  were  any  loss  in  the  layer.  The  results  of  Chapter  1  also  seem  to 
indicate  that  the  returning  rays  which  lie  close  to  the  reflected  critical  ray 
have  field  contributions  that  are  not  predictable  by  the  theory  of  geometrical- 
optics.,  The  field  in  this  region  could  be  obtained  by  an  asymptotic  analysis  of 
the  integral  representation  for  the  reflected  field. 

The  next  class  of  layers  to  be  considered  are  those  layers  with  n>2  . 

An  examination  of  the  critical  ray  in  this  class  of  layers  shows  that  it  is 
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similar  to  the  critical  ray  observed  in  Chapter  2.  There  the  critical  ray 

never  becomes  tangent  to  the  lower  interface  but  only  approaches  it 

asymptotically  i.  e. ,  as  z  L,  x  -«  +  ®  where  x  and  z  are  the  coordin- 

c  c  c  c 

ates  of  the  critical  ray.  Because  this  ray  never  becomes  tangent  to  the 
lower  interface,  the  geometrical  interpretation,  rendered  to  the  lateral 
wave  in  the  n  =  1  case,  does  not  seem  applicable  here.  However,  it  does 
appear  that  returning  rays,  emitted  close  to  the  critical  angle,  form  a 
diffraction  region. 

To  see  the  behavior  of  these  returning  rays  for  a  special  case,  we 
have  chosen 

e(z)  =  el  +  A  (z/L  +  1  )n  ,  n >3  .  (5.2.6) 

We  have  not  included  n  =  2  since  an  example  of  this  case  has  been  given 
in  Chapter  2.  The  ray  equation  (5.  2.  1) ,  when  rewritten  in  terms  of  the 
ray  parameter  p,  becomes 


x  =  p(z  +  z')/p2  +  2xt 


(5.2.7) 
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The  dominant  asymptotic  approximation  for  x^ 
given  by 


x 

t 


WeJn 


1M 


n-4 

n 


as 


|pj  |  -0  is  then 
n>3 


(5.2.  10) 


where  1^  is  a  convergent  integral  given  by 


CO 


We  see  as  |p^  |  -*0  ,  the  slope  of  the  ray  tends  toward  ./A/e ^  while  the 
intercept  at  the  x  axis  tends  toward  infinity.  As  in  the  n=  2  case,  these 
rays  have  trajectories  similar  to  the  lateral  wave  observed  in  the  linear 
layer  in  Chapter  1, 

We  see  that  rays  in  this  type  of  medium,  which  are  emitted  at  angles 
close  to  the  critical  angle,  contribute  to  the  reflected  field  at  large  distances  . 
from  the  source.  The  results  of  Chapter  2  would  also  indicate  that  the 
methods  of  geometrical- optics  are  not  generally  successful  in  computing 
the  amplitudes  of  these  ray  contributions. 

The  obvious  step  to  take  after  considering  profiles  with  n>  2,  is  to 

investigate  an  analytic  layer .  At  present  this  has  not  been  done  for  a  finite 

layer,  since  no  analytic  profile  could  be  found  which  has  known  wave 

functions  associated  with  it.  Recourse  was  taken  to  infinite  layers,  such  as 

the  symmetrical  Epstein  transition,  for  which  the  wave  functions  were  known. 

An  examination  of  the  critical  ray  in  an  infinite  layer  shows  that  as 

x  -»®  ,  z  ®.  This  differs  from  the  finite  layers  where  z  -♦  -  L  .  The  results 
c  c  c 

of  Chapters  3  and  4  have  indicated  that  returning  rays  emitted  close  to  the 
critical  angle  tend  to  act  in  a  similar  manner  as  the  recurning  rays  from 
finite  layers  with  n>2  ,  with  one  exception.  This  exception  is  that  as 
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Ipj  I  —  0  ,  the  turning  point  z^  tends  toward  infinity.  The  results  of 
Chapters  3  and  4  also  indicate  that  the  method  of  geometrical-optics 
does  not  break  down  as  the  angle  of  ray  emission  from  the  source 
approaches  the  critical  angle.  Therefore  no  diffraction  region  is  observed. 

5.  3  Evaluation  of  the  Reflected  Field  for  k  L  »1 

_ °  P 


5.  3.  1  General  Considerations 

In  this  section  we  will  investigate  the  reflected  field  from  an  arbitrary 

transition  when  k  L  »1  and  L  »L.  From  this  investigation  we  will 
o  p  p  s 

determine  the  nature  of  the  diffraction  field  at  large  distances  from  the 

source  when  k  L»1  ,  and  relate  this  diffraction  field  to  the  lateral  wave 
o 

contribution  v  hich  is  excited  on  an  abrupt  interface.  The  analytical  treat¬ 
ment  will  be  limited  to  finite  layers  ;  infinite  layers  will  be  discussed  by 
using  the  results  of  Chapters  3  and  4. 

5.  3.  2  Finite  Layers 

The  integral  representation  for  the  reflected  field  from  an  arbitrary 
finite  layer  is  given  by 

ik  [p  (z+z')  +  px] 

e  °  dpj  ,  z  > 0  (5.3.1) 


y* 


1 

4tt  i 


\ 

C 


ill 

pp. 


where 


r  *  -  »,  /4b 


(5.  3.2) 
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with 


J*;(-L)+ikP,  u-dI 

j  2  u  1  <.  j 

J»^(o)±ik)p2«2(o)j 


ikoP2 


The  integral  representation  is  the  same  as  the  integral  representation 
appearing  in  Eqt  (1.4.  1)  with  one  exception  ;  their  reflection  coefficients 
are  different.  The  reflection  coefficient,  appearing  in  Eq.  (5.  3.  1)  ,  has 
been  derived  in  Appendix  A.  In  this  Appendix,  we  have  defined  *^(z)  and 
♦^(z)  as  two  independent  solutions  to  the  reduced  wave  equation  in  the 
layer  region.  The  square  roots  p  =  -  p^  ,  p^  =  Jb  +  p^  and  the  path  C 

are  defined  in  Figs.  1.5  and  1.6. 


If  we  now  use  the  assumption  that  k  L«k  L  ,  then  the  reflection 

o  o  p 

coefficient  in  Eq.  (5.  3.  1)  is  slowly  varying  compared  to  the  exponential 
in  the  integrand,  and  the  standard  steepest  descents  techniques  can  be  used 
to  asymptotically  evaulate  the  integral  representation  for  large  .  The 

asymptotic  evaluation  parallels  the  evaluation  of  the  integral  in  Eq.  (1.4.  1) 
and,  as  a  result,  we  will  not  present  it  here.  The  results  of  this  asymptotic 
analysis  are  that  two  saddle  point  contributions  make  up  the  reflected  field 
when  k  L  »1  .  The  first  of  these  is  a  contribution  from  a  ray  reflected 


o  p 


-1/2. 


from  the  z  =  0  interface.  The  contribution  has  a  0(kQ  )  for  all  layer 
thicknesses.  When  kQL»l,  the  contribution  can  be  predicted  by  the  class¬ 
ical  theory  of  geometrical-optics. 


The  second  contribution  to  the  reflected  field  comes  from  the  saddle 
point  occurring  at  p^  =0,  as  shown  in  Chapter  1.  This  contribution,  which 
will  be  denoted  by  E.  is  given  by 
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where 


1/4  dr 

1  dPL  Pl  =  0 


2v/2tTa  (k  L  )' 
o  p 


ik^C^fz+z')  +  ,/ejx]  -  i tt/4 


(5.3.3) 


Ipa  =  0 


2  T  JE  [§2(o)  $'(o)  -  ^(oH^o)]2 

|[#2'  L)f/1(o)-$/1(-L)$^(o)]-ikoVS'[$/2(-L)$1(o)-$/1(-L)$2(o)]}2 


(5.3.4) 


t  =  k^L  and  the  prime  indicates  differentiation  with  respect  to  z/L..  The 
above  contribution  can  be  written  in  the  following  form 


the  lateral  wave 
contribution  on 
an  abrupt  interface. 


(5.3.5) 


where 


JK  dr 

2  dp, 


Ipj  =0 


(5.3.6) 


In  the  above  T(t)  is  a  transition  function  which  is  independent  of  the 

observation  point  if  we  keep  within  our  original  assumption  that  k  L»  1 

and  L  » L  . 

P 

It  would  be  interesting  to  learn  what  particular  characteristics  of  the 
layer's  profile  affect  the  transition  function  T(t)  ,  however  the  wave 
functions  $j(z)  and  $  (z)  are  not  known.  As  a  result,  we  will  limit  our 
investigation  of  T(t)  to  values  of  t  which  are  small  or  large  compared 
with  unity. 
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We  will  proceed  first  with  an  investigation  of  T(")  for  small  t  . 

An  approximate  expression  for  T(t)  when  t  is  small  can  be  obtained  by 
expressing  F  in  a  power  series  in  t  .  We  can  obtain  this  power  series  by 
considering  the  generalized  reflection  coefficient  F  (z)  in  the  layer.  This 
reflection  coefficient  obeys  the  Ricatti  equation  as  given  by  Brekhovskikh^  1 
i.  e: 

dF/dz  =  2ik  p  (z)T(z)  +  e(z)  [l-r  2(z)]/p2(z)  (5.3.7) 

o  z  z 

where  p  (z)  =  Jc(z)-pZ  .  This  reflection  coefficient  is  the  ratio  of  incident 
z 

to  reflection  wave  fields  at  any  point  z.  As  a  result  when  z  becomes  less 
then  -L,  the  reflected  wave  disappears  and  r(-L)  =  0.  We  use  this  condition 
as  the  one  initial  condition  necessary  to  make  the  solution  to  Eq.  (5.  3.  7) 
unique. 

Actually  the  function  F(z)  defined  by  this  equation  has,  in  general,  no 

physical  meaning  since  incident  and  reflected  waves  are  indistinguishable 

in  an  inhomogeneous  medium.  However,  when  the  medium  is  slowly 

varying,  compared  to  wavelength  or  z  is  located  in  an  homogeneous  region, 

this  difficulty  no  longer  occurs  and  the  reflection  coefficient  takes  on  the 

meaning  that  we  usually  ascribe  to  it.  In  our  particular  case,  F(z)  will  be 

evaluated  at  z  =  0+  ,  a  poht  where  T(z)  is  physically  meaningful.  Fora 

more  complete  discussion  of  Eq.  (5.  3.  7)  and  its  physical  interpretation  we 

(35) 

refer  the  reader  to  Schelkunoff 

Brekhovskikh^  ^  ^  and  Wait^^  have  developed  a  procedure  for  expres 
ing  F(z)  in  terms  of  an  ascending  power  series  in  k^  which  would  be  useful  for 
our  purposes.  Unfortunately  the  above  series  does  not  converge  for  p^  =  0 
and  it  is  at  this  point  that  we  require  the  series  development.  To  alleviate 
this  difficulty  we  have  slightly  modified  Brekhovskikh's  method  and  developed 
a  series  which  does  converge  at  p^  =  0  . 
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At  this  point  we  will  assume  e(z)  =  e(z/L).  This  assumption  is  not 
necessary  but  it  will  simplify  the  analysis  to  follow  and,  at  the  same  time, 
allows  for  a  sufficient  amount  of  generality.  The  transformations 


F(  z  ) 


_  t  \  \  /  —  \ 

FzV^;  vi*  ;  -  u\  *  ; 

p  (7)  v(7)  +  u(z  ) 

Z 


z  =  z/  L 


(5.  3.  8) 


where  p  ( z  )  will  now  be  used  in  Eq.  (5.  3.  7)  and  an  equation  for  v(  z  )  and 
z 

u(7)  will  be  obtained.  It  is 


uVu  -  vVv  = 


u _ 

Pz(z)  V 


p  (z  )v 
z 

u 


] 


(5.3.9) 


where  the  prime  indicates  differentiation  with  respect  to  ?  .  This  equation 
can  be  satisfied  by  subjecting  u  and  v  to  the  two  equations 

u/;=-iTp^(7)v  (5. 3. 10) 

z 

v'  =  -iTU  .  (5.  3.  11) 

❖ 

The  solutions  of  the  above  equations  with  the  boundary  conditions 

u(-l)  =  PL  ,  v(-l)  =  1 

yields  Eq.  (5.  3.  7)  with  the  boundary  condition  that  r|^  _  ^  =  0  .  Now 

combining  the  Eqs.  (5.  3.  10)  and  (5.  3.  11)  we  obtain  a  second  order  equation 
for  v(z) .  It  is 

v,/  =  -t^p^(T)v  (5. 3. 12) 

z 


The  boundary  conditions  are  imposed  on  the  negative  side  of  the 
discontinuity  at  z  =  -  L  . 
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with 

v(- 1 )  =  1  ,  v/(-l)=-iTp1  .  (5.3.13) 


We  will  now  convert  Eq.  (5.  3.  12)  into  an  integral  equation  by  considering 

the  Green's  function  problem  v"  =  -  6(z-z)  and  by  considering  the  right  hand 

(37) 

side  of  Eq.  (5.  3.  12)  as  a  source  term.  By  using  Friedman's  result  for 
the  Green's  function,  Eq.  (5.  3.  12)  becomes 


v(  z  )  =  1  -  i  T  p.  z  -  f  p2(t)  (  z  -  t)  v(t)  dt  .  (5.  3.  1 4) 

1  -J1  z 

(38) 

With  the  aid  of  Tricomi  ,  it  can  be  shown  that  the  method  of  successive 

approximations  with  1  -irp^z-  ,  as  the  zeroth  approximation,  converges 

to  the  solution  of  Eq.  (5.  3. 14).  The  only  requirements  are  that  1-  i  Tp  z 
2  2  — 

and  T  p  (t)(z-t)  be  square  integrable  in  the  interval  -1<Z<0  .  We  take 
z 

special  note  that  this  solution  will  converge  at  p^  =  0  . 


By  applying  the  method  of  successive  approximations  to  find  v("z  )  and 
then  by  using  Eq.  (5.  3.  11)  to  find  uCz  ),  we  can  obtain  an  approximation  to 
T(0)  for  small  t  .  It  is 


r  =  r(o)  = 


p2-pi 
p2  +  Pl 


2  i  t  p. 


r  0 

Hi 


(e(t)-  e^dt-p2] 


+  0(t  )  ,  (5.  3.  15) 


<p2  +  pi> 


If  we  now  take  the  derivative  of  T  with  respect  to  p^  ,  we  can  obtain  an 
approximate  expression  T(t)  when  t«1  .  This  expression  is 

T(t)  =  1  +  |e(7)  -  6lJ  +  0(r2)  (5.3.16) 

where  the  term  in  the  brackets  is  the  average  of  e  (z  )  -  ,  i.  e.  , 
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1.°  o 

G(z)-CJ  =  f  (C(t)-^)dt=  -  f  (e^D-ep dz  .  (5.3.17) 

'av  -1  -L 


If  we  refer  to  Eq.  (5.  3.  5)  we  see  that  as  t  -*0  ,  the  lateral  wave  contribu¬ 
tion  E  tends  toward  the  lateral  wave  excited  on  an  abrupt  interface.  The 
first  perturbing  term,  as  can  be  seen  from  Eq.  (5.  3.  16),  is  proportional 
to  the  average  amount  that  the  dielectric  constant  exceeds  . 

Since  we  have  obtained  the  approximate  form  of  T(t)  when  t«1  ,  we 

will  now  investigate  T(t)  when  t»1  .  The  general  procedure  for  doing 

this  will  be  to  asymptotically  approximate  the  function  $^(z)  and  ^(z) 

appearing  in  Eq.  (5.  3.  4).  The  asymptotic  approximations  to  these  functions 

(39) 

will  be  obtained  by  using  Langer’s  comparison  method.  This  method 
allows  us  to  find  the  asymptotic  solution  to  one  equation  in  terms  of  the 
solution  of  another.  For  the  procedure  to  apply,  the  two  equations  must 
have  turning  points  of  equal  order. 

The  unknown  functions  $j(z)  and  $ ^ (z)  obey  a  reduced  wave  equation 


2 

\^2  f  j2  [  ej  I  $i(z)  =  0  *  z  -  z/L 

'  Hz  ' 


(5.  3.  18) 


with  the  boundary  condition  that  $^(z)  and  must  he  outgoing  and 

incoming  waves  respectively,  when  k  L  is  large.  We  have  assumed,  as 

_  o 

before,  that  e(z)  =  e(  z  ) .  Since  Eq.  (5.  3.  18  )  has  one  turning  point  of  n  order, 
we  choose  as  a  comparison  equation 


■V  tV[v5>=o 

us  2 


(5.3.  19) 


where  ^(z)  must  obey  the  same  boundary  conditions  as  $ j(z).  The  above 

2  2 
equation  can  be  transformed  into  a  Bessel  equation  and  has  the  following 


-162- 


solutions 


4.(5)  =  JC 

1  ’  ‘ 

2 


1  ln+2 
n+2 


/c  7  otw 

J  •  4#  V  f 


By  applying  Langer's  method  with  Eq.  (5.  3.  19)  as  the  comparison  equation, 

the  asymptotic  approximations  i^(z)  are  given  by 

2 


J  (e(T)-  e^dT 

[e(J)-c,  ]'/4 


<2>  f  i  1 

H  1  T  |  (e(T)  -  G1)SdT  .  (5.  3.21) 

n+2  ^ 


We  shall  also  require  the  asymptotic  approximation  for  d$^(z)/dz  .  This 

2 

can  be  obtained  by  taking  the  formal  derivations  of  Eq.  (5.  2.  21)  and  retain¬ 
ing  only  the  higher  order  term.  The  result  is 


dS^z) 


-  T[«<7>-«if4[aT  fWTi-sAP  H-nJTf1(6(,)-£l,idT] 


(5.  3.  22) 


An  examination  of  Eq.  (5.  3.  4)  shows  that  we  will  require  ant* 

their  derivative  at  z  =  0  .  When  t  is  large,  the  argument  of  the  Hankel 
functions  in  Eqs.  (5.  3.  21)  and  (5.  3.  22)  becomes  large.  By  using  the 
asymptotic  expansion  for  Hankel  functions  with  large  argument  and  fixed 


orders,  we  obtain 


*i  * 

2 


±i  ^jWl-e^dT  -  4^2)n 


(5.  3.23) 
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*i  (o) 


/(n+2  )A 


frf  (c(r> 

1  -1 


i 

V  dr 


n+4 

4(n+2) 


(5.  3.24) 


The  evaluation  of  Eq.  (5.  3.  4)  will  also  require  knowledge  of  5^(-L)  .  This 

2 

can  be  obtained  by  using 


HZv(z)  *  +  -J-lSlil -  p(v)  e±TrV  ,  z«  1 


(5.3.25) 


in  Eq.  (5.  3.  22) .  The  result  is 


*;(-l)  -  + 


f  i  "ll/n+2  ± .  _ 

,  f  c(n)  ]  r  (  n+l\ 

T(n+2)[-^-j  J  r,  — )  e 


n+1 

ITT  - 7 

n+2 


(5.3.26) 


where  e(n)  =  d  |— __j.  If  we  now  use  Eqs.  (5.  3.  23) ,  (5.3.24) 

and  (5.  3.  26)  in  Eq.  (5.  3.  4)  we  find  that  T(t)  is  given  by 

r  _  inTT 


T(t)  ~At 


v  a 


(£.3.27) 


when  t»1  .  The  constant  A  is  given  by 


TT.yS” 


‘-a! - ln+2 

(n+2)n  e(n)J 


.2  /  n+l\ 
\  n+2/ 


(5.  3.  28) 


An  inspection  of  Eq.  (5.  3.  2  7)  shows  that  as  n  increases,  the  order  of 

T(t)  rises.  This  means  that  the  more  continuous  the  dielectric  profile  is 

at  z  =  -  L,  the  stronger  the  lateral  wave  is  excited.  In  fact,  as  n  tends 

1/2 

toward  infinity,  the  lateral  wave  contribution  approaches  an  0(kQ  )  .  This 
means  that  the  lateral  wave  contribution  will  be  as  strong  as  the  geometrical- 
optic  contribution. 
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It  is  interesting  to  compare  the  lateral  wave  field  with  the 

reflected  field  discussed  in  section  5.  2  .  There,  when  n  =  1 ,  a  lateral 

7/6 

wave  is  excited  and  has  an  0(k  '  )  .  As  the  lateral  distance,  L>  , 

o  p 

becomes  great  compared  to  wavelength,  the  lateral  wave  contribution 
found  in  section  5.  2  becomes  identical  with  the  lateral  wave  contribution 
given  by  Eq.  (5,  2.  3)  with  Eq.  (5.  3.  27)  and  n  =  1  .  We  also  note  in  section 
5.  2  that  when  n>2  ,  the  reflected  field  is  composed  of  two  contributions. 
One  of  these  contributions  can  not  be  obtained  by  the  methods  of  geomet¬ 
rical-optics  when  L  becomes  large  compared  with  wavelength.  In  this 
section  we  see  that  the  geometrical-optic  contribution,  observed  in 
section  5.2,  makes  a  transition  to  a  lateral  wave  type  contribution  as 

k  L  »1.  We  have  made  no  attempt  to  find  the  transition  function  between 
op 

these  two  contributions. 

5.  3.  3.  Infinite  Layers 

We  will  not  treat  the  class  of  infinite  layers  analytically,  as  has  been 

done  for  finite  layers  in  the  last  section,  but  instead,  we  will  discuss  the 

contributions  that  we  might  expe  t  on  the  basis  of  the  results  ox  Chapters 

3  and  4.  The  results  of  both  of  these  chapters  indicate  that,  as  in  the 

finite  case,  the  reflected  field  for  k  L  »1  consists  of  two  contributions. 

o  p 

One  of  these  contributions  appears  to  come  from  a  r ay  reflected  fvom  the 

z  =  0  interface  while  the  other  contribution  has  a  ray  trajectorv  similar  to 

“  3/2 

a  lateral  ray.  The  latter  contribution  has  an  algebric  decay  factor  of 
which  is  independent  of  k^L  .  This  behavior  is  also  observed  in  sections 
5.  3.  2.  The  excitation  coefficient  of  the  contribution  changes  with  t  .  For 
small  t  the  contribution  appears  to  reduce  to  a  lateral  wave  on  an  abrupt 
interface  while,  for  large  T  ,  the  contribution  reduces  to  a  geometrical- 
optic  ray.  It  is  interesting  to  note  that  the  lateral  wave  contribution, 
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investigated  in  the  last  section  for  arbitrary  n,  approaches  a  0(k^ 

as  n-.°°.  This  is  the  same  order  as  a  geometrical-optic  contribution.  The 

interesting  contrast  between  the  lateral  type  contributions  on  finite  and 

infinite  layers  for  k  L  »k  L»1  is  that  in  the  infinite  case  the  methods 

o  p  o 

of  geometrical- optics  can  be  used  to  obtain  the  lateral  contribution  while 
in  the  finite  cases  it  cannot. 


» 
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APPENDIX  A. 

ONE-DIMENSIONAL  GREEN'S  FUNCTION 
PROBLEM  FOR  AN  ARBITRARY  LAYER 


The  function  #(z,  p)  is  a  one-dimensional  Green's  function.  It  obeys 
the  equation 


j-~  +  Je(z)-p2]|  $(z,p)  =  -  6{z  -  z) 


(A-  1 ) 


and  must  satisfy  the  radiation  condition  as  z-±®  .  We  shall  choose  an 
C(z)  which  is  given  by 


,  z>  0 

e(z)  =  c(z)  ,  -L<z<0 

e!  .  z<-L 


(A- 2) 


where  l*(z)  is  an  arbitrary  function  of  z  .  The  formal  solution  to  Eq.  (A-l) 


is  well  known'  '  and  is  given  by 


$(z,p)  = 


(z<)  f(z>) 


-frit) 


(A- 3) 


where  z^  is  the  greater  of  z  or  z  ;  z<  is  the  lesser  of  z  or  z'  ;  $  (z) 
and^F(z)  are  independent  homogeneous  solutions  to  Eq.  (A-l)  which  satisfy 
the  boudnary  conditions  at  +  and  -  infinity  respectively;  and  the 
Wronskian  W  ,  ~T)  is  given  by 


<v>  =  f*(z)  4^  -  4^  • 


(A-4) 


i 
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At  this  point  we  will  limit  ourselves  to  finding  the  Green's  function 
for  z>0  and  z  >0  .  In  this  region  the  function  #  (z)  and  i(z)  are  given 
by 


i(z)  =  e 


lV2  Z 


z  >  0 


‘ik«P?Z  -  ikJ>7Z 

*  x  /  i  o  L  -p  O  L 

* (z)  =  e  +  i  e 


z  >0 


(A- 5) 

(A-6) 


.=£  2 


where P2=VC2-P  •  The  coefficient  T  will  be  determined  at  a  later  point. 
By  using  Eqs.  (A-5)  and  (A-6)  we  can  calculate  #(z,p)  .  It  is  given  by 

ik  p|z  -  z#|  __ 


"  -2ik 


ik  p  (z+z') 
o  c 


oP2 


2ik  n 
o  2 


(A- 7) 


All  that  remains  to  be  done  is  to  determine  I*  .  This  can  be  accom- 
plished  by  imposing  the  boundary  condition  on  $  (z)  that  it  must  be  out¬ 
going  at  -  00  .  This  gives 


1 


%)  = 


-ikp  z 

-  1  z 

ZT  O  C 

z  >  0 

e  + 

r  e 

af  (z)  + 

1  1 

8*2(z) 

-L<  z<0 

(A- 8) 

-ikoplz 

Y  e 

z  <  -  L 

where  $  (z)  and  $  (z)  are  two  independent  homogeneous  solutions  of 
1  2  / - 2 

Eq.  (A- 1 )  for  -L<z<0  and  p^-ye^-p  .  If  we  recall  that  the  one-dimen¬ 
sional  Green's  function,  $(z,  p)  is  related  to  the  electric  field  via  a  Fourier 
transform,  then  we  can  conclude  that  $(z,p)  and  (d/dz)$(z,p)  must  be  con¬ 
tinuous  across  the  discontinuities  of  e(z)  since  the  tangential  fields  E^  and 
H  are  continuous  across  these  discontinuities.  By  applying  these  conditions 
at  z  =  0,  we  obtain 


t 
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l  +  r  =  a«  j(o)  +  6  f2(o) 

-  ikp  +  ikp  I*  =  a{'(o)  +  9{'(o) 
o  L  o  Z  1  Z 


(A- 9) 

(A- !  0) 


By  imposing  the  same  conditions  at  z  =  -  L,  we  also  obtain 


ikp  L 

ye  =  a^t-L)  +  8i2(-L) 


ikoPlYe 


ikoPlL 


=  af'(-L)  +  9?'2(-L) 


(A- : ; ; 

(A- 12) 


where 


$!(z)  =  dt.(z)/dz  ,  i=l,2 


The  solution  of  these  four  equations  for  T  gives 


r  =-  f 

ab 


with 


(A- 13) 


t 

b 


j»'(-L)  +  ikjPl»2(.L)|  |*;<-u 


f'2(o)  ±  ikoP2#2(o) 


(o)  ±  ik  p  $  (o) 
i  O  Z  1 


It  must  be  emphasized  at  this  point  that  ^(z)  and  $ 2(z)  are  any  pair  of 
independent  solutions  to  the  homogeneous  equation.  For  each  different  pair 
of  solutions  we  obtain  a  different  representation  for  the  reflection  coefficient. 
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APPENDIX  B. 

INTEGRAND'S  SINGULARITIES 


The  pole  singularities  of  the  integrand  are  the  singularities  of  the 
reflection  coefficient,  (3.2.4)  ,  or  the  roots  of 


2 r ( v / 2  +  u/2  +  nr (y/2  -  u/2  +  1)  r_,  x  n_. 

'P2~'v/2  u/2  +1/2)  r(v/2-u/2  +  1  / 2 )  "L  (v+u)/2 J* < 1B > 


Since  the  square  root  (-:^-p^  )*  doesn’t  appear  in  the  above  equation,  a  root 
of  (IB)  will  correspond  to  a  pole  on  the  first  and  second  sheet  of  the  Riemann 
surface. 

The  poles  lying  on  the  integration  path,  C,  (Jigs. (1.5)  and  (1.6)  )  will  be 

* 

examined  first.  Consider  roots  of  (IB)  in  the  range  Rep^  =  0,  0  <Imp^<(i)  . 
The  right-hand  side  of  (IB)  is  real  and,  therefore,  no  roots  exist  in  this 
range.  In  the  range  0<  Re  p^  <  (-:  ^  )*  ,  Im  p^=  0,  ( IB )  can  be  simplified. 

Using  the  simplified  form 


2  T  (v/2+u/2+l) 
Tp2  ?2(v/2+u/2+l/2) 


-sinhrr  u  +  j  sin^v 
cos  TT  v  +  coshn  |u 


one  sees  the  right  hand  side  has  an  imaginary  part  when  sin(nv)^  0.  When 

sin(rv)  =  0,  the  right  hand  side  of  the  equation  is  less  than  zero.  Therefore, 

£ 

no  roots  exist  in  the  range  0<  Re  p^  <  (e^)  ,  Imp^=  0.  The  final  range  to 

* 

consider  is  Rep^=  0  ,  (i)  <Imp^<  00  .  Transform  ( 1 B )  to 

zrq/2  -  v/2  ■  lu|/2)  r  !»;;< 

'tIp2I.~(-v/2  +  |ul/2)r(v/2  +|u|/2+  1/2)  1  “ 

where  the  gamma  functions  appearing  in  (3B)  are  greater  than  zero  since 
lu|-v>0  for  the  range  of  p^  under  consideration.  As  a  result,  (3B)  has  no 
solutions  on  the  integration  path  C. 


I 
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Since  no  poles  lie  on  the  integration  path,  the  only  other  residue 

contributions  lie  between  the  original  and  deformed  path.  If  t  has  a  small  or 

moderate  value,  then  for  large  k  any  residue  contribution  will  be  a  rapidly 

o 

exponentially  decay  function  which  can  be  neglected  asymptotically.  When 
T  is  large  enough,  the  asymptotic  form  of  the  reflection  coefficient, 

R  -  e17"'*'*  +  Ipi>  1 

can  be  used.  Since  the  exponential  function  has  no  singularities  in  the  finite 
p^  plane,  the  reflection  coefficient  is  analytic  there.  As  a  result,  we  see 
no  residues  contribute  to  the  asymptotic  reflected  field  solution.  We  also 
note  that  since  the  saddle  points  lie  on  the  integration  path  C,  no  singularities 
approach  the  saddle  points  as  k  increases. 
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APPENDIX  C. 

DECAY  REGIONS 

The  decay  regions  of  the  integrand  appearing  in  (3.4.12)  are  those 


regions  where  Im  i (p ^ )  >  0  with 


-  -  - 1  £ 
*(Pj)  =  Vzz  +  px  +kQ  v(p1.  (A)  T)  • 


(1C) 


Since  the  integration  path  is  on  two  sheets  of  a  four-sheeted  Riemann  surface, 
the  decay  regions  must  be  found  for  both  of  these  sheets.  The  method  of 
finding  the  decay  regions  will  be:  first,  to  find  the  location  of  the  saddle  points 
and;  second,  to  find  the  constant  level  paths  through  these  saddle  points,  i.  e. 
paths  on  which  Im  * (p ^ )  =  0.  Since  the  constant  level  paths  separate  the  decay 
and  growth  regions,  once  these  paths  are  found,  the  boundaries  of  the  decay 
regions  are  known. 

The  location  of  the  saddle  point  can  be  found  in  the  main  part  of 

the  text.  On  the  top  sheet  of  the  Riemann  surface  three  saddle  points  were 

£ 

found.  T-.vo  were  on  the  positive  imaginary  axis  where  0  <  Im  p^<(A)  and 

£ 

one  was  on  the  negative  imaginary  axis  where  -(/')  <  Im  p^  <  0.  On  the 

second  sheet  there  was  only  one  saddle  point.  It  was  located  on  the  negative 

£ 

imaginary  axis  where  -(t'.t  <  Im  p^<  0. 

Two  constant  level  paths  pass  through  each  saddle  point  at  right 
angles  to  each  other.  The  imaginary  p^  axis  where  |Im  p^  I  <  (£)  is  one  of 
the  constant  level  paths.  The  other  paths  cross  the  imaginary  axis  at  the 
saddle  point  normal  to  the  axis.  Since  their  functional  behavior  is  complicated, 
only  approximate  path  locations  can  be  found.  The  asymptotes  of  these  paths 
as  |p^  |-°°  are  found  to  be 

.  ...*  .  . 

(2C) 


f(Pj.  (A)  T)  in  p  t  ,  p  »  1 


4 
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The  expansion  is  valid  on  both  the  upper  and  lower  sheet.  On  the  upper  sheet 
the  four  asymptotes  are: 


pu-[l 

;x-nL)/zJplr 

lSt  quadrant 

(3C; 

pli  =  - 

^(x-TTL)/zJ  plr 

_nd 

2  quadrant 

(4  C) 

pli  =  - 

[(x  +  rrL)/z!  p^ 

,rd 

3  quadrant 

(5C) 

1 

II 

•H 

a 

[(x  +  ttL)/zJ  p^ 

.th 

4  quadrant 

(6C) 

+  ipu 

.  An  examination  of  $  (p  )  for  large  k  shows  that  a 

1  o 

constant  level  path  goes  through  the  saddle  point  which  is  highest  on  the  imaginary 
p^  axis  and  the  path  is  asymptotic  to  (3 C)  and  (4C).  Similarly,  a  constant  lev-.l 
path  goes  through  the  saddle  point  which  is  lowest  on  the  imaginary  p^  axis 
and  the  path  is  asymptotic  to  (5C)  and  (6C).  There  are  no  asymptotes  on  the 
second  sheet.  As  a  result,  the  constant  level  path  always  remains  in  the  finite 
region  of  the  plane. 

There  is  one  constant  level  path  on  each  sheet  still  to  be  found.  An 

investigation  of  §(p^)  along  the  real  p^  axis  shows  that  Im  $>(pj)>  0  on  both 

*  i 

sheets  for  0  <  Re  <  (e^)  and  Im  §(p^)<0  on  both  sheets  for  (e  )  <  Re  p^<  0. 

£  . 

The  equation  for  Im  $(p^)  where  (e^)  <|Rep^  |  is 

ta[*<Pl>]  =  -  (Pj2-  CL)^  x  +  ko_1  Im  ^.(A^T)  (7 C) 

i 

which  is  valid  on  both  sheets.  This  equation  has  one  zero  for  Rep^  >  (e^) 

t 

and  another  for  Rep^<  -(e^)  when  t  small  or  moderate  and  k^large.  When 


T  becomes  large  (7C)  becomes 


Im  $(pL)  =  -  (pL  -  e  p  x  +  tt  p^L. 
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This  equation  has  a  similar  arrangement  of  zeros,  if  x  >  nL..  This  restriction 

does  not  hinder  us  since  the  caustic  lies  in  this  region.  With  the  above 

information,  the  trajectory  of  the  remaining  constant  level  paths  becomes 

clear.  The  path  forms  a  complete  circle  crossing  the  two  saddle  points,  one 

on  each  sheet,  at  right  angles  to  the  imaginary  axis  and  crossing  the 

£ 

Pi  =  ±  (£j)  branch  cuts  at  the  zeros  of  (7C).  The  constant  level  paths  and 
decay  regions  are  shown  in  Figs.  (3.3)  and  (3.4)  . 
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APPENDIX  D. 

EXTREMUM  OF  THE  LOCUS 
OF  TURNING  POINTS 


The  purpose  of  this  appendix  is  to  show  the  relationship  between 
the  locus  of  turning  points  and  the  caustic  formed  in  a  stratified  medium. 
Let  us  assume  that  there  exists  a  stratified  medium  whose  dielectric 
constant  varies  continuously  in  z  between  free  space  and  e^<l  .  Rays 
emitted  from  a  two-dimensional  point  source  located  at  (o,  z)  divide  into 
transmitted  and  returning  rays.  The  ray  equation  for  the  returning  rays 
after  they  have  turned  is  given  by 


* 

r  -  r 

J/  J 


'®(t)  “  p 


e(zt)  =  P 


(D-l) 


where  p  is  the  ray  parameter,  (x  ,  z  )  are  the  coordinates  of  turning  point 

v  t 

of  the  ray  and  (x,  z)  are  the  coordinates  of  ray.  The  caustic  formed  by  the 
rays  given  in  Eq.  (D-l)  can  be  obtained  by  solving  the  constraint  equation, 


(D-2) 


_d_  r  r  pfli- 

dz  )  J/  J  /  .  .  ’ 

t  I  L  z  ztJ  V  G  (t)-  p 


simultaneously  with  Eq.  (D-l)  .  The  constraint  equation  ,  D-2  ,  was  obtained 
by  taking  the  derivative  of  Eq.  (D-l)  with  respect  to  z^ .  We  see  Eq.  (D-2) 
will  give  p  or  z^  in  terms  of  z  which,  in  turn,  can  be  used  in  Eq.  (D-l)  to 
obtain  the  equation  for  the  caustic.  For  comparison  purposes  at  a  later 
time  we  will  evaluate  Eq.  (D-2)  at  z  =  z  .  The  result  is 


—  f 

dz^  J , 


(D-3) 


e(T)-p 
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The  locus  of  turning  points  for  the  turning  rays  is  given  by 


z 


pd  t 


G  (T)  -  p 


T 


(D-4) 


The  extremum  of  this  equation  is  given  by  dx^/dz^O  .  By  imposing 
this  condition  on  Eq.  (D-4) ,  the  resulting  equation  is  Eq.  (D-3)  .  This 
shows  that  extremum  of  the  locus  of  turning  points  correspond  to  points 
on  the  caustic  which  coincide  with  the  z  =  z'  line  . 


APPENDIX  E. 


The  above  transformation,  E-5  ,  is  a  relation  between  the  w-plane 
and  the  ?-plane  for  a  particular  value  of  argv  .  As  can  be  seen,  the 
transformation  is  independent  of  the  magnitude  of  v  .  Olver  has  made 
a  thorough  investigation  of  the  relationship  between  w  and  5  .  He  has 
shown  that  by  cutting  the  ?  plane  as  shown  in  Fig.  E-l  ,  the  asymptotic 
approximations  given  in  Eqs  .  (E-l)  -  (E -4)  are  valid  over  the  whole  w 
plane,  cut  as  shown  in  Fig,  E-2  with  |argv|<rr/2  . 


Fig.  E  -  1 
The  P  Plane 


Fig.  E  -  2 


The  w  Plane 
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The  branch  cuts  in  Fig.  E-l  are  constant  level  curves  of 
3/2  3  2/3 

exp(-2/3  ?  )  whose  level  is(  —  tt)  .  When  these  branch  cuts  are 

transformed  into  the  w  plane,  they  fall  on  top  of  one  another  and  form 
a  section  of  the  w  plane  branch  cut.  The  parameter  equations  of  this 
section  of  the  w  plane  branch  cut  are 


w  =  se ch  t  0<t<® 


o  -  tanh  n 


=  -  in  +  ,ei(3/2'-argv) 


(E-6) 


It  is  interesting  to  note  that  the  branch  cuts  change  position  as  the 
argv  is  varied,  however  this  will  not  concern  us  since  the  branch  cuts 
will  always  remain  out  of  the  region  of  our  interest. 

In  his  paper  Olver  has  obtained  the  asymptotic  approximations  to  the 
Bessel  functions  for  |argv|  <  tt/2  .  We  will  now  show  that  these  formulae 
are  valid  when  arg V  =- it/2  and  0<  w<°.  To  show  this  fact  we  will  find 
the  asymptotic  approximations  of  J^fvw)  in  the  region  |n  -argv|  <  tt/2  . 

If  the  asymptotic  approximations  from  the  right  and  left  hand  sides  of 
the  V  plane  yield  the  same  results  when  evaluated  at  argV  =  -  tt/2  ,  then 
Olver's  formula  is  correct  at  argv=  -  tt/2  .  If  the  two  formulae  differ  by 
more  than  an  exponentially  small  term  along  some  sector  of  the  Rew  axis, 
this  will  mean  that  the  particular  sector  is  an  anti-Stokes  line  'and  an 
oscillating  term  will  be  neglected  if  Olver's  original  formula  is  used  at 
argV  =  -  tt/2  .  We  should  not  make  the  faulty  conclusion,  however  that  if 
the  two  foimulae  agree,  this  sector  of  the  Rew  axis  is  not  an  anti-Stokes 
line. 

First,  evaluating  Olver's  asymptotic  formula,  E-l,  at  arg  v  =  -  rr/2  , 
w  real  and  0  < w  <  1  ,  we  have 
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J  (vw)  ~  e 
v 


in/2  Ml.lalrhw  Ai<|v|  I-' 


M|,l 


V  ,-w2 


(E-7) 


Next,  we  want  to  find  an  expression  for  J^(vw)  in  the  region  | tt  -  argv|<*r/2 

This  can  be  done  by  making  use  of  the  analytic  continuation  formula  given 
(42) 

by  Watson  .  The  result  is 


J  (vw)  =  J  (vwe*  )  =  e  V  1  J  —  (vw) 
v  -v  -v 


(£-8) 


where  v=  v  e*  and  then  using  the  definition  of  J  _(Vw)  ,  we  obtain 


V  TT  i  f  _  _  _  _  "| 

J^(vw)  =  e  Icos  vti  J-(vw)  -  sinvT  Y-  (vw)l  , 


(E-9) 


This  will  give  us  the  asymptotic  approximation  of  J^(vw)  in  terms  of 
functions  whose  asymptotic  approximations  are  valid  in  the  right  half  of 
the  v  plane.  The  asymptotic  approximation  to  Y—  (vw)  has  not  been 
given  as  yet  but  can  be  obtained  from  Eq.  (E-2)  since  Yv(z)  =  H^(z)-  H^(z)  . 
We  now  evaluate  Eq.  (E-9)  at  argv  =  -  tt/2  ,  0<  w<  1  and  neglect  exponentially 
small  terms.  The  result  is 


Jv(vw)  =  |  [j-(vw)  -  i  Y— ( Vw)  j 


(E- 10) 


Substituting  the  asymptotic  approximations  given  in  Eqs.(E-l)  and  (E-2), 
in  the  right  hand  side  of  the  above  equation,  gives  us 

‘-w  J  'v|  (E- 11) 


+  ein/3)]  . 


[ 
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Using  the  connection  formula  between  Airy  functions  shows  us  that 
the  asymptotic  approximation  for  J^(vw)  is  the  same  as  Eq.  (E-7)  . 

The  same  procedure  can  be  repeated  foT  arg  v  =  -  "12  ,  1<  w  <®  , 
which  shows  that  Olver’s  formula  can  be  extended  to  argv  -  -~/2  for 


this  section  of  the  w  plane  also.  In  addition  it  can  be  shown  that  the 

(1)  ,  (1)' 

asymptotic  formula  for  ^2  (vw)  ,  J  (vw)  and  r2  (vw)  are  also  valid 

Hv  V  H  v 

for  arg v= -tt/2  ,  0<w<®. 
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*,,T"*c^In  the  past  the  lateral  wave  has  been  investigated  for  the  case  of  a  sharplv 
(bounded  transition  layer  and  a  layer  with  a  linear  velocity  variation.  The  interpre’  - 
|tion  of  the  lateral  wave  given  in  these  instances  cannot  be  extended  to  an  arbitrar 
smooth  layer,  and  in  many  cases,  the  exact  nature  of  the  lateral  wave  contribution  be¬ 
comes  unclear.  It  is  our  purpose  to  clarify  these  matters  and  to  present  the  charac¬ 
teristics  of  lateral  wave  propagation  on  a  layer  of  arbitrary  variation.  The  models 
employed  can  simulate  an  inhomogeneous  plasma  having  a  number  densitv  profile,  n^, 
(which  varies  continuously  between  two  homogeneous  half  spaces.  An  integral  repre¬ 
sentation  for  the  scattered  field  in  the  optically  denser  half  space  is  found  and  evaluai- 
ed  asymptotically  in  the  high  frequency  limit.  This  asymptotic  evaluation  is  carried 
out  in  two  parameter  ranges:  first,  when  the  layer  is  thick  compared  with  wavelength; 
and  second,  when  the  layer  thickness  is  arbitrary  but  the  observation  point's  distance 
along  the  interface  is  large  compared  with  layer  thickness. 

When  the  layer  thickness  is  large  compared  with  wavelength,  the  asymptotic  analy¬ 
sis  of  the  scattered  field  shows  that  the  interpretation  of  the  lateral  wave  depends 
markedly  upon  the  gradient  of  n(z)  at  the  junction  with  the  optically  rarer  homogeneous 
palf  space.  It  is  found  that  when  a  finite  gradient  of  n(z)  exists,  the  conventional  inter 
pretation  of  the  lateral  wave  contribution  is  correct;  however,  the  lateral  wave  mech¬ 
anism  is  different  in  the  case  of  a  zero  gradient.  For  observation  points  situated  at  a 
Jlarge  distance  .along  a  layer  of  arbitrary  thickness,  the  asymptotic  expression  for  the 
ateral  wave  contribution  has  an  amplitude  dependence  on  distance  identical  with  that 

Eot  the  lateral  wave  on  an  abrupt  interface.  In  addition,  the  lateral  wave  expression 
•educes  to  the  thick  layer  result  for  large  layer  thickness  and  it  reduces  to  the  abrupt 
nterface  result  for  thin  layer  thickness. 
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